Frazer Sinclair
2011-Jan-22 21:39 UTC
[R] Using plotLMER.fnc for lmer models with crossed random effects
Hi I am attempting to use plotLMER.fnc from the languageR library to plot partial effects of models fited with lmer. This appears to be straightforward for some of my models, e.g:> mod<-lmer(NoBranches~PCA1+SP+(1|ProvCode)+(1|SoilZone),REML=FALSE,data=treedata1)> plotLMER.fnc(mod)However, when a model contains crossed random effects, I get an error message. e.g:> mod<lmer(NoBranches~PCA1+SP+(1|ProvCode)+(1|SoilZone/ParcelleID),REML=FALSE,data=treedata1)> plotLMER.fnc(mod)Error in 1:(which(predictors =ranefnames[1]) - 1) : argument of length 0 Any suggestions greatly appreciated. Frazer Sinclair [[alternative HTML version deleted]]