Hi Ben,
I've attached the data and Tinn-R file I've been using to run it,
thank you for the fast response. Please let me know if you see
something odd, I have no idea where to begin looking
-Rob
On Mon, Sep 5, 2011 at 5:30 PM, bbolker [via R]
<ml-node+3792187-1531262488-213...@n4.nabble.com> wrote:
> 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
>
> ______________________________________________
> [hidden email] 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.
>
>
> ________________________________
> If you reply to this email, your message will be added to the discussion
> below:
> http://r.789695.n4.nabble.com/p-values-greater-than-1-from-lme4-tp3791526p3792187.html
> To unsubscribe from p values greater than 1 from lme4, click here.
--
View this message in context:
http://r.789695.n4.nabble.com/p-values-greater-than-1-from-lme4-tp3791526p3792237.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
