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 (将#修改为@）