Dear All,I am trying to fit one compartment IV bolus model on pharmacokinetic
data using phenoModel function in nlme package, based on a tutorial entitled
"Development of population PK model using R- Case study I". The codes are as
fallowing:
library(nlme)mydata.grp <- groupedData(CONC~TIME|CID,data=mydata)mydata.fit <
-nlme(CONC~phenoModel(CID,TIME,AMT,lCl,lV),fixed=lCl+lV~1,random=pdDiag(lCl+lV~1),
data=mydata.grp,start=c(-5,0),weight=varConstPower(const=1,fixed=list(power=1)),na.action=na.include,naPattern=~!is.na(CONC))
when I run above codes, I get the following error:Error in nlme.formula(CONC ~
phenoModel(CID, TIME, AMT, lCl, lV), fixed = lCl + : object 'na.include' not
found
could you please tell me, how should I correct the codes?Thanks in advanceElham
Haem
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
