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