쿡의 베타 Regrssion 거리

Question

Cook의 거리를 수동으로 계산하고 베타 회귀를 사용하여 cooks.distance 함수를 사용하여 두 가지 다른 결과를 얻었습니다. 누군가가 이유를 이해하도록 도와 줄 수 있습니까?

아래는 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)

여기에 recamended대로 Cook의 거리를 수동으로 계산하는 방법 나는 동일한 문제에도 3불구하고 시도했습니다 2.

당신의 도움에 미리 감사드립니다!

Answer 1

여기서 주요 용도는 잔차 유형입니다. 나는 사용 Response residual했지만 betareg패키지가 사용되었습니다 Pearson residual.