[R] an error in nlme package

[Patty Haaem]

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]]