Consider the following:
> library(nlme)
Loading required package: lattice
> DF <- data.frame(x=rep(c("a", "b"), 4), y=1:8)
> fit <- lme(y~1, random=~1|x, data=DF)
> summary(fit)
Linear mixed-effects model fit by REML
Data: DF
AIC BIC logLik
40.48753 40.32526 -17.24376
Random effects:
Formula: ~1
Dear R experts,
I want to extract standard error from the output of lme, but I found
fix.effects() does not include SE of the coefficients. Many thanks in
advance.
Best regards,
Yu-Kang
_
現在就上 MSN 會員目錄：在線上結交新朋友，找到興趣相投的夥伴
__
