windows Vista
R 2.10.1
(1) How can I get the complete table of for the fixed effects from lmer. As can
be seen from the example below, fixef(fit2) only give the estimates and not the
SE or t value
> fit3<- lmer(y~time + (1|Subject) + (time|Subject),data=data.frame(data))
> summary(fit3)
Linear mixed model fit by REML
Formula: y ~ time + (1 | Subject) + (time | Subject)
Data: data.frame(data)
AIC BIC logLik deviance REMLdev
-126.2 -116.4 70.1 -152.5 -140.2
Random effects:
Groups Name Variance Std.Dev. Corr
Subject (Intercept) 2.9311e+01 5.41396385
Subject (Intercept) 0.0000e+00 0.00000000
time 0.0000e+00 0.00000000 NaN
Residual 8.1591e-07 0.00090328
Number of obs: 30, groups: Subject, 10
Fixed effects:
Estimate Std. Error t value
(Intercept) 14.998216 1.712046 9
time -0.999779 0.000202 -4950
Correlation of Fixed Effects:
(Intr)
time -0.001
> fixef(fit3)
(Intercept) time
14.9982158 -0.9997793
(2) lmer does not give p values or confidence intervals for the fixed effects.
How we are to interpret the estimates given that no p value or CI is given for
the estimates?
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10 North Greene Street
GRECC (BT/18/GR)
Baltimore, MD 21201-1524
(Phone) 410-605-7119
(Fax) 410-605-7913 (Please call phone number above prior to faxing)
Confidentiality Statement:
This email message, including any attachments, is for th...{{dropped:6}}
______________________________________________
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.
