Health impacts of waterborne microbes: Difference between revisions

Latest revision as of 10:10, 11 September 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.store(Wateruse)
cat("Ovariable Wateruse saved. \n")

+ 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)
cat("Ovariable RawConsumption saved. \n")

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

+ 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", "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)
                                }


)

objects.store(MicrobeLogDecrease)
cat("Ovariable MicrobeLogDecrease saved. \n")

Pathogen concentration in drinking water

+ Show code - Hide code

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

TreatedConcentration <- Ovariable("TreatedConcentration", dependencies = data.frame(
                                  Name=c("RawConcentration", "MicrobeLogDecrease"),
                                  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)
                                  })

objects.store(TreatedConcentration)

cat("Ovariable TreatedConcentration saved. \n")

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(...){
                        out <- TreatedConcentration * RawConsumption
                        colnames(out@output)[colnames(out@output) == "Pathogen"] <- "Exposure_agent"
                        return(out)
                      }
)

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

Health impacts of waterborne microbes: Difference between revisions

Latest revision as of 10:10, 11 September 2019

Contents

Question

Answer

Rationale

Data

Calculations

RCode

See also

Navigation menu