Väestön kohdekohtainen ympäristöperäisen haitta-ainealtistumisen arviointi
Opasnet Suomista
Versio hetkellä 5. tammikuuta 2012 kello 09.10 – tehnyt Virpi Kollanus (keskustelu | muokkaukset) (→R-koodi)
R-koodi
library(OpasnetBaseUtils) library(xtable) Bw <- c(Aikuinen=Bw.adult, Lapsi=Bw.child) # Haitta-aineen päivittäinen saanti juomaveden kautta (mg/kd/d): Ir.dw <- c(Aikuinen=Ir.dw.adult, Lapsi=Ir.dw.child) Intake.dw <- (C.dw*Ir.dw*Ef.dw)/(Bw*365) # Haitta-aineen päivittäinen saanti maaperän kautta (mg/kg/d) Ir.soil <- c(Aikuinen=Ir.soil.adult, Lapsi=Ir.soil.child) Intake.soil <- (C.soil*Ir.soil/1000000*Ef.soil)/(Bw*365) # Haitta-aineen päivittäinen saanti ravintokasvin kautta (mg/kg/d) Ir.food.plant <- c(Aikuinen=Ir.food.plant.adult, Lapsi = Ir.food.plant.child) Intake.food.plant <- (C.food.plant*Ir.food.plant*Ef.food.plant)/(Bw*365) cat("Keskimääräinen päivittäinen altistuminen haitta-aineelle (µg/kg)\n") Intake <-data.frame(Juomavesi=Intake.dw, Maaperä=Intake.soil, Ravinto=Intake.food.plant, Kokonaissaanti=(Intake.dw+Intake.soil+Intake.food.plant)) print(xtable(Intake*1000), type = 'html') cat("Eri altistumisreittien osuus kokonaissaannista (%)\n") Intake.fraction <- Intake/Intake$Kokonaissaanti print(xtable(Intake.fraction*100), type = 'html') cat("Vaaraosamäärä\n") HQ.adult <- Intake[1,4]*1000/TDI HQ.child <- Intake[2,4]*1000/TDI HQ <- data.frame(Aikuinen=HQ.adult, Lapsi=HQ.child) print(xtable(HQ), type = 'html') |