使用 {betareg} 在 R 中引导 Beta 回归模型时出错

Question

我需要 bootstrap 一个 beta 回归模型来检查它的稳健性——因为一个数据点有一个很大的厨师距离——和引导包（欢迎其他建议）。

我有以下错误：

Error in t.star[r, ] <- res[[r]] : 
  incorrect number of subscripts on matrix

这是一个可重现的例子：

library(betareg)
library(boot)

fake_data <- data.frame(diet = as.factor(c(rep("A",10),rep("B",10))),
                            fat = c(runif(10,.1,.5),runif(10,.4,.9)) )
    
plot(fat~diet, data = fake_data)
   
my_beta_reg <- function(data,i){
        data_i <- data[i,]
        mod <- betareg(data_i[,"fat"] ~ data_i[,"diet"])
        return(mod$coef)
      }
    
b = boot(fake_data, statistic = my_beta_reg, R= 50)

Error in t.star[r, ] <- res[[r]] : 
incorrect number of subscripts on matrix

有什么问题？
提前致谢。

Answer 1

问题是 mod$coef 是一个列表：

betareg(fat ~ diet, data = fake_data)$coef
#$mean
#(Intercept)       dietB 
#  -1.275793    2.490126 
#
#$precision
#   (phi) 
#20.59014

您需要 unlist 它或者最好使用您应该用于提取系数的函数：

my_beta_reg <- function(data,i){
  mod <- betareg(fat ~ diet, data = data[i,])
  #unlist(mod$coef)
  coef(mod)
}

b = boot(fake_data, statistic = my_beta_reg, R= 50)
print(b)
#ORDINARY NONPARAMETRIC BOOTSTRAP
#
#
#Call:
#boot(data = fake_data, statistic = my_beta_reg, R = 50)
#
#
#Bootstrap Statistics :
#     original       bias    std. error
#t1* -1.275793 -0.019847377   0.2003523
#t2*  2.490126  0.009008892   0.2314521
#t3* 20.590142  8.265394485  17.2271497

使用 {betareg} 在 R 中引导 Beta 回归模型时出错

Error when Bootstraping a Beta regression model in R with {betareg}

regression

r