Hi folks, I am trying to generate 95% confidence intervals for a gls model using predict.nlme with R version 2.5.1 (2007-06-27) . nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1-83. I have looked in help, and I can do it for lm and glm models, and I can generate simple predictions for lme models with various levels -- I am familiar with the basics. Is there a way to get prediction intervals for gls models? My "best" model uses varPower(), so I am reluctant to fall back on lm predictions. Thank you, Hank Dr. Hank Stevens, Associate Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/~stevenmh/ http://www.muohio.edu/ecology/ http://www.muohio.edu/botany/ "E Pluribus Unum"
Martin,
Have you checked ?intervals.gls
This intervals are approximate, but this would be the obvious starting
point to me.
Joris
"Martin Henry H.
Stevens"
<HStevens at muohio. To
edu> R-Help <r-help at
stat.math.ethz.ch>
Sent by: cc
r-help-bounces at st
at.math.ethz.ch Subject
[R] lme or gls prediction intervals
24/07/2007 19:22
Hi folks,
I am trying to generate 95% confidence intervals for a gls model
using predict.nlme
with
R version 2.5.1 (2007-06-27)
. nlme: Linear and Nonlinear Mixed Effects Models. R package version
3.1-83.
I have looked in help, and I can do it for lm and glm models, and I
can generate simple predictions for lme models with various levels --
I am familiar with the basics.
Is there a way to get prediction intervals for gls models? My "best"
model uses varPower(), so I am reluctant to fall back on lm predictions.
Thank you,
Hank
Dr. Hank Stevens, Associate Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056
Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.cas.muohio.edu/~stevenmh/
http://www.muohio.edu/ecology/
http://www.muohio.edu/botany/
"E Pluribus Unum"
______________________________________________
R-help at 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.
Hi Joris, Thank you for the reply. I also realize that I meant confidence intervals, not prediction intervals! I am trying to do something analogous to predict(model.gls, newdata, interval="confidence") but predict.gls does not have an interval argument. I am guessing that this is because it is subtle and complicated to figure out what the appropriate hypothesis is, but I was hoping to get some feedback. Cheers, Hank On Jul 25, 2007, at 3:21 AM, joris.dewolf at cropdesign.com wrote:> > > Martin, > > > Have you checked ?intervals.gls > > > This intervals are approximate, but this would be the obvious starting > point to me. > > > Joris > > > > > > > > > > > > "Martin Henry H. > Stevens" > > <HStevens at muohio. To > edu> R-Help <r- > help at stat.math.ethz.ch> > Sent > by: cc > r-help-bounces at st > at.math.ethz.ch > Subject > [R] lme or gls prediction > intervals > > 24/07/2007 19:22 > > > > > > > > > Hi folks, > I am trying to generate 95% confidence intervals for a gls model > using predict.nlme > with > R version 2.5.1 (2007-06-27) > . nlme: Linear and Nonlinear Mixed Effects Models. R package version > 3.1-83. > > I have looked in help, and I can do it for lm and glm models, and I > can generate simple predictions for lme models with various levels -- > I am familiar with the basics. > > Is there a way to get prediction intervals for gls models? My "best" > model uses varPower(), so I am reluctant to fall back on lm > predictions. > > Thank you, > > Hank > > > Dr. Hank Stevens, Associate Professor > 338 Pearson Hall > Botany Department > Miami University > Oxford, OH 45056 > > Office: (513) 529-4206 > Lab: (513) 529-4262 > FAX: (513) 529-4243 > http://www.cas.muohio.edu/~stevenmh/ > http://www.muohio.edu/ecology/ > http://www.muohio.edu/botany/ > > "E Pluribus Unum" > > ______________________________________________ > R-help at 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. > >Dr. Hank Stevens, Associate Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/~stevenmh/ http://www.muohio.edu/ecology/ http://www.muohio.edu/botany/ "E Pluribus Unum" If you send an attachment, please try to send it in a format anyone can read, such as PDF, text, Open Document Format, HTML, or RTF. Please try not to send me MS Word or PowerPoint attachments- Why? See: http://www.gnu.org/philosophy/no-word-attachments.html