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

Mitchell Maltenfort mmalten at gmail.com
Fri Sep 10 05:04:29 CEST 2010


Try coef(summary(fit3))


On Thu, Sep 9, 2010 at 11:00 PM, John Sorkin
<jsorkin at grecc.umaryland.edu> wrote:
> 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 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.
>



More information about the R-help mailing list