R help - May 2006 - Incomplete Output from lmer{lme4}

I'm still relatively new to R, so my apologies if this is covered
somewhere.  I've been running some mixed-effect models in R using
lme{nlme}, but read in Faraway's recent book, Extending the Linear Model
with R, that lmer in package lme4 is a much improved version.  I tried
using this approach, but the output for the fixed effects doesn't report
a p-value or the degrees of freedom (unlike the examples listed in
Faraway's book, which I tried and got the same incomplete output).  Any
idea how I can get the complete output?

In a related problem, I've been having trouble with the lines and abline
in the graphics package.  Occasionally, they fail to add the requested
line to the figure, but R doesn't return any error (it's as if executing
the command does nothing).  I've replicated this problem with code and
data directly from R texts such as Crawley's Statistics: An Introduction
Using R, so I don't think it is operator error (but it could be).  Any
ideas?

Thanks,

Jeff--
________________________________________________________
 
Jeffrey J. Lusk, Ph.D.
Postdoctoral Research Associate
Department of Forestry & Natural Resources
Purdue Climate Change Research Center
Purdue University
715 West State Street; Pfendler Hall
West Lafayette, IN 47907
Phone:  765-494-9701
E-mail:  jjlusk at purdue.edu
Web: http://www.geocities.com/jefflusk2002

You may find it helpful to read the posting guide (http://www.R- 
project.org/posting-guide.html)

On May 18, 2006, at 11:40 AM, Lusk, Jeffrey J wrote:
> I'm still relatively new to R, so my apologies if this is covered
> somewhere.  I've been running some mixed-effect models in R using
> lme{nlme}, but read in Faraway's recent book, Extending the Linear  
> Model
> with R, that lmer in package lme4 is a much improved version.  I tried
> using this approach, but the output for the fixed effects doesn't  
> report
> a p-value or the degrees of freedom (unlike the examples listed in
> Faraway's book, which I tried and got the same incomplete output).   
> Any
> idea how I can get the complete output?
P-values for model coefficients are based on a Monte Carlo resampling  
using mcmcsamp. Do ?mcmcsamp. Results are best summarized with the  
summary function in the coda package.
P-values for anova are done with stepwise model building
anova(larger.model, smaller.model)
>
> In a related problem, I've been having trouble with the lines and  
> abline
> in the graphics package.  Occasionally, they fail to add the requested
> line to the figure, but R doesn't return any error (it's as if  
> executing
> the command does nothing).  I've replicated this problem with code and
> data directly from R texts such as Crawley's Statistics: An  
> Introduction
> Using R, so I don't think it is operator error (but it could be).  Any
> ideas?
Couldn't tell you without and example.
Cheers,
Hank
>
> Thanks,
>
> Jeff--
> ________________________________________________________
>
> Jeffrey J. Lusk, Ph.D.
> Postdoctoral Research Associate
> Department of Forestry & Natural Resources
> Purdue Climate Change Research Center
> Purdue University
> 715 West State Street; Pfendler Hall
> West Lafayette, IN 47907
> Phone:  765-494-9701
> E-mail:  jjlusk at purdue.edu
> Web: http://www.geocities.com/jefflusk2002
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting- 
> guide.html
Dr. M. Hank H. Stevens, Assistant Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056

Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.cas.muohio.edu/~stevenmh/
http://www.muohio.edu/ecology/
http://www.muohio.edu/botany/
"E Pluribus Unum"

Lusk, Jeffrey J <jjlusk <at> purdue.edu> writes:
> 
> using this approach, but the output for the fixed effects doesn't
report
> a p-value or the degrees of freedom (unlike the examples listed in
> Faraway's book, which I tried and got the same incomplete output).  Any
> idea how I can get the complete output?
This is a temporary state of lmer(lme4): currently p-values are not printed, 
because Douglas Bates is re-evaluating degrees of freedom estimation in his 
model, and as long he is not sure, he takes the "better missing than
wrong"
stand. 

Hope I quoted him right, the original citation is so hidden that I remember it, 
but cannot find it any more.

Dieter