Cooks Entfernung von Beta Regrssion

Question

Ich habe die Entfernung von Cook manuell und mit der Funktion cooks.distance with Beta Regression berechnet und zwei unterschiedliche Ergebnisse erhalten. Kann mir bitte jemand helfen zu verstehen warum?

Im Folgenden wird die Entfernung (en) von Cook berechnet:

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)

Wie hier empfohlen. Wie man Cooks Entfernung manuell berechnet, habe ich versucht, 3trotz des 2gleichen Problems.

Vielen Dank im Voraus für deine Hilfe!

Answer 1

Hier ist die Hauptausgabe die Art des Residuums. Ich habe Response residualaber betaregPaket benutzt Pearson residual.