当我尝试在使用 GAMLSS 拟合的具有随机效应的模型上使用挖泥机时,我收到以下错误。
data(usair)
usair$rand<-c(rep("A",11),rep("B",10),rep("C",10),rep("D",10))
fit1<-gamlss(y~x1+x2+x3+x4+x5+re(random=~1|rand), data=usair, family = "GA",na.action="na.fail" )
dredge(fit1, rank = "BIC")
矩阵错误(NA_real_,ncol = 3L,nrow = length(x),dimnames = list(coefNames,:'dimnames' [1] 的长度不等于数组范围
谁能告诉我为什么?有办法解决吗?
此外,当我尝试在模型上使用函数 ggpredict 时,我收到以下错误:
错误:无法计算置信区间。
- 原因:“可变长度不同(为 're(random = ~1 | rand)' 找到)”
- 来源:model.frame.default(object, data, xlev = xlev)
这表明问题出在随机效应上。当我移除随机效应时,这两个功能都有效。