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@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.
