我正在尝试使用R分析某些行为数据。我的因变量是一个非常适合对数正态分布的连续变量吗?对这些数据执行GLMM的理想方法应该是什么?
我已经使用以下代码尝试过lme4
软件包。我不知道这是否正确(在某个地方的博客中阅读)
mod1<-glmer (time~ A+B+A*B+ (1|Fish), data=df1, family = Gamma(link = "log"))
固定效果模型矩阵的秩不足,因此删除12列/系数
警告信息: 在checkConv(attr(opt,“ derivs”),opt $ par,ctrl = control $ checkConv,中: 模型无法与max | grad |收敛= 0.0160032(tol = 0.001,第1部分)
ss <- getME(mod1,c("theta","fixef"))
mod1.1 <- update(mod1,start=ss)
我尝试在glmmPQL
包中使用MASS
,但出现错误:
迭代1 getGroups.data.frame(dataMix,groups)中的错误: 组的无效公式