Cook's Distance of Beta Regrssion

Aug 23 2020

Ho calcolato la distanza di Cook manualmente e con la funzione cooks.distance con Beta Regression e ho ottenuto due risultati differenti. Qualcuno può aiutarmi a capire perché?

Di seguito è riportato come calcolo le distanze di 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)

Come consigliato qui Come calcolare manualmente la distanza di Cook Ho provato 3nonostante lo 2stesso problema.

Grazie in anticipo per il vostro aiuto!

Risposte

1 Angel Aug 24 2020 at 10:06

Qui il problema principale è il tipo di residuo. Stavo usando Response residualma betaregpacchetto usato Pearson residual.