[R] Accessing standard errors for coefficients in nlme

[Rui Barradas]

Hello,

What you want is component tTable. Using the first example in ?lme,

library(nlme)

fm1 <- lme(distance ~ age, data = Orthodont)
s <- summary(fm1)
s$tTable[,2]


Hope this helps,

Rui Barradas

Em 31-08-2013 18:31, Dennis Fisher escreveu:
> R 3.0.1
> OS X
>
> Colleagues,
>
> When I run
> 	summary(lme(...))
> in the nmle package, I get output of this sort:
>
> Linear mixed-effects model fit by REML
>   Data: NULL
>        AIC     BIC   logLik
>    50.2764 59.6771 -16.1382
>
> Random effects:
>   Formula: ~1 | SEQIDS
>          (Intercept)  Residual
> StdDev:   0.4806963 0.1730584
>
> Fixed effects: LOG ~ SEQ + PERIOD + TREAT
>                   Value  Std.Error DF   t-value p-value
> (Intercept)  2.1775070 0.13909328 11 15.655012  0.0000
> SEQ1        -0.0294207 0.22605462 10 -0.130149  0.8990
> SEQ2         0.0406411 0.22605462 10  0.179784  0.8609
> SEQ3        -0.0261052 0.22605462 10 -0.115482  0.9103
> PERIOD1     -0.3075528 0.06702522 11 -4.588613  0.0008
> PERIOD2      0.1627682 0.04468348 11  3.642693  0.0039
> TREAT1      -0.0114106 0.03351261 11 -0.340486  0.7399
>
> I want to access the Std.Error values.
> I tried:
> 	summary(lme(…))$coeff
> but this provides only the "Value" column.
>
> I also tried:
> 	str(summary(lme(…)))
> to see if it provided any insights (none, unfortunately).
>
> How can I access those values?  Any help would be greatly appreciated.
>
> Dennis
>
> Dennis Fisher MD
> P < (The "P Less Than" Company)
> Phone: 1-866-PLessThan (1-866-753-7784)
> Fax: 1-866-PLessThan (1-866-753-7784)
> www.PLessThan.com
>
> ______________________________________________
> R-help at r-project.org 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.
>