【问题标题】:Robust ANCOVA in R with package WRS2R 中的鲁棒 ANCOVA,带有 WRS2 包
【发布时间】:2025-11-22 13:05:02
【问题描述】:

WRS 包中的函数 ancboot 使用 20% 修剪均值和自举执行单向 ANCOVA,其中一个协变量:

ancboot(formula, data, tr = 0.2, nboot = 599, sm = FALSE, fr1 = 1, fr2 = 1,
        pr = TRUE, pts = NA)

该命令返回一个 ancova 类的对象,其中包含 5 个设计点处的修剪均值:

[1] "Note: confidence intervals are adjusted to control FWE"
[1] "But p-values are not adjusted to control FWE"
[1] "Taking bootstrap samples. Please wait."
$output
X n1 n2 DIF TEST ci.low ci.hi p.value
[1,] 2 15 30 -1.277778 -1.846494 -3.234577 0.6790212 0.0855
[2,] 4 38 50 -1.066667 -2.139515 -2.476449 0.3431158 0.0440
[3,] 4 38 50 -1.066667 -2.139515 -2.476449 0.3431158 0.0335
[4,] 5 42 50 -1.092308 -2.256915 -2.460882 0.2762670 0.0320
[5,] 6 25 41 -1.333333 -2.341798 -2.943341 0.2766743 0.0300
$crit
[1] 2.8277

但是,我曾期望报告稳健的 ANCOVA,很像传统(参数)ANCOVA,包括综合 p 值、95% CI、F 和 [ df]。如何使用稳健的 ANCOVA(即,使用修整均值和自举)获得所需的输出?

【问题讨论】:

  • WRS 还是 WRS2 包?
  • @Roman Luštrik 我刚刚仔细检查过,我实际上使用的是 WRS 版本 0.27.5。

标签: r


【解决方案1】:

我检查了raw source file 是否有任何方法的痕迹,但结果似乎没有特殊的类,这意味着不能定义特殊的方法。要求包作者/维护者包含此功能或提供源代码。

【讨论】: