[R] lmer fixed effects, SE, t . . . and p

[John Sorkin]

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}}