[R] lme or gls prediction intervals

[Martin Henry H. Stevens]

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