r语言 - 具有 2 个以上组/级别的稳健 ANCOVA 选项?



我正在尝试使用具有 1 个分组变量/因子的 4 个水平和一个连续协变量进行稳健ANCOVA。

有没有人知道可以执行此功能的功能?

WRS2 的ancova()只需要 2 个级别,当我尝试安装较旧的 WRS 时(如此处 https://stats.stackexchange.com/questions/276697/how-to-implement-a-robust-ancova-with-more-than-2-groups-in-r 所示),我收到以下警告:

安装包中的警告: 软件包"WRS"不可用(对于 R 版本 3.6.0)"。

如果您仍然想与两个以上的组进行健壮的ANCOVA,请按照以下步骤操作:

# load the latest WRS package
> source("https://dornsife.usc.edu/assets/sites/239/docs/Rallfun-v38.txt")
# pairwise comparisons for independent groups with one covariate:
ancmg1(x,y, op=3, tr=0.2)

其中 x 是包含每个组的协变量值的列表,y 是每个组的 DV 值列表,OP 为:

op=1 对修剪均值使用综合检验,由 tr 给出修整。

op=2 对中位数使用综合检验。(当存在平值时不建议使用,使用 op=4)

op=3 使用修剪和百分位数自举的多重比较。此方法似乎最适合一般用途。

OP=4 使用中位数和百分位数自举的多重比较

如果你的数据是长格式的,准备x和y的最简单方法是使用split

x <- split(yourCovariate, list(Group))
y <- split(yourDV, list(Group))