今天收到加拿大Sherbrooke大学一个博士生的来写咨询为什么模型里面带有scale标准化处理固定效应项目和交互项无法使用glmm.hp包，原文如下：

她咨询的glmm模型是用glmmTMB来处理：Numeric1 ~ Categoric1 + scale(Numeric2) + scale (Numeric3) + scale (Numeric4) + Categoric3 + Categoric4 + Categoric3*Categoric4  + (1 | Random1) + (1 |Random2) , family = "poisson",data = SubsetX

当前的版本glmm.hp包，不能在模型里面出现scale函数，也不能有交互作用项目。如果把模型改为

Numeric1 ~ Categoric1 + Numeric2 + Numeric3 + Numeric4 + Categoric3 + Categoric4  + (1 | Random1) + (1 |Random2) , family = "poisson",data = SubsetX

这样应该是可以运行的。如果要标准化变量，可以通过提前处理来解决，因为glmm.hp是从原始模型提取文本变量，scale函数会带来很大的麻烦。至于交互项目，如果是定性变量，目前是无法处理的，如果是定量变量，可以通过交互作用项目的乘积提前设为一个新的变量来解决。

转载本文请联系原作者获取授权，同时请注明本文来自赖江山科学网博客。

链接地址：https://wap.sciencenet.cn/blog-267448-1374989.html?mobile=1

收藏

分享