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