0

当我尝试在使用 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)

这表明问题出在随机效应上。当我移除随机效应时,这两个功能都有效。

4

0 回答 0