[R] p values greater than 1 from lme4
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
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:
> 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)
> attach (x)
> 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
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.
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