50 Exposure-Response Modeling: Methods and Practical Implementation
floxacin PK data. The function nlmeODE(.) is an interface from the NLMM
function nlme(.) to the ODE solver, taking the ODE and initia l value specifi-
cations from Onecomp.
bdata=bdata[order(bdata$Subject,bdata$Time),]
gdata=groupedData(conc~Time|Subject,bdata)
OneComp <- list(DiffEq=list(
dy1dt = ~ -ka*y1 , #Absoption dynamic
dy2dt = ~ ka*y1-ke*y2), #Central comp. with elimination rate ke
ObsEq=list(
c1 = ~ 0, #No measure for absorption
c2 = ~ y2/Cl*ke), #Observe concentration c2=y2/volume
#volume=ke/Cl, y2=amount
Parms=c("ka","ke","Cl"),
States=c("y1","y2"),
Init=list(0,0)) #Initial states for the ODE
moximodel=nlmeODE(OneComp,gdata) #generate function for nlme() call
jk=nlme(conc~moximodel(k ...