r - 'gamma' 系列不允许使用非正值

Question

我的glmer模型有一些问题。

我的数据集由一个从 200 到 2000 的连续因变量 (y) 组成。测量以毫秒 (ms) 为单位。我有三个预测变量，A、B、C。前两个预测变量有两个级别，A1、A2、B1、B2，最后一个有三个级别 C1、C2 和 C3。

我检查我的因变量（y）的分布

descdist(mydata$y)

从我上传的图像中，我假设我的因变量是作为伽马分布分布的。

我也试过

fit.gamma <- fitdist(mydata$RT, distr = "gamma", method = "mme")

当我绘制拟合时，它似乎很合适

所以我尝试了glmer一个family=Gamma：

model<-glmer(y~A+B+C+A*B*+(1|Subject), family= Gamma(), data=mydata)

但我收到错误消息：

Error in eval(family$initialize, rho) : 
  non-positive values not allowed for the 'gamma' family

从图像 Fit_gamma 中可以看出，我所有的值都是正数。

所以，我的问题是为什么我会收到该消息，或者我做错了什么？

标签： r