# (1) load data file
data("dataCompartment4")
# (2) prepare data set
dataPBK_C4 <- dataPBK(
object = dataCompartment4,
col_time = "temps",
col_replicate = "replicat",
col_exposure = "condition",
col_compartment = c("intestin", "reste", "caecum", "cephalon"),
time_accumulation = 7)
# (3) run Bayesian fitting: <5 sec to be executed
# 1 chain and 10 iterations is fast to run but provide
# bad goodness-of-fit
fitPBK_C4_FASTbadGOF <- fitPBK(dataPBK_C4, chains = 1, iter = 10)
# \donttest{
# (3) run Bayesian fitting: > 5 sec to be executed
# 4 chains and 2000 iterations provides better estimates
fitPBK_C4 <- fitPBK(dataPBK_C4, chains = 4, iter = 2000)
# }
Run the code above in your browser using DataLab