Väestön kohdekohtainen ympäristöperäisen haitta-ainealtistumisen ja terveysriskin arviointi

Tämä sivu on tietokide alatyypiltään metodi. Sivutunniste: Op_fi2814
Moderaattori:Ei ole (katso kaikki) Kuinka ryhtyä moderaattoriksi? Sivun edistymistä ei ole arvioitu. Arvostuksen määrää ei ole arvioitu (ks. peer review).
Lisää dataa {{#opasnet_base_link:Op_fi2814}}

Kysymys

Kuinka arvioidaan kohdekohtaisesti väestön ympäristöperäinen altistuminen haitta-aineelle sekä altistumisesta aiheutuva terveysriski?

Vastaus

Ennakolta määritellyt laskentaparametrit

Ero sivun ”Väestön kohdekohtainen ympäristöperäisen haitta-ainealtistumisen ja terveysriskin arviointi” versioiden välillä(useita)

Obs	Parametri	Lyhenne	Yksikkö	Henkilö	Arvo
1	Altistumisperiodin pituus	At	vrk	Aikuinen	365
2	Altistumisperiodin pituus	At	vrk	Lapsi	365
3	Kehon paino	Bw	kg	Aikuinen	70
4	Kehon paino	Bw	kg	Lapsi	15
5	Hengitystiheys	Ir.air	m3/h	Aikuinen	0.83
6	Hengitystiheys	Ir.air	m3/h	Lapsi	0.32
7	Oleskeluaika sisätiloissa (keskiarvo altistumsperiodille)	Ef.air.indoor	h/vrk	Aikuinen	22.86
8	Oleskeluaika sisätiloissa (keskiarvo altistumsperiodille)	Ef.air.indoor	h/vrk	Lapsi	21.14
9	Juomaveden kulutus	Ir.dw	l/vrk	Aikuinen	0.7
10	Juomaveden kulutus	Ir.dw	l/vrk	Lapsi	0.2
11	Altistumistiheys, juomavesi	Ef.dw	vrk/vuosi	Aikuinen	365
12	Altistumistiheys, juomavesi	Ef.dw	vrk/vuosi	Lapsi	365
13	Niellyn maa-aineksen määrä	Ir.soil	mg/vrk/ton	Aikuinen	50
14	Niellyn maa-aineksen määrä	Ir.soil	mg/vrk	Lapsi	150
15	Maa-aineksen kiinnittyminen iholle	Soil.ad.skin	mg/cm2/vrk	Aikuinen	0.25
16	Maa-aineksen kiinnittyminen iholle	Soil.ad.skin	mg/cm2/vrk	Lapsi	0.25
17	Maa-ainekselle altistuvan ihon osuus	F.skin.soil		Aikuinen	0.05
18	Maa-ainekselle altistuvan ihon osuus	F.skin.soil		Lapsi	0.2
19	Altistumistiheys, maa-aines	Ef.soil	vrk/vuosi	Aikuinen	255
20	Altistumistiheys, maa-aines	Ef.soil	vrk/vuosi	Lapsi	255
21	Vihannesten kokonaiskulutus (keskiarvo altistumsperiodille)	Ir.lfc	kg/vrk	Aikuinen	0.110
22	Vihannesten kokonaiskulutus (keskiarvo altistumsperiodille)	Ir.lfc	kg/vrk	Lapsi	0.055
23	Juuresten kokonaiskulutus (keskiarvo altistumsperiodille)	Ir.tfc	kg/vrk	Aikuinen	0.200
24	Juuresten kokonaiskulutus (keskiarvo altistumsperiodille)	Ir.tfc	kg/vrk	Lapsi	0.100
25	Marjojen kokonaiskulutus (keskiarvo altistumsperiodille)	Ir.berries	kg/vrk	Aikuinen	0.044
26	Marjojen kokonaiskulutus (keskiarvo altistumsperiodille)	Ir.berries	kg/vrk	Lapsi	0.035
27	Paikallisten vihannesten osuus kokonaiskulutuksesta (keskiarvo altistumsperiodille)	F.lfc.loc		Aikuinen	0.1
28	Paikallisten vihannesten osuus kokonaiskulutuksesta (keskiarvo altistumsperiodille)	F.lfc.loc		Lapsi	0.1
29	Paikallisten juuresten osuus kokonaiskulutuksesta (keskiarvo altistumsperiodille)	F.tfc.loc		Aikuinen	0.1
30	Paikallisten juuresten osuus kokonaiskulutuksesta (keskiarvo altistumsperiodille)	F.tfc.loc		Lapsi	0.1
31	Paikallisten marjojen osuus kokonaiskulutuksesta (keskiarvo altistumsperiodille)	F.berries.loc		Aikuinen	0.5
32	Paikallisten marjojen osuus kokonaiskulutuksesta (keskiarvo altistumsperiodille)	F.berries.loc		Lapsi	0.5
33	Peseytymiseen kuluva aika	Ed.bw	h/vrk	Aikuinen	0.25
34	Peseytymiseen kuluva aika	Ed.bw	h/vrk	Lapsi	0.25
35	Altistumistiheys, pesuvesi	Ef.bw	vrk/vuosi	Aikuinen	365
36	Altistumistiheys, pesuvesi	Ef.bw	vrk/vuosi	Lapsi	365
37	Pesuvedelle altistuvan ihon osuus	F.skin.bw		Aikuinen	1
38	Pesuvedelle altistuvan ihon osuus	F.skin.bw		Lapsi	1
39	Uidessa niellyn pintaveden määrä	Ir.sw	l/uintikerta	Aikuinen	0.05
40	Uidessa niellyn pintaveden määrä	Ir.sw	l/uintikerta	Lapsi	0.05
41	Uintikerran kesto	Ed.sw	h/uintikerta	Aikuinen	0.25
42	Uintikerran kesto	Ed.sw	h/uintikerta	Lapsi	0.25
43	Altistumistiheys, pintavesi	Ef.sw	Uintikerta/vuosi	Aikuinen	30
44	Altistumistiheys, pintavesi	Ef.sw	Uintikerta/vuosi	Lapsi	30
45	Pintavedelle altistuvan ihon osuus	F.skin.sw		Aikuinen	1
46	Pintavedelle altistuvan ihon osuus	F.skin.sw		Lapsi	1

Mallifunktio

+ Näytä koodi- Piilota koodi library(OpasnetBaseUtils) library(xtable) page <- "Op_fi2814" C.air.outdoor <- 0.118/1000 C.air.indoor <- 0.118/1000 C.soil <- 2645 C.dw <- 0.06 C.lfc <- 2.5 C.tfc <- 0.9 C.berries <- 1.9 C.bw <- 0.06 C.sw <- 0.003 Abs.skin.soil <- 0.003 Kp <- 0.001 #Riskinarviointi <- function(Whatever){ data <- op_baseGetData("opasnet_base", page) data <- tidy(data) #############Riskinarvioinnin laskenta #Aikuinen, haitta-aineen päivittäinen saanti data.a <- data[data$Henkilö == "Aikuinen",] Bw.a <- data.a[data.a$Lyhenne == "Bw",][,5] At.a <- data.a[data.a$Lyhenne == "At",][,5] ##Hengitystiet ###Ilma Ir.air.a <- data.a[data.a$Lyhenne == "Ir.air",][,5] Ef.air.indoor.a <- data.a[data.a$Lyhenne == "Ef.air.indoor",][,5] Air.inh.a <- ((C.air.outdoor * (24-Ef.air.indoor.a/24)) + (C.air.indoor * (Ef.air.indoor.a/24)) * Ir.air.a) / (Bw.a * At.a) ##Ruuansulatuskanava ###Juomavesi Ir.dw.a <- data.a[data.a$Lyhenne == "Ir.dw",][,5] Ef.dw.a <- data.a[data.a$Lyhenne == "Ef.dw",][,5] Dw.ing.a <- (C.dw * Ir.dw.a * Ef.dw.a) / (Bw.a * At.a) ###Maaperä Ir.soil.a <- data.a[data.a$Lyhenne == "Ir.soil",][,5] Ef.soil.a <- data.a[data.a$Lyhenne == "Ef.soil",][,5] Soil.ing.a <- (C.soil * (Ir.soil.a / 1000000) * Ef.soil.a) / (Bw.a * At.a) ###Vihannekset Ir.lfc.a <- data.a[data.a$Lyhenne == "Ir.lfc",][,5] F.lfc.loc.a <- data.a[data.a$Lyhenne == "F.lfc.loc",][,5] Lfc.ing.a <- (C.lfc * Ir.lfc.a * F.lfc.loc.a) / Bw.a ###Juurekset Ir.tfc.a <- data.a[data.a$Lyhenne == "Ir.tfc",][,5] F.tfc.loc.a <- data.a[data.a$Lyhenne == "F.tfc.loc",][,5] Tfc.ing.a <- (C.tfc * Ir.tfc.a * F.tfc.loc.a) / Bw.a ###Marjat Ir.berries.a <- data.a[data.a$Lyhenne == "Ir.berries",][,5] F.berries.loc.a <- data.a[data.a$Lyhenne == "F.berries.loc",][,5] Berries.ing.a <- (C.berries * Ir.berries.a * F.berries.loc.a) / Bw.a ###Pintavesi Ir.sw.a <- data.a[data.a$Lyhenne == "Ir.sw",][,5] Ef.sw.a <- data.a[data.a$Lyhenne == "Ef.sw",][,5] Sw.ing.a <- (C.sw * Ir.sw.a * Ef.sw.a) / (Bw.a * At.a) ##Iho Skin.sa.a <- (4 * Bw.a + 7) / (Bw.a + 90) * 10000 #Lasketaan ihon kokonaispinta-ala ###Maaperä Soil.ad.skin.a <- data.a[data.a$Lyhenne == "Soil.ad.skin",][,5] F.skin.soil.a <- data.a[data.a$Lyhenne == "F.skin.soil",][,5] Soil.derm.a <- (C.soil * Soil.ad.skin.a * 0.000001 * Skin.sa.a * F.skin.soil.a *Abs.skin.soil * Ef.soil.a) / (Bw.a * At.a) ###Pesuvesi Ef.bw.a <- data.a[data.a$Lyhenne == "Ef.bw",][,5] Ed.bw.a <- data.a[data.a$Lyhenne == "Ef.bw",][,5] F.skin.bw.a <- data.a[data.a$Lyhenne == "F.skin.bw",][,5] Bw.derm.a <- (C.bw * 0.001 * Kp * Ef.bw.a * Ed.bw.a * F.skin.bw.a) / (Bw.a * At.a) ###Pintavesi Ef.sw.a <- data.a[data.a$Lyhenne == "Ef.sw",][,5] Ed.sw.a <- data.a[data.a$Lyhenne == "Ef.sw",][,5] F.skin.sw.a <- data.a[data.a$Lyhenne == "F.skin.sw",][,5] Sw.derm.a <- (C.sw * 0.001 * Kp * Ef.sw.a * Ed.sw.a * F.skin.sw.a) / (Bw.a * At.a) Tot.intake.a <- Air.inh.a + Dw.ing.a + Soil.ing.a + Lfc.ing.a + Tfc.ing.a + Berries.ing.a + Sw.ing.a + Soil.derm.a + Bw.derm.a + Sw.derm.a #Lapsi, haitta-aineen päivittäinen saanti data.c <- data[data$Henkilö == "Lapsi",] Bw.c <- data.c[data.c$Lyhenne == "Bw",][,5] At.c <- data.c[data.c$Lyhenne == "At",][,5] ##Hengitystiet ###Ilma Ir.air.c <- data.c[data.c$Lyhenne == "Ir.air",][,5] Ef.air.indoor.c <- data.c[data.c$Lyhenne == "Ef.air.indoor",][,5] Air.inh.c <- ((C.air.outdoor * (24-Ef.air.indoor.c/24)) + (C.air.indoor * (Ef.air.indoor.c/24)) * Ir.air.c) / (Bw.c * At.c) ##Ruuansulatuskanava ###Juomavesi Ir.dw.c <- data.c[data.c$Lyhenne == "Ir.dw",][,5] Ef.dw.c <- data.c[data.c$Lyhenne == "Ef.dw",][,5] Dw.ing.c <- (C.dw * Ir.dw.c * Ef.dw.c) / (Bw.c * At.c) ###Maaperä Ir.soil.c <- data.c[data.a$Lyhenne == "Ir.soil",][,5] Ef.soil.c <- data.c[data.c$Lyhenne == "Ef.soil",][,5] Soil.ing.c <- (C.soil * (Ir.soil.c / 1000000) * Ef.soil.c) / (Bw.c * At.c) ###Vihannekset Ir.lfc.c <- data.c[data.c$Lyhenne == "Ir.lfc",][,5] F.lfc.loc.c <- data.c[data.c$Lyhenne == "F.lfc.loc",][,5] Lfc.ing.c <- (C.lfc * Ir.lfc.c * F.lfc.loc.c) / Bw.c ###Juurekset Ir.tfc.c <- data.c[data.c$Lyhenne == "Ir.tfc",][,5] F.tfc.loc.c <- data.c[data.c$Lyhenne == "F.tfc.loc",][,5] Tfc.ing.c <- (C.tfc * Ir.tfc.c * F.tfc.loc.c) / Bw.c ###Marjat Ir.berries.c <- data.c[data.c$Lyhenne == "Ir.berries",][,5] F.berries.loc.c <- data.c[data.c$Lyhenne == "F.berries.loc",][,5] Berries.ing.c <- (C.berries * Ir.berries.c * F.berries.loc.c) / Bw.c ###Pintavesi Ir.sw.c <- data.c[data.c$Lyhenne == "Ir.sw",][,5] Ef.sw.c <- data.c[data.c$Lyhenne == "Ef.sw",][,5] Sw.ing.c <- (C.sw * Ir.sw.c * Ef.sw.c) / (Bw.c * At.c) ##Iho Skin.sa.c <- (4 * Bw.c + 7) / (Bw.c + 90) * 10000 #Lasketaan ihon kokonaispinta-ala ###Maaperä Soil.ad.skin.c <- data.c[data.c$Lyhenne == "Soil.ad.skin",][,5] F.skin.soil.c <- data.c[data.c$Lyhenne == "F.skin.soil",][,5] Soil.derm.c <- (C.soil * Soil.ad.skin.c * 0.000001 * Skin.sa.c * F.skin.soil.c *Abs.skin.soil * Ef.soil.c) / (Bw.c * At.c) ###Pesuvesi Ef.bw.c <- data.c[data.c$Lyhenne == "Ef.bw",][,5] Ed.bw.c <- data.c[data.c$Lyhenne == "Ef.bw",][,5] F.skin.bw.c <- data.c[data.c$Lyhenne == "F.skin.bw",][,5] Bw.derm.c <- (C.bw * 0.001 * Kp * Ef.bw.c * Ed.bw.c * F.skin.bw.c) / (Bw.c * At.c) ###Pintavesi Ef.sw.c <- data.c[data.c$Lyhenne == "Ef.sw",][,5] Ed.sw.c <- data.c[data.c$Lyhenne == "Ef.sw",][,5] F.skin.sw.c <- data.c[data.c$Lyhenne == "F.skin.sw",][,5] Sw.derm.c <- (C.sw * 0.001 * Kp * Ef.sw.c * Ed.sw.c * F.skin.sw.c) / (Bw.c * At.c) Tot.intake.c <- Air.inh.c + Dw.ing.c + Soil.ing.c + Lfc.ing.c + Tfc.ing.c + Berries.ing.c + Sw.ing.c + Soil.derm.c + Bw.derm.c + Sw.derm.c ###########Ilmoitetaan tulokset cat("Laskennassa käytetty haitta-aineen pitoisuus väliaineessa (mg/kg) \n") Pitoisuus <- data.frame(Ulkoilma=C.air.outdoor, Sisäilma=C.air.indoor, Juomavesi=C.dw, Maaperä=C.soil, Vihannekset=C.lfc, Juurekset=C.tfc, Marjat=C.berries, Pesuvesi=C.bw, Pintavesi=C.sw) print(xtable(Pitoisuus), type = 'html') cat("Keskimääräinen päivittäinen altistuminen haitta-aineelle (µg/kg), aikuinen \n") Kokonaissaanti.a <- data.frame(Ilma=Air.inh.a, Juomavesi=Dw.ing.a, Maaperä.rsk=Soil.ing.a, Vihannekset=Lfc.ing.a, Juurekset=Tfc.ing.a, Marjat=Berries.ing.a, Pintavesi.rsk=Sw.ing.a, Maaperä.iho=Soil.derm.a, Pesuvesi.iho=Bw.derm.a, Pintavesi.iho=Sw.derm.a, Kokonaissaanti=Tot.intake.a)*1000 print(xtable(Kokonaissaanti.a), type = 'html') cat("Keskimääräinen päivittäinen altistuminen haitta-aineelle (µg/kg), lapsi \n") Kokonaissaanti.c <- data.frame(Ilma=Air.inh.c, Juomavesi=Dw.ing.c, Maaperä.rsk=Soil.ing.c, Vihannekset=Lfc.ing.c, Juurekset=Tfc.ing.c, Marjat=Berries.ing.c, Pintavesi.rsk=Sw.ing.c, Maaperä.iho=Soil.derm.c, Pesuvesi.iho=Bw.derm.c, Pintavesi.iho=Sw.derm.c, Kokonaissaanti=Tot.intake.c)*1000 #Kokonaissaanti.c <- data.frame(Altistusreitti = c("Ilma", "Juomavesi", ...), Henkilö = , Arvo = c(Air.inh.c, Dw.ing.c, Soil.ing.c...) #, Vihannekset=Lfc.ing.c, Juurekset=Tfc.ing.c, Marjat=Berries.ing.c, Pintavesi.rsk=Sw.ing.c, Maaperä.iho=Soil.derm.c, Pesuvesi.iho=Bw.derm.c, Pintavesi.iho=Sw.derm.c, Kokonaissaanti=Tot.intake.c)*1000 print(xtable(Kokonaissaanti.c), type = 'html') cat("Osuus kokonaissaannista, aikuinen \n") Saanti.osuus.a <- Kokonaissaanti.a / Kokonaissaanti.a$Kokonaissaanti print(xtable(Saanti.osuus.a), type = 'html') cat("Osuus kokonaissaannista, lapsi \n") Saanti.osuus.c <- Kokonaissaanti.c / Kokonaissaanti.c$Kokonaissaanti print(xtable(Saanti.osuus.c), type = 'html') #return(Kokonaissaanti.c) #}

Vaihtoehtokoodi

Tämä koodi käyttää funktiota haluttujen rivien slaissaamiseen.

+ Näytä koodi- Piilota koodi cat("Ladataan funktioita. \n") library(OpasnetBaseUtils) library(xtable) ########################## tf is a function that takes one large table, finds matching rows and #### produces a standard-format vector that can be operated with easily. Parameters: # condition = a vector with conditions for each condition column. The length must be ncol(col.condition). ############ External parameters (these must be defined as objects before the function is run; ##### they are not given in the function but they are used by it): # data = a data frame that contains the data # col.condition = vector with names of columns that contain the indices that must match in the data # col.common = vector with names of columns that are used as indices in calculations. # col.result = the column that contains the actual values. tf <- function(condition, ...) { # Take the part of data that fits the condition. Drop all columns but col.common and col.result. out <- data[data[col.condition] == condition, c(col.common, col.result)] # Merge the selected data with all col.common and drop possible extra rows so that the outcome always has the same size. out <- merge(unique(data[col.common]), out, all = TRUE) out <- out[rownames(unique(out[col.common])), ] return(out[, col.result]) } #################### exposure is a function that calculates exposure based on environmental concentrations. exposure <- function(C.air.outdoor, C.air.indoor, data, ...) { Air.inh <- ((C.air.outdoor * (24 - tf("Ef.air.indoor") / 24)) + (C.air.indoor * tf("Ef.air.indoor") / 24) * tf("Ir.air")) / (tf("Bw") * tf("At")) ### TÄHÄN VÄLIIN LOPPUOSA LASKENNASTA ERI ALTISTUSREITEILLÄ # Laskennan jälkeen lisätään yhteiset sarakkeet col.common ja pinotaan kaikki altistusreitit allekkain data.frameen. common <- merge(unique(data[col.common]), data.frame(temp=1))[-(length(col.common)+1)] #colnames(common) <- col.common # Seuraavassa toistokerroin 2 on eri altistusreittien määrä. out <- data.frame(common, Altistusreitti = rep(c("Inhalaatio", "Etc"), each = 2), Result = c(Air.inh, 0, 0)) return(out) } cat("Ladataan dataa. \n") page <- "Op_fi2814" C.air.outdoor <- 0.118/1000 C.air.indoor <- 0.118/1000 C.soil <- 2645 C.dw <- 0.06 C.lfc <- 2.5 C.tfc <- 0.9 C.berries <- 1.9 C.bw <- 0.06 C.sw <- 0.003 Abs.skin.soil <- 0.003 Kp <- 0.001 data <- tidy(op_baseGetData("opasnet_base", page)) ### Määritellään ulkoiset parametrit kyseisen taulukon muokkaamista varten. col.common <- "Henkilö" col.result <- "Result" col.condition = "Yksikkö" # Tämä on vain tilapäinen demonstraatio siitä että eri sarakkeita voi käyttää. tf("vrk") col.condition <- "Lyhenne" # Varsinaisesti meitä kiinnostavat sarakkeen Lyhenne määrittämät asiat. tf("Bw") tf("Ir.air") exposure(C.air.outdoor = 10, C.air.indoor = 50, data)

Arvioinnin tulos

Haitta-aineen pitoisuus ulkoilmassa (µg/m3): Haitta-aineen pitoisuus sisäilmassa (µg/m3): Haitta-aineen pitoisuus juomavedessä (mg/l): Haitta-aineen pitoisuus maaperässä (mg/kg): Haitta-aineen pitoisuus vihanneksissa (mg/kg tuorepainoa): Haitta-aineen pitoisuus juureksissa (mg/kg tuorepainoa): Haitta-aineen pitoisuus marjoissa (mg/kg tuorepainoa): Haitta-aineen pitoisuus pesuvedessä (mg/l): Haitta-aineen pitoisuus pintavedessä (mg/l): Turvallisen päivittäisen saannin viitearvo (µg/kg/d): + Näytä koodi- Piilota koodi library(OpasnetBaseUtils) library(xtable) #juomavesi(page, C.dw, Ir.dw, Ef.dw, Bw, At) #Intake.dw <- juomavesi(page, C.dw, Ir.dw, Ef.dw, Bw, At) #print(xtable(Intake.dw), type = 'html') C.dw