我正在研究一个包含 44000 名新生儿的数据库，只有两个变量：
gestational age at birth in weeks和bmi. 没有空格。
我运行一个gamlss：
gamlss.males <- gamlss(bmi2clean ~ cs(sapp), family = BCPE, data = bmimales)
一切似乎都很顺利：我得到了百分位曲线和百分位预测。后来我尝试获取参数的edf，如下所示：
edf(gamlss.males, what = c("mu", "sigma", "nu", "tau"), print=TRUE)
或像这样：
edfAll(gamlss.males)
并且总是出现以下错误：
Error in edf[i] <- x[[i]]$edf : replacement has length zero
无法弄清楚为什么。我真的很感谢你的帮助！