我对比较独立方法的不同稳健方法感到困惑。我在统计教科书中找到了很好的解释。例如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?非常感谢你
答案 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)