[R-sig-ME] nlme library

[Martin Henry H. Stevens]

Hi Roberta,
The VarCorr extractor function should help:
 > VarCorr(m1)
influent = pdLogChol(1)
Variance StdDev
(Intercept) 63.323   7.9576
Residual    42.658   6.5313
 > as.numeric(VarCorr(m1)[1,2])
[1] 7.9576

Hank

On Jun 25, 2008, at 5:22 AM, roberta varriale wrote:
> Hello,
> I'm Roberta.
> I'm working with R since some months, but I still have some problems.
> At the moment I'm using the library nlme, in particular the function
> lme and I have a question.
>
> Just as an example I write the code:
> library(nlme)
> library(SASmixed)
> data(Mississippi)
> m1 <- lme(y ~ 1, data = Mississippi, random = ~ 1|influent)
> summary(m1)
>
>
> and I obtain these results:
>
> Linear mixed-effects model fit by REML
>  Data: Mississippi
>        AIC      BIC    logLik
>   258.3511 263.1017 -126.1756
>
> Random effects:
>  Formula: ~1 | influent
>         (Intercept) Residual
> StdDev:    7.957598 6.531317
>
> Fixed effects: y ~ 1
>                Value Std.Error DF  t-value p-value
> (Intercept) 21.22312  3.429034 31 6.189241       0
>
> Standardized Within-Group Residuals:
>         Min          Q1         Med          Q3         Max
> -1.91443594 -0.53645956 -0.03217348  0.83713515  1.95823685
>
> Number of Observations: 37
> Number of Groups: 6
>
>
> Now, I would like to recall the parameter "Random effects StdDev
> (Intercept)", equal to 7.957598, but I'm not able to find its "R- 
> name".
> I read the Documentation on R web site and I found the "name" of other
> parameters (such as ...$coeff$fixed and ...$sigma) but I really  
> need that
> one.
> Someone can help me, please?
>
> Thanks in advance,
> best regards,
> Roberta Varriale
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

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.users.muohio.edu/harkesae/
http://www.cas.muohio.edu/ecology
http://www.muohio.edu/botany/
"E Pluribus Unum"

"I love deadlines. I love the whooshing noise they make as they go by."
                                             (Douglas Adams)


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.
Why? See:  http://www.gnu.org/philosophy/no-word-attachments.html