Question

我在R中有以下数据集。

df<- data.frame(x = c(rnorm(50, mean=1), rnorm(50, mean=2)), 
                group = rep(1:2, each = 50))

我的目标是使用一些非参数测试来比较两组中x的平均值。除了wilcox.test(x ~ group , data=df, paired= FALSE)之外，我的老板建议我在Stata

中做一些类似于这个命令的事情

mean x, over(group)
 test [x]1 = [x]2

1-如何在R？

中复制类似的结果

2-这项测试做了什么，这项测试是否可以作为Wilcoxon测试的替代非参数测试？

Answer 1

如果这不是您正在寻找的，我道歉。我没有使用Stata。但是，如果您只是尝试采用群组方法，那么使用dplyr包实际上非常简单。

尝试：

library(dplyr)
df %>% 
  group_by(group) %>% 
  summarise(mean.x = mean(x))

返回每组的平均值。

# A tibble: 2 x 2
  group   mean.x
  <int>    <dbl>
1     1 1.155232
2     2 1.778388

如果你想按组计算x的其他统计数据，只需用另一个函数替换mean（x）。