我的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 中可以看出，我所有的值都是正数。
所以，我的问题是为什么我会收到该消息，或者我做错了什么？