Distância de Cook de Beta Regrssion

Question

Calculei a distância de Cook manualmente e com a função cooks.distance com Beta Regression e obtive dois resultados diferentes. Alguém pode me ajudar a entender por quê?

Abaixo está como eu calculo a (s) distância (s) de Cook:

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)

Como recomendado aqui. Como calcular manualmente a distância de Cook, tentei 3apesar do 2mesmo problema.

Agradeço antecipadamente por sua ajuda!

Answer 1

Aqui, o uso principal é um tipo de resíduo. Eu estava usando, Response residualmas o betaregpacote foi usado Pearson residual.