クックのベータ登録の距離

Question

クックの距離を手動で計算し、ベータ回帰を使用して関数cooks.distanceを使用すると、2つの異なる結果が得られました。誰かが私が理由を理解するのを手伝ってくれませんか？

以下は、クックの距離を計算する方法です。

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)

ここで推奨されているように、クックの距離を手動で計算する方法3は2、同じ問題にもかかわらず試しました。

よろしくお願いします！

Answer 1

ここでの主な問題は、残差のタイプです。使用してResponse residualいましたが、betaregパッケージを使用しましたPearson residual。