Re: [R] augPred plot in nlme library

From: Afshartous, David <afshart_at_exchange.sba.miami.edu>
Date: Fri 08 Sep 2006 - 21:18:13 GMT

 

Deepayan,

Thanks for your suggestion. Here are more details:

I have a grouped data object for repeated measures data just like the Pixel grouped data object on p.42 of Pinheiro and Bates (2000).

comp.adj.UKV.3 <- groupedData(adj.UKV ~ Time | Patient_no/Lisinopril,

	data = comp.adj.UKV.frm, order.groups = F
	#labels = list(x = "Hour", y = "adj.UKV") 
)

i.e., the response is continuous, Time is not treated as a factor, and there exists two factors, one nested within the other (Lisinopril nested

witin patient, similar to Side within Dog on p.42).

I also fit a model very similar to their model:

fm1comp = lme(adj.UKV ~ Time + Time.sq, data = comp.adj.UKV.3, random = list(Patient_no = ~ 1 , Lisinopril = ~ 1) )

However, the command below does not produce the fitted curves from this model,
but rather it seems to be the fitted curves from a linear model.

plot(augPred(fm3comp))

Possibly augPred behaves differently in R than in S, but reading the R help and
trying various other approaches has not solved this.

Thanks!
Dave

-----Original Message-----
From: Deepayan Sarkar [mailto:deepayan.sarkar@gmail.com] Sent: Friday, September 08, 2006 4:37 PM To: Afshartous, David
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] augPred plot in nlme library

On 9/8/06, Afshartous, David <afshart@exchange.sba.miami.edu> wrote:
> All,
>

> I've solved part of the problem below by making sure that the formula
> in the grouped data object is the same as the formula specified within

> lme (this isn't the case in the cited example from Pinheiro & Bates).
>
> However, augPred seems to plot only a linear model instead of the
> polynomial model. Does anyone know how to make sure that augPred
> plots the same model as that specified in the model (as below)?

You are unlikely to get any helpful answers unless you give us more information, as every r-help message asks you to do:

> R-help@stat.math.ethz.ch mailing list
> 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.

-Deepayan



R-help@stat.math.ethz.ch mailing list
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. Received on Sat Sep 09 07:28:06 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 08 Sep 2006 - 22:30:04 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.