Displaying 1 result from an estimated 1 matches for "rinln".
Did you mean:
rinen
2011 Sep 05
1
p values greater than 1 from lme4
...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...