[R] p values greater than 1 from lme4

Bogaso Christofer bogaso.christofer at gmail.com
Tue Sep 6 20:29:42 CEST 2011

Hello Bolker, Hope you will make available at least the problem and
reasoning to this list. I am also very much interested to see the problem

-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of Ben Bolker
Sent: 06 September 2011 02:58
To: r-help at stat.math.ethz.ch
Subject: Re: [R] p values greater than 1 from lme4

RTSlider <rob.t.slider <at> gmail.com> writes:

> Hello,
> I'm running linear regressions using the following script where I have 
> separated out species using the "IDtotsInLn" identifier
> x<-read.csv('tbl02TOTSInLn_ENV.csv', header=T)  
> x                                          
> attach (x)                                  
> library(lme4)
> rInLn<-lmList(InLn~pMoist | IDtotsInLn, x, pool=F)
> write.table(summary(rInLn)$coefficients, "rInLnPlots.csv")   
> write.table(summary(rInLn)$r.squared, append=T, "rInLnPlots.csv") 
> write.table(summary(rInLn)$df, append=T, "rInLnPlots.csv")
> The script seems to be working for most of the species, but for some 
> it is returning a p value of greater than 1 (e.g. 20). I thought this 
> might be for the few cases where the independent variable remained 
> constant, but found other species where this was not the case and the 
> p value was still much greater than 1.
> Any help would be appreciated
> -RTS

  This is very interesting but practically impossible to solve because it's
not reproducible; is there any chance that you can make the data available?
You can send it directly to me (Ben Bolker -- my e-mail is pretty easy to
find on the web) if you like.

   Ben Bolker

R-help at r-project.org mailing list
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