我已经在 R 中运行了 beta 回归,并希望评估进行残差诊断。我使用了 plot 函数并获得了图,但是,没有标记潜在的异常值。怎么才能给异常值对应的标签呢?
breg.full <- betareg(Percentage ~ Total_testscore + Campus + Programme +
Gender + SE_track + Hours_Math_SE, data = starters, # [-c(53, 24, 35), ]
link = "logit") # , , link.phi = NULL, type = "ML"
summary(breg.full)
par(mfrow = c(2,3))
plot(breg.full, which = 1:6)
编辑:
作者提供了此代码的链接(http://www.de.ufpe.br/~cribari/betareg_example.zip。)但是它不再工作......