[R-sig-ME] nlme library
Martin Henry H. Stevens
HStevens at muohio.edu
Wed Jun 25 11:39:14 CEST 2008
The VarCorr extractor function should help:
influent = pdLogChol(1)
(Intercept) 63.323 7.9576
Residual 42.658 6.5313
On Jun 25, 2008, at 5:22 AM, roberta varriale wrote:
> 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:
> m1 <- lme(y ~ 1, data = Mississippi, random = ~ 1|influent)
> 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-
> 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
> 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
Dr. Hank Stevens, Associate Professor
338 Pearson Hall
Oxford, OH 45056
Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
"E Pluribus Unum"
"I love deadlines. I love the whooshing noise they make as they go by."
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
More information about the R-sig-mixed-models