使用lme4 / MASS

时间:2019-04-04 06:41:55

标签: r glm lme4

我正在尝试使用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)中的错误:     组的无效公式

0 个答案:

没有答案