ERF of waterborne microbes

This page is a knowledge crystal of subtype variable. The page identifier is Op_en7948
Moderator:Heta (see all)
Upload data {{#opasnet_base_link:Op_en7948}}

Question

What are the dose-response functions of pathogens in drinking water?

Answer

+ Show code- Hide code library(OpasnetUtils) objects.latest("Op_en7948", code_name = "variable") oprint(AnnosVaste)

Rationale

Data

ERF(-)

Obs	Exposure agent	Response	Exposure	Exposure unit	ER function	Scaling	Threshold	ERF	Description
1	campylobacter	campylobacter infection	ingestion	?	beta poisson approximation	None	0.011	0.024
2	rotavirus	rotavirus infection	ingestion	?	exact beta poisson	None	0.191	0.167
3	norovirus	norovirus infection	ingestion	?	exact beta poisson	None	0.055	0.04
4	cryptosporidium	cryptosporidium infection	ingestion	?	exact beta poisson	None	0.176	0.115
5	giardia	giardia infection	ingestion	?	exponential	None	0	0.0199
6	E.coli O157:H7	E. coli infection	ingestion	?	exact beta poisson	None	9.16	0.157	E. coli O157:H7 strain

Probability of illness, given infection(probability)

Obs	Response	Illness	Subgroup	Result
1	giardia infection	gastroenteritis		1
2	cryptosporidium infection	gastroenteritis		0.71
3	cryptosporidium infection	death	Age:age 0-4	0.00000070645
4	cryptosporidium infection	death	Age:age 5-9	0
5	cryptosporidium infection	death	Age:age 10-14	0
6	cryptosporidium infection	death	Age:age 15-64	0.000014839
7	cryptosporidium infection	death	Age:age 65-79	0.0011644
8	cryptosporidium infection	death	Age:age 80+	0.0011644
9	norovirus infection	gastroenteritis		0.7
10	norovirus infection	death	Age:age 0-4	0.000002058
11	norovirus infection	death	Age:age 5-9	0
12	norovirus infection	death	Age:age 10-14	0
13	norovirus infection	death	Age:age 15-64	0
14	norovirus infection	death	Age:age 65-79	0.000168
15	norovirus infection	death	Age:age 80+	0.000168
16	rotavirus infection	gastroenteritis		0.9
17	rotavirus infection	death	Age:age 0-4	0.00001917
18	rotavirus infection	death	Age:age 5-9	0
19	rotavirus infection	death	Age:age 10-14	0
20	rotavirus infection	death	Age:age 15-64	0
21	rotavirus infection	death	Age:age 65-81	0
22	rotavirus infection	death	Age:age 80+	0
23	campylobacter infection	gastroenteritis		0.323
24	campylobacter infection	clinical Guillian-Barré syndrome		0.000066
25	campylobacter infection	residual Guillian Barré syndrome		0.000066
26	campylobacter infection	reactive arthritis		0.0066
27	campylobacter infection	death		0.000131
28	E. coli infection	gastroenteritis		0.4346
29	E. coli infection	hemorrhagic colitis		0.3854
30	E. coli infection	haemolytic uraemic syndrome		0.0082
31	E. coli infection	end stage renal disease		0.0009676
32	E. coli infection	death		0.00062238

⇤--arg6688: . Find out the units of the parameters to understand the functions precisely. --Jouni (talk) 12:54, 11 July 2019 (UTC) (type: truth; paradigms: science: attack)

----arg2087: . The ERF and Threshold numbers in the tables correspond to parameters called alpha and beta respectively in the papers cited below in the cases of campylo, rota, noro and crypto. The equations alpha and beta are used in in the papers don't (seem to) match the equations below. At least for campylo the equations are significantly more complicated. For giardia the value in ERF does seem to be used in the paper it's from in the way the exponential equation shows. In the paper it's marked as r, the probability of infection. The equation in the paper uses "mean number per portion" where Dose is in this equation. E.coli is the biggest mystery: I wasn't able to find the numbers in the table above from the paper sited for it at all. --Heta (talk) 10:32, 12 July 2019 (UTC) (type: truth; paradigms: science: comment)

In these equations, Param1 and Param2 are in Threshold and ERF columns, respectively.

• Beta Poisson approximation: 1-(1+Dose/Param2)^-Param1
• Exact beta Poisson: 1-exp(-(Param1/(Param1+Param2))*Dose)
• Exponential: 1-exp(-Param1*Dose)

Pathogen	Reference
Campylobacter
Rotavirus
Norovirus
Cryptosporidium
Giardia
E.coli O157:H7

Calculations

+ Show code- Hide code # This is code Op_en7948/ERF_micr2 on page [[ERF of waterborne microbes]] # Note! This version has ERF and threshold in the same ovariable. library(OpasnetUtils) ERF_micr <- Ovariable("ERF_micr", ddata = "Op_en7948", subset="ERF") colnames(ERF_micr@data) <- gsub(" ", "_", colnames(ERF_micr@data)) objects.store(ERF_micr) cat("Ovariable ERF_micr stored.\n")

+ Show code- Hide code # This is code Op_en7948/ERF_micr on page [[ERF of waterborne microbes]] library(OpasnetUtils) d <- opbase.data("Op_en7948")[-1] colnames(d) <- gsub(" ", "_", colnames(d)) d$Result <- ifelse(d$Result == "", "0", as.character(d$Result)) ERF_micr <- Ovariable("ERF_micr", data = d[d$Observation == "ERF", colnames(d) != "Observation"]) threshold_micr <- Ovariable("threshold_micr", data = d[d$Observation == "Threshold", colnames(d) != "Observation"]) objects.store(ERF_micr, threshold_micr) cat("Ovariables ERF_micr, threshold_micr stored.\n")

⇤--arg5268: . The code below is old and does not work with the new ERF table. --Jouni (talk) 12:54, 11 July 2019 (UTC) (type: truth; paradigms: science: attack)

+ Show code- Hide code library(OpasnetUtils) f.exact.beta.poisson <- function(Dose, Param1, Param2, ...) 1 - exp(-(Param1 / (Param1 + Param2)) * Dose) f.beta.poisson.approx <- function(Dose, Param1, Param2, ...) 1 - (1 + Dose / Param2)^(-Param1) f.exponential <- function(Dose, Param1, ...) 1 - exp(-Param1 * Dose) f.alt.exponential <- function(Dose, Param1, Param2, ...) 1 - (1 + (Dose / Param1))^(Param2) EFR <- Ovariable("ERF", data = tidy(opbase.data("Op_fi2653"))) objects.store(ERF, f.exact.beta.poisson, f.beta.poisson.approx, f.exponential, f.alt.exponential) cat("Ovariable ERF and functions f.exact.beta.poisson, f.beta.poisson.approx, f.exponential and f.alt.exponential saved.\n") oprint(ERF)

See also

• Water guide
• HIA
• http://wiki.camra.msu.edu/index.php?title=Dose_response_assessment

References