3

我的glmer模型有一些问题。

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

我检查我的因变量(y)的分布

descdist(mydata$y)

decdist_image

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

我也试过

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

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

Fit_gamma

所以我尝试了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 中可以看出,我所有的值都是正数。

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

4

0 回答 0