进行两个样本t检验，并输出R中两组矩阵的t、p值

我正试图进行这样的两样本t检验

wt1

，

wt2

，

wt3

，

mut1

，

mut2

，

mut3

是3x3矩阵。运行t检验后，我想得到

t.stat

和

p.value

矩阵，其中

t.stat[i,j] <- the t value from t.test(c(wt1[i,j],wt2[i,j],wt3[i,j]),c(mut1[i,j],mut2[i,j],mut3[i,j]))
p.value[i,j] <- the p-value from t.test(c(wt1[i,j],wt2[i,j],wt3[i,j]),c(mut1[i,j],mut2[i,j],mut3[i,j]))

---

t.stat[i，j]您可以这样做：

test<- sapply(1:9, function(x) t.test(c(wt1[x], wt2[x], wt3[x]), 
                                      c(mut1[x], mut2[x], mut3[x])))
t.stat<- matrix(test["statistic", ], nrow = 3)
p.value<- matrix(test["p.value", ], nrow = 3)

您可以这样做：

test<- sapply(1:9, function(x) t.test(c(wt1[x], wt2[x], wt3[x]), 
                                      c(mut1[x], mut2[x], mut3[x])))
t.stat<- matrix(test["statistic", ], nrow = 3)
p.value<- matrix(test["p.value", ], nrow = 3)

差不多。您还可以将
sapply
+
t.test
的结果保存一次，然后提取零件<代码>lappy（结果，`[`，c（“统计”，“p.value”））

或类似内容。@latemail非常好的建议！我相应地编辑了答案！非常感谢Yannis的帮助！我更新了问题，你能再帮我一次吗？@latemail非常好。你也可以将

sapply

+

t.test的结果保存一次，然后提取部分。
lappy（results`[`，c（“statistic”，“p.value”）
或类似内容。@thelatemail非常好的建议！我据此编辑了答案！非常感谢Yannis的帮助！我更新了问题，你能再次帮助我吗？@thelatemail