Health impacts of waterborne microbes: Difference between revisions

Revision as of 06:09, 16 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/logdecrease" on page [[Health impacts of waterborne microbes]]
library(OpasnetUtils)

	
MicrobeLogDecrease <- Ovariable("MicrobeLogDecrease",
                                dependencies=data.frame(
                                  Name=c("TreatmentEfficiency", "Disinfection", "ChlorineEfficiency"),
                                  Ident=c("Op_en7954/variable", "Op_en7955/variable", "Op_en7956/efficiency")
                                ),
                                formula=function(...){
                                  # Calculate the total log decrease of pathogen concentration from disinfection and other treatment
                                  # methods, then sum them to get the total not including chlorination
                                  MicrobeLogDecrease <- oapply(Disinfection, cols="TreatmentMethod", FUN=sum) +
                                    oapply(TreatmentEfficiency, cols="TreatmentMethod", FUN=sum)
                                  # Then sum that with chlorination to get the total microbe log decrease
                                  MicrobeLogDecrease <- MicrobeLogDecrease+ChlorineEfficiency
                                  
                                  return(MicrobeLogDecrease)
                                }


)

objects.store(MicrobeLogDecrease)

+ Show code - Hide code

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



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)

@@ Line 38: / Line 38: @@
 oprint(Wateruse)
 </rcode>
 <rcode name="rawconsumption">
@@ Line 57: / Line 56: @@
 Concentrations of pathogens in drinking water
+<rcode
+	name="logdecrease"
+	label="Initialize microbe log decrease"
+	graphics="1"
+>
+#This is code "Op_en7957/logdecrease" on page [[Health impacts of waterborne microbes]]
+library(OpasnetUtils)
+MicrobeLogDecrease <- Ovariable("MicrobeLogDecrease",
+                                dependencies=data.frame(
+                                  Name=c("TreatmentEfficiency", "Disinfection", "ChlorineEfficiency"),
+                                  Ident=c("Op_en7954/variable", "Op_en7955/variable", "Op_en7956/efficiency")
+                                ),
+                                formula=function(...){
+                                  # Calculate the total log decrease of pathogen concentration from disinfection and other treatment
+                                  # methods, then sum them to get the total not including chlorination
+                                  MicrobeLogDecrease <- oapply(Disinfection, cols="TreatmentMethod", FUN=sum) +
+                                    oapply(TreatmentEfficiency, cols="TreatmentMethod", FUN=sum)
+                                  # Then sum that with chlorination to get the total microbe log decrease
+                                  MicrobeLogDecrease <- MicrobeLogDecrease+ChlorineEfficiency
+                                  return(MicrobeLogDecrease)
+                                }
+)
+objects.store(MicrobeLogDecrease)
+</rcode>
 <rcode
 	name="pathconcentration"
-	label="Initialize variable"
+	label="Initialize microbe log decrease"
 	graphics="1"
 >
@@ Line 66: / Line 98: @@
 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)

Health impacts of waterborne microbes: Difference between revisions

Revision as of 06:09, 16 July 2019

Contents

Question

Answer

Rationale

Data

Calculations

RCode

See also

Navigation menu