Distancia de Cook de regresión beta

Question

Calculé la distancia de Cook manualmente y con la función cooks.distance con Beta Regression y obtuve dos resultados diferentes. ¿Alguien puede ayudarme a entender por qué?

A continuación se muestra cómo calculo la (s) distancia (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 se recuerda aquí Cómo calcular manualmente la distancia de Cook Intenté a 3pesar del 2mismo problema.

¡Gracias de antemano por su ayuda!

Answer 1

Aquí el uso principal es el tipo de residuo. Estaba usando Response residualpero el betaregpaquete usó Pearson residual.