||
glmm.hp函数输入对象很非常简单,即lme4包和nlme包运行(G)LMM得到对象均可。但lme4包允许输入的数据框带有NA值,但计算过程自动会把带有NA的行全部去掉,但输出的对象的data还是原始带有NA的data,(而nlme是从一开始就不允许带有NA值的数据输入)。glmm.hp是从lme4输出对象获取原始的data, 这个data里面如果带有NA,就会导致在计算各种解释变量的组合的情况下所使用数据的总行数不一致的情况,那最终平均分配共享R2显然是不对的。现在已经修正,在计算模型组合之前先把带有NA的行去掉,这样就保证每种情况下数据均为一致。
请大家更新包,CRAN和github均已经更新,如果您之前的数据解释变量(固定效应)带有NA,且是用lme4包计算的,请您更新包后重新计算。
R包总是在用户不断尝试反馈过程中完善,这是R开源的优势和魅力所在!
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 22:21
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社