标签: r statistics
我希望以比单个t检验更有效的方式对所有因子水平之间的差异进行t检验(实际数据有50个因子水平)。
示例数据:
dat <- cbind(rnorm(10), c("x","x","x","y","y","y","z","z","a","a"))
我希望输出均值差异迭代的t-stat(x-y,y-z,z-x,a-x,a-y,a-z)
我已经看到如何使用申请差异0的测试来做到这一点,但无法联系如何利用这种方法来区分手段。可以通过for循环破解它,但必须有更优雅的方式。