On Mon, Jul 13, 2015 at 3:04 AM, Patty Haaem via R-help <r-help@r-project.org> wrote: > 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
First, you should not post in HTML. Second, you should put the argument in quotes: na.action="na.include" Sarah > > > [[alternative HTML version deleted]] > -- Sarah Goslee http://www.functionaldiversity.org ______________________________________________ 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.