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