#Example based on Pearce et al. (2017):
# Set up means:
params <- parameterize_pbtk(chem.name="zoxamide")
# Nothing changes:
monte_carlo(params)
vary.params <- NULL
for (this.param in names(params)[!(names(params) %in%
c("Funbound.plasma", "pKa_Donor", "pKa_Accept" )) &
!is.na(as.numeric(params))]) vary.params[this.param] <- 0.2
# Most everything varies with CV of 0.2:
monte_carlo(
parameters=params,
cv.params = vary.params)
censored.params <- list(Funbound.plasma = list(cv = 0.2, lod = 0.01))
# Fup is censored below 0.01:
monte_carlo(
parameters=params,
cv.params = vary.params,
censored.params = censored.params)
Run the code above in your browser using DataLab