如何解决R中多个组的威尔科克森测试错误
我想测试五个不同组的多样性的统计差异。
我的数据如下所示:
样品名称 | 治疗 | 巢 | 观察到 |
---|---|---|---|
a1 | 控制 | 1 | 36 |
b2 | 啶虫脒 | 1 | 25 |
c3 | 噻虫啉 | 1 | 36 |
d4 | 草酸 | 1 | 36 |
e5 | 控制 | 1 | 36 |
我想进行 Wilcoxon 检验,我尝试了两件事:
foo <- pairwise.wilcox.test(data$Observed,data$Treatment,p.adjust.method="bonferroni")
foo
当我运行第一行代码时,出现此错误:
在 wilcox.test.default(xi,xj,paired = paired,...) 中:不能用关系计算精确的 p 值
但是当我打印 foo 对象时,我得到了组之间每个比较的 p 值,我应该相信这些结果吗?
我尝试过的另一件事是:
data %>%
filter(Treatment != "Control") %>%
group_by(Treatment) %>%
summarise(p_value = wilcox.test(Data$Observed,Observed,exact = FALSE)$p.value,p.adjust.methods="bonferroni")
在这种情况下,它有效,但 bonferroni 校正无效。
您认为哪种方式最好?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。