data(MultBolus)
onecompIV <- list(DiffEq=list( #Differential equations
dy1dt = ~ -ke*y1), #Compartment 1
ObsEq=list( #Observation equation
c1 = ~ y1/Vd),
States=c("y1"), #The names of the states in the sequence of DiffEq
Parms=c("ke","Vd"), #Parameter names
LogParms=TRUE, #Estimate the logarithm of the parameters
Init=list(0), #Estimate the Initial states
JAC=TRUE, #Use the Jacobian
SEQ=FALSE) #Use sensitivity equations
MultBolusModel <- nlmeODE(onecompIV,MultBolus)
fit <- nlme(Conc ~ MultBolusModel(ke,Vd,Time,ID),
data = MultBolus, fixed=ke+Vd~1, random = pdDiag(ke~1),
start=c(ke=log(0.1),Vd=log(1)),
control=list(msVerbose=TRUE))
plot(augPred(fit,level=0:1))
Run the code above in your browser using DataLab