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

时间:2016-02-28 15:59:50

标签: r bootstrapping robust hypothesis-test

我对比较独立方法的不同稳健方法感到困惑。我在统计教科书中找到了很好的解释。例如yuen(),如果样本量相同。我的样本相当不平等,因此我想尝试一种bootstrap-t方法(来自Wilcox的书:Robust Estimation and Hypothesis Testing简介,第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-function?非常感谢你

1 个答案:

答案 0 :(得分:2)

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

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

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)