我正在研究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
无法弄清楚原因。 我真的很感谢你的帮助！