Dennis Fisher
2013-Aug-31 17:31 UTC
[R] Accessing standard errors for coefficients in nlme
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
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. >