我正在使用R中glmer
包中的lme4
来进行多变量分析。
我的响应变量定义为
rate <- cbind(A, B)
其中A
是同类群组中的阳性病例数,B
是该群组的总数减去A
。
我使用的公式是
glmer(rate ~ C + D + E + F + G + (1 | Z), family=poisson())
其中C
- F
为1或0,G
为计数,Z
为6个等级的因子。
当我运行此模型时,我收到错误updateXwts: dimension mismatch
,但据我所知,glmer
支持多变量分析。
关于如何解决这个问题的任何想法?