我正在尝试使用具有 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))