* Stata
ttest y, by(group) unequal       // Welch t 检验（不假定方差齐）
ttest y1 == y2                   // 配对 t 检验
oneway y group, bonferroni       // 单因素 ANOVA + Bonferroni 两两比较
anova y group                    // 方差分析表

## R
t.test(y ~ group, data = df)              // 默认 Welch t 检验
t.test(df$y1, df$y2, paired = TRUE)       // 配对 t 检验
fit <- aov(y ~ group, data = df)
summary(fit); TukeyHSD(fit)               // ANOVA + Tukey 两两比较