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

Question

我遇到了glmer模型的一些问题。

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

我用

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

descdist(mydata$y)

从我上传的图片中我假设我的dependend变量作为伽玛分布分发。

我也试过

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中看到。

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