如何为t-test压缩多行代码？

Question

我正在处理一个简单的数据集，它最好使用t- best进行分析，但至少有24个t-best，这会创建很多代码行和许多表！我想知道是否有一种方法可以让一行代码包含所有的y变量，因为我的分组变量保持不变，然后它们显示在一个表中。谢谢你的帮助。另外，如果有一种简单的方法可以将结果以表格格式导出到word或R，请让我知道。我会很感激的。

以下是我的代码示例。

t.test(XC[tissue =="Ent"], XC[tissue =="Stom"])
t.test(XN[tissue =="Ent"], XN[tissue =="Stom"])
t.test(CN[tissue =="Ent"], CN[tissue =="Stom"])

Answer 1

您可以使用apply函数

dta=data.frame(XC,XN,CN,tissue)
res=data.frame(apply(dta[,1:3],2,function(z)(unlist(t.test(z~dta$tissue))[1:8])),stringsAsFactors = F)
res[]=sapply(res,as.numeric)
round(res,6)

您将拥有类似以下内容的内容：

# > round(res,6)
#                                      XC         XN         CN
# statistic.t                    -109.63424 -10.842582 -45.532057
# parameter.df                    101.66755 105.801883 103.753688
# p.value                           0.00000   0.000000   0.000000
# conf.int1                       -20.25112  -2.304988 -10.219741
# conf.int2                       -19.53135  -1.592335  -9.366679
# estimate.mean in group Ent       20.16557  10.044373  20.172782
# estimate.mean in group Stom      40.05680  11.993035  29.965992
# null.value.difference in means    0.00000   0.000000   0.000000