Cook 的 Beta 回归距离

Question

我手动计算了库克的距离，并使用函数 cooks.distance 和 Beta 回归，我得到了两个不同的结果。有人可以帮我理解为什么吗？

以下是我计算库克距离的方法：

require(betareg)
df<-data("ReadingSkills")
y<-ReadingSkills$accuracy
n<-length(y)

bfit<-betareg(accuracy ~ dyslexia + iq, data = ReadingSkills)
yhat<-fitted(bfit)
cook<-cooks.distance(bfit)
hatv<-hatvalues(bfit)
res<-residuals(bfit, type = "response")

RHS<-hatv/(1-hatv)

s23<-var(res)###
p23<-res^2/s23
Response<-(1/2)*p23*RHS

cbind(Response, cook)

正如这里推荐的How to manually calculate Cook's distance 我尝试了3 尽管2 但同样的问题。

提前感谢您的帮助！

Answer 1

这里的主要问题是残差类型。我使用的是Response residual，但使用了betareg 包Pearson residual。