On Mon, Jul 13, 2015 at 11:19 AM, Bert Gunter <bgunter.4...@gmail.com> wrote: > No, Sarah. > > na.action must be a function, not a character string. But you're > close: there is no na.include function, as the message says. It > should be na.exclude .
My apologies: I was assuming it worked like options(): > class(options()$na.action) [1] "character" > options()$na.action [1] "na.omit" Sarah > > Cheers, > Bert > > > Bert Gunter > > "Data is not information. Information is not knowledge. And knowledge > is certainly not wisdom." > -- Clifford Stoll > > > On Mon, Jul 13, 2015 at 7:35 AM, Sarah Goslee <sarah.gos...@gmail.com> wrote: >> 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 ______________________________________________ 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.
