r - 在 R-PIRLS step-halvings 中执行 Gamma glmer 时出现错误消息，无法减少 pwrssUpdate 中的偏差

Question

我正在尝试使用 Gamma 错误系列在 R 中执行 glmer。我收到错误消息:“错误:(maxstephalfit) PIRLS 步减半未能减少 pwrssUpdate 中的偏差”

我的响应变量是花质量。我的固定效应是基础质量、F1 处理和施肥方法。我的随机效应是线和嵌套在线内的母亲 ID。

当我使用整数作为响应(即花号)执行相同的分析时，不会发生此错误。

这是我的数据示例:

LINE MATERNAL_ID  F1TREAT  SO   FLWR_MASS  BASE_MASS
17   81           stress   s    2.7514     9.488
5    41           control  o    0.3042     1.809
37   89           control  o    2.3749     6.694
5    41           stress   s    3.6140     9.729
9    5            control  s    0.5020     7.929
13   7            stress   s    0.4914     0.969
35   88           stress   s    0.4418     1.840
1    57           control  o    2.1531     6.673
13   7            stress   s    3.0191     7.131

这是我使用的代码:

library(lme4)
m <- glmer(data=mydata,
           FLWR_MASS~BASE_MASS*F1TREAT*SO+(1 |LINE/MATERNAL_ID),family=Gamma)

(我在 windows 上使用 r 3.0.3)

Answer 1

@HongOoi 在评论中回答了这个问题，但我会在这里为遇到此问题的其他人重复一遍。他建议改变

family=Gamma

到

family=Gamma(link=log)