У меня есть некоторые проблемы с моделью glmer.
Мой набор данных состоит из непрерывной зависимой переменной (y), которая изменяется от 200 до 2000. Мера берется в миллисекундах (мс).У меня есть три предиктора, 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.
Итак, мой вопрос: почему я получаю это сообщение или что я делаю неправильно?