3

我正在尝试使用 Gamma 错误系列在 R 中执行 glmer。我收到错误消息：“错误：(maxstephalfit) PIRLS step-halvings failed to reduce deviance in 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）

4

1 回答 1

1

@HongOoi 在评论中回答了这个问题，但我会在这里重复给其他有这个问题的人。他建议改变

family=Gamma

至

family=Gamma(link=log)
于 2014-08-27T02:05:15.373 回答