Health impacts of waterborne microbes: Difference between revisions

From Opasnet
Jump to navigation Jump to search
 
(6 intermediate revisions by the same user not shown)
Line 31: Line 31:
       name="consumption"
       name="consumption"
       label="Initialize water use"
       label="Initialize water use"
      embed=1
>
>
# This is code Op_en7957/consumption on page [[Health impacts of waterborne microbes]]
# This is code Op_en7957/consumption on page [[Health impacts of waterborne microbes]]
Line 45: Line 46:
         name="rawconsumption"
         name="rawconsumption"
         label="Initialize raw consumption"
         label="Initialize raw consumption"
embed=1
>
>
# This is code Op_en7957/rawconsumption in page [[Health impacts of waterborne microbes]]
# This is code Op_en7957/rawconsumption in page [[Health impacts of waterborne microbes]]
Line 58: Line 60:
objects.store(RawConsumption)
objects.store(RawConsumption)
cat("Ovariable RawConsumption saved. \n")
cat("Ovariable RawConsumption saved. \n")
</rcode>
The raw water default class: Surface water - high contamination
<rcode
name="rawclass"
label="Initialize raw water class"
>
# This is code "Op_en 7957/rawclass" on page [[Health impact of waterborne microbes]]
library(OpasnetUtils)
RawClass <- Ovariable("RawClass", data=data.frame(RawWaterClass = "Surface water - high contamination", RawClassResult = 1))
objects.store(RawClass)
cat("Ovariable RawClass saved. \n")


</rcode>
</rcode>
Line 83: Line 69:
label="Initialize microbe log decrease"  
label="Initialize microbe log decrease"  
graphics="1"
graphics="1"
embed=1
>
>
#This is code "Op_en7957/logdecrease" on page [[Health impacts of waterborne microbes]]
#This is code "Op_en7957/logdecrease" on page [[Health impacts of waterborne microbes]]
Line 90: Line 77:
MicrobeLogDecrease <- Ovariable("MicrobeLogDecrease",
MicrobeLogDecrease <- Ovariable("MicrobeLogDecrease",
                                 dependencies=data.frame(
                                 dependencies=data.frame(
                                   Name=c("TreatmentEfficiency", "Disinfection", "ChlorineEfficiency"),
                                   Name=c("TreatmentEfficiency", "ChlorineEfficiency", "OzoneEfficiency", "UVEfficiency"),
                                   Ident=c("Op_en7954/variable", "Op_en7955/variable", "Op_en7956/efficiency")
                                   Ident=c("Op_en7954/variable", "Op_en7956/efficiency", "Op_en7955/ozefficiency", "Op_en7955/uvefficiency")
                                 ),
                                 ),
                                 formula=function(...){
                                 formula=function(...){
                                   # Calculate the total log decrease of pathogen concentration from disinfection and other treatment
                                   # Calculate the total log decrease of pathogen concentration from treatment methods
                                  # methods, then sum them to get the total not including chlorination
                                   MicrobeLogDecrease <- oapply(TreatmentEfficiency, cols="TreatmentMethod", FUN=sum)
                                   MicrobeLogDecrease <- oapply(Disinfection, cols="TreatmentMethod", FUN=sum) +
                                   # Then sum that with the disinfection methods to get the total microbe log decrease
                                    oapply(TreatmentEfficiency, cols="TreatmentMethod", FUN=sum)
                                   MicrobeLogDecrease <- MicrobeLogDecrease+ChlorineEfficiency+OzoneEfficiency+UVEfficiency
                                   # Then sum that with chlorination to get the total microbe log decrease
                                   MicrobeLogDecrease <- MicrobeLogDecrease+ChlorineEfficiency
                                    
                                    
                                   return(MicrobeLogDecrease)
                                   return(MicrobeLogDecrease)
Line 118: Line 103:
label="Initialize pathogen concentration"  
label="Initialize pathogen concentration"  
graphics="1"
graphics="1"
embed=1
>
>
#This is code "Op_en7957/pathconcentration" on page [[Health impacts of waterborne microbes]]
#This is code "Op_en7957/pathconcentration" on page [[Health impacts of waterborne microbes]]
Line 123: Line 109:


TreatedConcentration <- Ovariable("TreatedConcentration", dependencies = data.frame(
TreatedConcentration <- Ovariable("TreatedConcentration", dependencies = data.frame(
                                   Name=c("RawConcentration", "MicrobeLogDecrease", "RawClass"),
                                   Name=c("RawConcentration", "MicrobeLogDecrease"),
                                   Ident=c("Op_en7953/variable", "Op_en7957/logdecrease", "Op_en7957/rawclass")
                                   Ident=c("Op_en7953/variable", "Op_en7957/logdecrease")
                                   ),
                                   ),
                                   formula = function(...){
                                   formula = function(...){
                                    # Use the rows of RawConcentration with the right type of raw water
 
                                    RawConcentration <- RawConcentration*RawClass
                                   
                                     # Calculate the concentration of bacteria after treatment from log decrease
                                     # Calculate the concentration of bacteria after treatment from log decrease
                                     out <- RawConcentration * 10 ^(-1*(MicrobeLogDecrease))
                                     out <- RawConcentration * 10 ^(-1*(MicrobeLogDecrease))
Line 147: Line 131:
         name="exposure"
         name="exposure"
         label="Initialize exposure"
         label="Initialize exposure"
embed=1
>
>
#This is code "Op_en7957/exposure" on page [[Health impacts of waterborne microbes]]
#This is code "Op_en7957/exposure" on page [[Health impacts of waterborne microbes]]
library(OpasnetUtils)
library(OpasnetUtils)


Exposure <- Ovariable("Exposure",  
exposure <- Ovariable("exposure",  
                       dependencies = data.frame(
                       dependencies = data.frame(
                         Name = c("TreatedConcentration", "RawConsumption"),
                         Name = c("TreatedConcentration", "RawConsumption"),
Line 157: Line 142:
                         ),
                         ),
                       formula = function(...){
                       formula = function(...){
                         return(TreatedConcentration * RawConsumption)
                         out <- TreatedConcentration * RawConsumption
                        colnames(out@output)[colnames(out@output) == "Pathogen"] <- "Exposure_agent"
                        return(out)
                       }
                       }
)
)


objects.store(Exposure)
objects.store(exposure)
cat("Ovariable Exposure saved. \n")
cat("Ovariable exposure saved. \n")


</rcode>
</rcode>

Latest revision as of 10:10, 11 September 2019



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:

Calculations

RCode

Daily cold tap water use of Finns

+ Show code

+ Show code

Log decrease of concentration of microbes in drinking water due to water treatment

+ Show code

Pathogen concentration in drinking water

+ Show code

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

+ Show code

See also