Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/64.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 库克';贝塔回归距离_R_Regression_Outliers - Fatal编程技术网

R 库克';贝塔回归距离

R 库克';贝塔回归距离,r,regression,outliers,R,Regression,Outliers,我手动计算了Cook距离,并使用函数cooks.distance和Beta回归计算了Cook距离,得到了两个不同的结果。有人能帮我理解为什么吗 下面是我如何计算库克距离的: require(betareg) df这里的主要问题是残余类型。我使用的是Response resistance,但是betareg使用的包Pearson resistance看看函数本身:betareg:::cooks.distance.betareg,它只有几行代码 require(betareg) df<-dat

我手动计算了Cook距离,并使用函数cooks.distance和Beta回归计算了Cook距离,得到了两个不同的结果。有人能帮我理解为什么吗

下面是我如何计算库克距离的:

require(betareg)

df这里的主要问题是残余类型。我使用的是
Response resistance
,但是
betareg
使用的包
Pearson resistance

看看函数本身:
betareg:::cooks.distance.betareg
,它只有几行代码
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)