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

我需要用引导包引导一个beta回归模型来检查它的健壮性——因为数据点有很大的库克距离（欢迎其他建议）

我有以下错误：

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

---

t.star[r，]中的

错误问题在于
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

my_beta_reg它工作得很好，谢谢！
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