[R] Accessing standard errors for coefficients in nlme
Dennis Fisher
fisher at plessthan.com
Sat Aug 31 19:31:36 CEST 2013
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
More information about the R-help
mailing list