为什么下面创建的 gamlss() 模型预测超出响应值范围？我是否错误地使用了 predict.gamlss() 函数，或者模型只是那么差？
我的响应变量是具有零和一个的比例数据，所以我使用的是 gamlss() 内的 BEINF、零和一个膨胀的 beta 分布。下面是一个数据样本，但是这个问题出现在完整的数据集上，并且也出现在比下面的更复杂的 gamlss() 模型形式中。
library(gamlss)
data<-data.frame(response=c(0.33, 0.43, 0, 0, 0.90, 0.76, 0.33, 0.76, 0.26, 0.80, 0.33, 0), predictor=c(14.9, 541.5, -80.4, -63.9, 1082.5, 709.7, 368.8, 376.7, 107.0, 885.8, 11.0, -33.5))
mod<-gamlss(formula=response~predictor, nu.formula=~predictor, tau.formula=~predictor, family=BEINF, data=data)
mod.pred<-predict(mod, newdata=NULL)
plot(mod.pred, data$response, xlab="predicted", ylab="actual", main = "predicted vs actual values for gamlss model")
谢谢