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

R-koodi

Kehon paino, aikuinen (kg): Kehon paino, lapsi (kg): Haitta-aineen pitoisuus juomavedessä (mg/l): Päivittäinen juomaveden kulutus, aikuinen (l): Päivittäinen juomaveden kulutus, lapsi (l): Päivien määrä vuodessa, jolloin juodaan haitta-ainetta sisältävää vettä: Haitta-aineen pitoisuus maaperässä (mg/kg): Päivittäin ruuansulatuskanavaan kulkeutuvan maa-aineksen määrä, aikuinen (mg): Päivittäin ruuansulatuskanavaan kulkeutuvan maa-aineksen määrä, lapsi (mg): Päivien määrä vuodessa, jolloin haitta-ainetta sisältävälle maaperälle altistutaan: Haitta-aineen pitoisuus ravintokasvissa (mg/kg): Päivittäinen ravintokasvin kulutus, aikuinen (kg): Päivittäinen ravintokasvin kulutus, lapsi (kg): Päivien määrä vuodessa, jolloin syödään haitta-ainetta sisältävää ravintokasvia: Turvallisen päivittäisen saannin viitearvo (µg/kg/d): + Näytä koodi- Piilota 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 altistuminenn 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("Vaaraosamäärä\n") HQ.adult <- Intake[1,4]*1000/TDI HQ.adult <-round(HQ.adult,2) HQ.child <- Intake[2,4]*1000/TDI HQ.child <- round(HQ.child,2) HQ <- data.frame(Aikuinen=HQ.adult, Lapsi=HQ.child) print(xtable(HQ), type = 'html')