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@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.
______________________________________________
R-help@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.
