使用数据分层进行wilcoxon检验

Question

我有一个非常基本的问题。我将一种化学品的浓度存储在第一列中，将研究参与者的性别存储在第二列中。 wilcoxon测试的代码是什么，以查看男孩中发现的浓度与女孩中发现的浓度之间是否存在差异？对代码的一些解释对我了解其工作原理也很有用。谢谢！

我得到了用于ANOVA测试的代码，这也很好。谁能告诉我它是否满足我的要求？

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方法        wilcox.test（公式，数据，子集，na.action等）
  ...
  公式：格式为“ lhs〜rhs”的公式，其中“ lhs”是数字             变量赋予数据值和“ rhs”一个因子，其中有两个             给出相应组的级别。

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