# (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)
# (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) check nesting
nested_model(dataPBK_C4)
# (2bis)
dataPBK_C42 <- dataPBK(
object = dataCompartment4,
col_time = "temps",
col_replicate = "replicat",
col_exposure = "condition",
col_compartment = c("intestin", "reste", "caecum", "cephalon"),
time_accumulation = 7,
ku_nest = c(1,1,1,1), # No Change here
ke_nest = c(1,1,1,1), # No Change here
k_nest = matrix(c(
c(0,1,1,1),
c(0,0,1,1),
c(0,0,0,0),
c(0,0,0,0)),
ncol=4,nrow=4,byrow=TRUE) # Remove
)
# (3bis) re-checking nesting
nested_model(dataPBK_C42)
Run the code above in your browser using DataLab