Health impacts of waterborne microbes: Difference between revisions

From Opasnet
Jump to navigation Jump to search
mNo edit summary
 
(17 intermediate revisions by the same user not shown)
Line 28: Line 28:
Daily cold tap water use of Finns  
Daily cold tap water use of Finns  


<rcode name="consumption">
<rcode
      name="consumption"
      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]]
library(OpasnetUtils)
library(OpasnetUtils)


Wateruse = Ovariable("Wateruse", data = data.frame(WateruseResult = 1153), save = TRUE)
Wateruse = Ovariable("Wateruse", data = data.frame(WateruseResult = 1153))
 
objects.store(Wateruse)
cat("Ovariable Wateruse saved. \n")


oprint(Wateruse)
</rcode>
</rcode>


 
<rcode
<rcode name="rawconsumption">
        name="rawconsumption"
        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]]
library(OpasnetUtils)
library(OpasnetUtils)


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


function = function(...){
objects.store(RawConsumption)
return(2 / 3 * Wateruse / 1000) # 2/3 consumed unboiled. Converted into liters.
cat("Ovariable RawConsumption saved. \n")
}
 
</rcode>
 
Log decrease of concentration of microbes in drinking water due to water treatment
 
<rcode
name="logdecrease"
label="Initialize microbe log decrease"
graphics="1"
embed=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", "ChlorineEfficiency", "OzoneEfficiency", "UVEfficiency"),
                                  Ident=c("Op_en7954/variable", "Op_en7956/efficiency", "Op_en7955/ozefficiency", "Op_en7955/uvefficiency")
                                ),
                                formula=function(...){
                                  # Calculate the total log decrease of pathogen concentration from treatment methods
                                  MicrobeLogDecrease <- oapply(TreatmentEfficiency, cols="TreatmentMethod", FUN=sum)
                                  # Then sum that with the disinfection methods to get the total microbe log decrease
                                  MicrobeLogDecrease <- MicrobeLogDecrease+ChlorineEfficiency+OzoneEfficiency+UVEfficiency
                                 
                                  return(MicrobeLogDecrease)
                                }
 


RawConsumption <- Ovariable("RawConsumption", dependencies = dependencies, formula = function, save = TRUE)
)


oprint(RawConsumption)
objects.store(MicrobeLogDecrease)
</rcode>
cat("Ovariable MicrobeLogDecrease saved. \n")


Concentrations of pathogens in drinking water
</rcode>
Pathogen concentration in drinking water


<rcode  
<rcode  
name="pathconcentration"  
name="pathconcentration"  
label="Initialize variable"  
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]]
library(OpasnetUtils)
library(OpasnetUtils)


dependencies <- data.frame(
TreatedConcentration <- Ovariable("TreatedConcentration", dependencies = data.frame(
Name = c("RawConcentration", "Disinfection", "TreatmentEfficiency", "ChlorineEfficiency"),
                                  Name=c("RawConcentration", "MicrobeLogDecrease"),
Ident = c("Op_en7953/variable", "Op_en7955/variable", "Op_en7954/variable", "Op_en7956/efficiency")
                                  Ident=c("Op_en7953/variable", "Op_en7957/logdecrease")
)
                                  ),
                                  formula = function(...){
 
                                    # Calculate the concentration of bacteria after treatment from log decrease
                                    out <- RawConcentration * 10 ^(-1*(MicrobeLogDecrease))
                                   
                                    return(out)
                                  })


funktio <- function(...) {
objects.store(TreatedConcentration)
#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, save = TRUE)
cat("Ovariable TreatedConcentration saved. \n")


oprint(TreatedConcentration)
</rcode>
</rcode>


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


<rcode name="exposure">
<rcode
        name="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)


dependencies = data.frame(
exposure <- Ovariable("exposure",
Name = c("TreatedConcentration", "RawConsumption"),
                      dependencies = data.frame(
Ident = c("Op_en7957/pathconcentration", "Op_en7957/rawconsumption")
                        Name = c("TreatedConcentration", "RawConsumption"),
                        Ident = c("Op_en7957/pathconcentration", "Op_en7957/rawconsumption")
                        ),
                      formula = function(...){
                        out <- TreatedConcentration * RawConsumption
                        colnames(out@output)[colnames(out@output) == "Pathogen"] <- "Exposure_agent"
                        return(out)
                      }
)
)


funktio = function(...){
objects.store(exposure)
return(TreatedConcentration * RawConcumption)
cat("Ovariable exposure saved. \n")
}


Exposure <- Ovariable("Exposure", dependencies = dependencies, formula = funktio, save = TRUE)
oprint(Exposure)
</rcode>
</rcode>


== See also ==
== See also ==
Line 148: Line 158:
*http://permanent.access.gpo.gov/lps35390/cfpub.epa.gov/ncea/cfm/recordisplay.cfm-deid=55145.htm
*http://permanent.access.gpo.gov/lps35390/cfpub.epa.gov/ncea/cfm/recordisplay.cfm-deid=55145.htm
*Finley, B., Proctor, D., Scott, P., Harrington, N., Paustenbach, D., Price, P. 1994. [http://www.ncbi.nlm.nih.gov/pubmed/7972957 Recommended distributions for exposure factors frequently used in health risk assessment.]  Risk Analysis 14(4), 533-553.
*Finley, B., Proctor, D., Scott, P., Harrington, N., Paustenbach, D., Price, P. 1994. [http://www.ncbi.nlm.nih.gov/pubmed/7972957 Recommended distributions for exposure factors frequently used in health risk assessment.]  Risk Analysis 14(4), 533-553.
*Havelaar, A. and Melse, J. M. (2003). Quantifying public health risk in the WHO Guidelines for Drinking‐ Water Quality: a burden of disease approach. RIVM raport 734301022. Bilthoven, RIVM.

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