首页 > 解决方案 > 从硬币包中的对象 QuadTypeIndependenceTest 和 ScalarIndependenceTest 中提取 p 值

问题描述

使用Aids2package 中的数据集MASS,我正在应用 Ansari-Bradley Non-Parametric Test 通过此片段测试组独立性

library(coin)
library(MASS)
a1 <- ansari_test(Aids2$diag~Aids2$state)
a2 <- ansari_test(Aids2$diag~Aids2$sex)

如果我看一下 a1 和 a2 对象,它会给出如下摘要:

Asymptotic K-Sample Ansari-Bradley Test

data:  Aids2$diag by Aids2$state (NSW, Other, QLD, VIC)
chi-squared = 4.4829, df = 3, p-value = 0.2138

Asymptotic Two-Sample Ansari-Bradley Test

data:  Aids2$diag by Aids2$sex (F, M)
Z = -1.6765, p-value = 0.09364
alternative hypothesis: true ratio of scales is not equal to 1

与相应的对象类是:

> class(a1)
[1] "QuadTypeIndependenceTest"
attr(,"package")
[1] "coin"
> class(a2)
[1] "ScalarIndependenceTest"
attr(,"package")
[1] "coin"

但是,由于我只对 p-value (p-value = 0.2138) 感兴趣,我只想从对象中提取它,但我找不到如何提取它..

标签: rp-valuestatistical-test

解决方案


由于对象喜欢"QuadTypeIndependenceTest""ScalarIndependenceTest"是从coin包的结果中创建的,因此有特定的函数来提取 pvalue,使用coin::pvalue(obj),特别感谢您指出@AntoniosK


推荐阅读