test(){ covariate(dose_grp) # taking amount of dose from input tau = 12 # dose interval ndoses = 6 # number of doses sequence{N = 10000} # initialization of bacterial growth # Simulation of multiple doses Ndose = floor(t/tau+1) < ndoses ? floor(t/tau+1) : ndoses # actual number of doses md_time = t - tau * (Ndose - 1) # actual time after multiple doses D = dose_grp * (1 - exp(Ndose * -Ke * tau)) / (1 - exp( -Ke * tau)) * exp(-Ke * md_time) #PK prediction secondary(accum_factor = 1 /( 1- exp(-Ke*tau))) # accumulation factor deriv(N = Kg * N - Kk * D * N) # Turnover model error(NEps = 0.1) # initial estimate of residual error observe(NObs = N + N*NEps) # multiplicative residual error model fixef(Ke = c(, 0.552466, )) # initial estimate of elimination rate constant fixef(Kg = c(, 0.101837, )) # ditto for growth rate constant fixef(Kk = c(, 0.824477, )) # ditto for kill rate constant }