# \donttest{
library(httk)
# dosing schedule of 1 mg/kg BW/day for 40 weeks
# return solution by hour
out <- solve_full_pregnancy(chem.name = "fipronil",
daily.dose = 1,
doses.per.day = 1,
time.course = seq(0, 40*7, 1/24))
# return solution in chemical amounts for fetal compartments + placenta
maternal_compts <- c('gutlumen', 'gut', 'liver', 'kidney', 'lung', 'ven', 'art',
'adipose','thyroid', 'rest')
fetal_compts <- c(maternal_compts[! maternal_compts %in% c('adipose', 'gutlumen') ],
"brain")
amt.out <- solve_full_pregnancy(chem.name = "fipronil",
daily.dose = 1,
doses.per.day = 1,
time.course = seq(0, 40*7, 1),
track.vars = c(paste0("Af", fetal_compts), "Aplacenta"))
# return solution in concentrations for fetal compartments + placenta
conc.out <- solve_full_pregnancy(chem.name = "fipronil",
daily.dose = 1,
doses.per.day = 1,
time.course = seq(0, 40*7, 1),
track.vars = c(paste0("Cf", fetal_compts), "Cplacenta"))
# }
Run the code above in your browser using DataLab