limit_tox <- 10.96
N <- 42
param_pk <- c(2,10,100) # PK parameters ka,CL,V
omega2 <- 0.7 # Inter-individual
omega_a <- 0
preal <- c(0.001,.05,0.1,.2,0.35,0.45)
time_sampling <- seq(0,24,length.out=48)
sigma2 <- c(rep(0.2,length(time_sampling)),rep(0.2,length(time_sampling)))
# sigma2: Additive or proportional error
gen.scen <- scenarios(param_pk,omega2,omega_a,sigma2,preal,limit_tox,time_sampling, N)
gen.scen
Run the code above in your browser using DataLab