赖江山的博客分享 http://blog.sciencenet.cn/u/laijiangshan 生态、统计与R语言

博文

glmm.hp包不允许带有scale函数和交互作用项(回答一个加拿大博士生的问题)

已有 2175 次阅读 2023-2-5 22:28 |个人分类:glmm|系统分类:科研笔记

今天收到加拿大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函数会带来很大的麻烦。至于交互项目,如果是定性变量,目前是无法处理的,如果是定量变量,可以通过交互作用项目的乘积提前设为一个新的变量来解决。

image.png



https://blog.sciencenet.cn/blog-267448-1374989.html

上一篇:rdacca.hp文章引用到了100次(researchgate统计数据)
下一篇:更新rdacca.hp包到1.0-9版本,升级permu.hp函数,让tibble数据类型也能使用
收藏 IP: 222.129.1.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-24 21:29

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部