Mann-whitney tesy 错误：组因子 2 水平？答案

【问题标题】：Mann-whitney tesy error: group factor 2 levels?Mann-whitney tesy 错误：组因子 2 水平？
【发布时间】：2018-07-18 09:30:19
【问题描述】：

我不断收到以下错误：

y <- as.factor(data$Pos.44)
x <- as.factor(data$Neg.44)
wilcox.test( x ~ y,  
                 alternative='two.sided', 
                 paired=FALSE, var.equal = TRUE, 
               conf.level = 0.95)

wilcox.test.formula 中的错误（x ~ y，alternative = "two.lateral"，paired = FALSE，：分组因子必须正好有 2 个级别

我尝试将数据转换为因子，如数字，并告诉 r 标题为 TRUE。有任何想法吗？？

这是一些数据。

Neg 44      Pos 44
13.8228455  7.126694925
13.8228455  8.402648457
13.8228455  8.402648457
6.85288518  12.27621627
9.232427575 
11.31157309 
195.5120218 
82.61317544 
9.751455813 
4.490663867 
7.732335943 
3.72520288  
17.51041865 
8.56912795  
17.59120857 
12.3571618

【问题讨论】：

您能否以可复制的“复制和粘贴”格式共享数据？使用dput 或代码生成示例数据。还可以使用编辑工具正确格式化代码、数据和文本。
据我所知，你不需要把你的数值变量变成因子，你需要的是Reshaping data.frame from wide to long format。

标签： r

【解决方案1】：

您的格式错误。 Mann-Whitney 检验的输入需要一列包含值，另一列包含两水平因子。您尝试向其提供两个转换为因子的数值向量。请看下面的例子：

x=c(rnorm(50,0,1),rnorm(50,1,1))
y=rep(c('A','B'),c(50,50))
wilcox.test(x~y, alternative='two.sided',  paired=FALSE, var.equal = TRUE,  conf.level = 0.95)

【讨论】：

谢谢。那么剧本会是什么？ data

【解决方案2】：

你可以试试

library(tidyverse)
d %>% 
  gather(k,v) %>% 
  with(.,wilcox.test(v ~ k, data=., exact =F))

Wilcoxon rank sum test with continuity correction

data:  v by k
W = 45, p-value = 0.2367
alternative hypothesis: true location shift is not equal to 0

或使用箱线图可视化您的数据

library(ggbeeswarm)
library(ggsignif)
d %>% 
  gather(k,v) %>% 
  ggplot(aes(k,v)) + 
   geom_boxplot() +
   geom_beeswarm() + 
   ggsignif::geom_signif(comparisons = list(c("Neg44", "Pos44")))

您的数据

d <- read.table(text="Neg44 Pos44
13.8228455  7.126694925
13.8228455  8.402648457
13.8228455  8.402648457
6.85288518  12.27621627
9.232427575 
11.31157309 
195.5120218 
82.61317544 
9.751455813 
4.490663867 
7.732335943 
3.72520288  
17.51041865 
8.56912795  
17.59120857 
12.3571618",header=T,  fill=T)

【讨论】：