Health impacts of waterborne microbes: Difference between revisions

From Opasnet
Jump to navigation Jump to search
 
(4 intermediate revisions by the same user not shown)
Line 60: 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"
embed=1
>
# 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 94: 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 128: 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 157: Line 136:
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 163: 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

Retrieved from "https://en.opasnet.org/index.php?title=Health_impacts_of_waterborne_microbes&oldid=43381"