如何複製塔塔的沃爾德檢驗（後評估）中的R

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）。