Bootstrap-t 在 R 中比较修剪均值的方法

Bootstrap-t Method for Comparing Trimmed Means in R

我对比较独立均值的不同稳健方法感到困惑。我在统计教科书中找到了很好的解释。例如 yuen() 在样本大小相等的情况下。我的样本相当不平等,因此我想尝试 bootstrap-t 方法(来自 Wilcox 书:稳健估计和假设检验简介,第 163 页)。它说 yuenbt() 是一个可能的解决方案。

但是所有教科书都说我可以在这里使用向量:

yuenbt(x,y,tr=0.2,alpha=0.05,nboot=599,side=F)

如果我查看本地描述,它会说:

yuenbt(formula, data, tr = 0.2, nboot = 599)

我的试用有什么问题:

x <- c(1,2,3)
y <- c(5,6,12,30,2,2,3,65)
yuenbt(x,y)

为什么我不能对我的两个向量使用 yuenbt 函数?非常感谢

查看 yuenbt 的帮助(对于那些想知道的人,yuenbt 来自包 WRS2...),它需要一个公式和一个数据框作为参数。我的印象是它需要长格式的数据。

使用您的示例数据,我们可以这样实现:

library(WRS2)

x <- c(1,2,3)
y <- c(5,6,12,30,2,2,3,65)

dat <- data.frame(value=c(x,y),group=rep(c("x","y"), c(length(x),length(y))))

然后我们可以使用函数:

yuenbt(value~group, data=dat)