[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