数据分层的wilcoxon检验

Question

我有一个非常基本的问题。我把一种化学物质的浓度储存在一栏中，而研究参与者的性别则储存在第二栏中。做wilcoxon检验的代码是什么，看看男孩的浓度和女孩的浓度是否有差异？对代码的一些解释对于我理解它是如何工作的也是有用的。谢谢!

我得到了这个方差分析测试的代码，这也很好。有人能告诉我它是否做了我需要的事吗？

av <- aov(UC_MEHP ~ BQF05C1, data=data)
av
summary(av)

输出如下所示

> av <- aov(UC_MEHP ~ BQF05C1, data=data)
> av
Call:
   aov(formula = UC_MEHP ~ BQF05C1, data = data)

Terms:
                  BQF05C1 Residuals
Sum of Squares     0.3445 2917.4564
Deg. of Freedom         1       151

Residual standard error: 4.395555
Estimated effects may be unbalanced
21 observations deleted due to missingness
> summary(av)
             Df Sum Sq Mean Sq F value Pr(>F)
BQF05C1       1    0.3   0.344   0.018  0.894
Residuals   151 2917.5  19.321               
21 observations deleted due to missingness

抱歉，我知道这不是一个很高级的问题.

Answer 1

来自?wilcox.test

.‘>类的S3方法(公式，数据，子集，na.action，.) 公式:形式‘lhs~rhs’的公式，其中‘lhs’是给出数据值的数值变量，‘rhs’是给出相应组的两个级别的因子。

因此，wilcox.test(UC_MEHP ~ BQF05C1, data=data)应该有效(假设BQF05C1是指定性别的列，UC_MEHP是集中的列)。