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]]
______________________________________________
[email protected] 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.