||
今天收到加拿大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函数会带来很大的麻烦。至于交互项目,如果是定性变量,目前是无法处理的,如果是定量变量,可以通过交互作用项目的乘积提前设为一个新的变量来解决。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-24 21:29
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社