Health impacts of waterborne microbes: Difference between revisions

Revision as of 10:44, 15 July 2019

This page is a knowledge crystal of subtype variable. The page identifier is Op_en7957
Moderator:Heta (see all)

Upload data {{#opasnet_base_link:Op_en7957}}

Question

What are the health impacts of waterborne microbes in drinking water?

Answer

Rationale

Data

Health impacts of waterborne microbes depends on the following parameters:

Exposure to pathogens
- Drinking water consumption
- Concentrations of pathogens in drinking water

Calculations

RCode

Daily cold tap water use of Finns

+ Show code - Hide code

# This is code Op_en7957/consumption on page [[Health impacts of waterborne microbes]]
library(OpasnetUtils)

Wateruse = Ovariable("Wateruse", data = data.frame(WateruseResult = 1153))

objects.save(Wateruse)

oprint(Wateruse)

+ Show code - Hide code

# This is code Op_en7957/rawconsumption in page [[Health impacts of waterborne microbes]]
library(OpasnetUtils)

RawConsumption <- Ovariable("RawConsumption",
                            dependencies = data.frame(Name = c("Wateruse"), Ident = c("Op_en7957/consumption")),
                            formula = function(...){
                              return(2 / 3 * Wateruse / 1000)  # 2/3 consumed unboiled. Converted into liters.
                            }
)

objects.store(RawConsumption)

oprint(RawConsumption)

Concentrations of pathogens in drinking water

+ Show code - Hide code

#This is code "Op_en7957/pathconcentration" on page [[Health impacts of waterborne microbes]]
library(OpasnetUtils)

dependencies <- data.frame(
		Name = c("RawConcentration", "Disinfection", "TreatmentEfficiency", "ChlorineEfficiency"),
		Ident = c("Op_en7953/variable", "Op_en7955/variable", "Op_en7954/variable", "Op_en7956/efficiency")
)

funktio <- function(...) {
	
	#MicrobeLogDecrease <- combine(Disinfection, TreatmentEfficiency, ChlorineEfficiency, name = "MicrobeLogDecrease")
	
	#MicrobeLogDecrease@output[["TreatmentMethod"]][ 
	#		MicrobeLogDecrease@output[["MicrobeLogDecreaseSource"]] == "ChlorineEfficiencyF"
	#] <- "Chlorination"
	
	#MicrobeLogDecrease <- oapply(MicrobeLogDecrease, cols = c("TreatmentMethod", "ChlorineDoseSource"), FUN = sum, na.rm = TRUE)
	# Varo laskemasta kloorausta monesti
	
	#MicrobeLogDecrease@output <- fillna(MicrobeLogDecrease@output, colnames(MicrobeLogDecrease@output)[MicrobeLogDecrease@marginal])
	
	MicrobeLogDecrease <- orbind(Disinfection, TreatmentEfficiency)
	MicrobeLogDecrease <- orbind(MicrobeLogDecrease, ChlorineEfficiency)
	
	MicrobeLogDecrease <- fillna(MicrobeLogDecrease, grep("Source$", colnames(MicrobeLogDecrease)))
	
	MicrobeLogDecrease$DisinfectionResult[is.na(MicrobeLogDecrease$DisinfectionResult)] <- 0
	MicrobeLogDecrease$TreatmentEfficiencyResult[is.na(MicrobeLogDecrease$TreatmentEfficiencyResult)] <- 0
	MicrobeLogDecrease$ChlorineEfficiencyResult[is.na(MicrobeLogDecrease$ChlorineEfficiencyResult)] <- 0
	
	MicrobeLogDecrease$MicrobeLogDecreaseResult <- MicrobeLogDecrease$DisinfectionResult + 
			MicrobeLogDecrease$TreatmentEfficiencyResult + MicrobeLogDecrease$ChlorineEfficiencyResult
	
	MicrobeLogDecrease <- Ovariable(
		"MicrobeLogDecrease", 
		dependencies = data.frame(
			Name = c("Disinfection", "TreatmentEfficiency", "ChlorineEfficiency")
		),
		output = MicrobeLogDecrease
	)
	
	MicrobeLogDecrease <- CheckMarginals(MicrobeLogDecrease, verbose = FALSE)
	
	MicrobeLogDecrease <- oapply(MicrobeLogDecrease, cols = c("TreatmentMethod"), FUN = sum, na.rm = TRUE)
	
	
	
	
	
	out <- RawConcentration * 10 ^(-1*(MicrobeLogDecrease))
	
	return(out@output)
}

TreatedConcentration <- Ovariable("TreatedConcentration", dependencies = dependencies, formula = funktio)

objects.store(TreatedConcentration)

The amount of pathogens an average water consumer is exposed to in a day

+ Show code - Hide code

#This is code "Op_en7957/exposure" on page [[Health impacts of waterborne microbes]]
library(OpasnetUtils)

Exposure <- Ovariable("Exposure", 
                      dependencies = data.frame(
                        Name = c("TreatedConcentration", "RawConsumption"),
                        Ident = c("Op_en7957/pathconcentration", "Op_en7957/rawconsumption")
                        ),
                      formula = function(...){
                        return(TreatedConcentration * RawConcumption)
                      }
)

objects.store(Exposure)

Health impacts of waterborne microbes: Difference between revisions

Revision as of 10:44, 15 July 2019

Contents

Question

Answer

Rationale

Data

Calculations

RCode

See also

Navigation menu