使用ezPredict可视化GLM结果

长期潜伏者，第一次海报。通过搜索旧问题，我得到了很多帮助，但这次我的谷歌fu让我失望了。我正在对一种大型哺乳动物进行栖息地选择研究，使用逻辑回归模型将使用过的地点与可用栖息地的随机样本进行比较。使用lme4的我的模型：

GLM <- glmer(u ~ hab * season + hab * sex + hab * location + feeding * hab + feeding:season:hab + location:season:hab + hab:season:location:sex + Z.dist_road * location + (1|bison), family=binomial, data=subset(bison_KDE,bison_KDE$dataquality!="Poor"))

这些因素包括：

u=1表示已使用位置，0表示可用位置

hab=栖息地，6层

季节=夏季或冬季

位置=3个不同的研究区域

喂养=是否在冬季对个体进行补充喂养

性别=男性还是女性

Z.dist_road=标准化的最近道路的距离（以米为单位）。唯一的连续变量

模型运行良好，但由于各种因素之间的相互作用，输出量巨大。下一步我想做的是生成带有标准偏差的条形图，显示根据其他因素选择的栖息地的差异，例如生成类似于以下内容的图形： !

有人建议我尝试使用ez包和ezPredict函数。首先，我建立了一个包含各种可能因素组合的数据框架

tp<-expand.grid(unique(bison_KDE$sex), unique(bison_KDE$location), unique(bison_KDE$season), unique(bison_KDE$feeding), unique(bison_KDE$hab))

然后我把它输入ezPredicts来预测争论

ezPredictGLM，to_predict=tp mm%*%f中出错：参数不一致

这就是我目前的处境。我对R比较陌生，在统计方面也不是很有经验，也许有更合理的方法来实现我的目标吗？非常感谢您的帮助

另外，我不能发布数据集，因为它有500000行，而且是机密的。谢谢你抽出时间！

---

/C

你的一部分数据怎么样？你看过lme4中使用预测函数的例子吗？尽管我之前浏览了常见问题解答，但我还没有看到这个例子！似乎正是我想要的，我会试试看。非常感谢。如果我不能让它工作，我会考虑发布一个匿名的数据子集。