R 组间比较Mann-Whitney检验
在我的数据集中,我有以下列: 组哪个是分组变量R 组间比较Mann-Whitney检验,r,dplyr,data.table,R,Dplyr,Data.table,在我的数据集中,我有以下列: 组哪个是分组变量 var哪个是范围变量 我只需要比较1-2组和1-3组。 如何使用wilcox.test进行测试? 这里是数据 mydat=structure(list(group = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L), var = c(23L, 24L, 24L, 23L, 23L, 24L, 24L, 23L, 23L, 24L, 24L, 23L)), .Names = c("group",
var哪个是范围变量 我只需要比较1-2组和1-3组。 如何使用wilcox.test进行测试?
这里是数据
mydat=structure(list(group = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L,
3L, 3L, 3L), var = c(23L, 24L, 24L, 23L, 23L, 24L, 24L, 23L,
23L, 24L, 24L, 23L)), .Names = c("group", "var"), class = "data.frame", row.names = c(NA,
-12L))
您可以通过
group
对数据进行split
frame,然后lappy
对所需的组对进行测试
lst <- split(mydat, mydat$group)
lapply(lst[-1], function(DF) wilcox.test(DF$var, lst[[1]]$var, exact = FALSE))
lst您可以split
data.frame按group
进行分组,然后lappy
对所需的组对进行测试
lst <- split(mydat, mydat$group)
lapply(lst[-1], function(DF) wilcox.test(DF$var, lst[[1]]$var, exact = FALSE))
lst您也可以尝试一种可视化解决方案
library(tidyverse)
library(ggsignif)
ggplot(mydat, aes(factor(group), var)) +
geom_boxplot() +
geom_signif(comparisons = list(c(1,3), c(1,2)), test = "wilcox.test",
test.args = list(exact = FALSE),step_increase = 0.1)
您也可以尝试视觉解决方案
library(tidyverse)
library(ggsignif)
ggplot(mydat, aes(factor(group), var)) +
geom_boxplot() +
geom_signif(comparisons = list(c(1,3), c(1,2)), test = "wilcox.test",
test.args = list(exact = FALSE),step_increase = 0.1)
很抱歉,但也许你能在这篇相关的文章中帮助我?对不起,也许你能帮我写这篇相关的文章?