# \donttest{
# 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"))
# plot solution based on output
plt.out <- solve_full_pregnancy(chem.name = "genistein",
dose = 1, plt = TRUE)
# }
Run the code above in your browser using DataLab