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 -- View this message in context: http://r.789695.n4.nabble.com/p-values-greater-than-1-from-lme4-tp3791526p3791526.html Sent from the R help mailing list archive at Nabble.com.
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 > -RTSThis 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