[R] lmer p-values for fixed effects missing

Steven McKinney smckinney at bccrc.ca
Wed Nov 19 02:29:48 CET 2008


> -----Original Message-----
> From: r-help-bounces at r-project.org on behalf of UKaraoz
> Sent: Tue 11/18/2008 2:16 PM
> To: r-help at r-project.org
> Subject: [R]  lmer p-values for fixed effects missing
>  
> 
> I am trying to replicate the repeated measures example from Dr.Faraway's book
> (Extending the linear model with R) as follows:
> 
> data(vision)
> vision$npower <- rep(1:4,14)
> 
> mmod <-lmer(acuity~power+(1|subject)+(1|subject:eye),vision)
> 
> 
> When I look at the fixed effects p-value, it is missing. Am I missing
> something here?
> Thanks.
> 
> Fixed effects:
>             Estimate Std. Error t value
> (Intercept) 112.6429     2.2348   50.40
> power6/18     0.7857     1.5400    0.51
> power6/36    -1.0000     1.5400   -0.65
> power6/60     3.2857     1.5400    2.13
> 



See the discussion titled "lmer, p-values and all that" at
https://stat.ethz.ch/pipermail/r-help/2006-May/094765.html
and a collection of related discussions at
http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests

In a nutshell:
The complexity of fully crossed mixed effects models
makes degrees of freedom and other model attributes
more difficult to calculate compared to fixed effects 
models.  

Research and debate is ongoing as to the best way to
handle model assessment, not all classical summary
statistics and associated p-values map to lmer models 
in a straightforward fashion.


Steven McKinney

Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre

email: smckinney +at+ bccrc +dot+ ca

tel: 604-675-8000 x7561

BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C. 
V5Z 1L3
Canada



> 
> -- 
> View this message in context: http://www.nabble.com/lmer-p-values-for-fixed-effects-missing-tp20569114p20569114.html
> Sent from the R help mailing list archive at Nabble.com.
> 
> ______________________________________________
> 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