使用 R，我有兩個向量 a 和 b：
a <- rnorm(3, mean = 6, sd = 2)
b <- rnorm(3, mean = 8, sd = 2)
我想生成a和b指定次數，每次執行t.test比較平均值a和平均值b
t.test(a, b, paired = FALSE)
我如何編碼重復rnorm和t.test函式指定的次數并檢索完整的輸出t.test()或只是p.value
uj5u.com熱心網友回復：
要擴展您的評論：
pvals <- replicate(20, t.test(rnorm(3, mean = 6, sd = 2), rnorm(3, mean = 8, sd = 2))$p.value)
sum(pvals < .05)
# [1] 2
hist(pvals)
獲得所有結果
results <- replicate(20, t.test(rnorm(3, mean = 6, sd = 2), rnorm(3, mean = 8, sd = 2)), simplify=FALSE)
results[[1]]
#
# Welch Two Sample t-test
#
# data: rnorm(3, mean = 6, sd = 2) and rnorm(3, mean = 8, sd = 2)
# t = 0.052882, df = 2.6225, p-value = 0.9616
# alternative hypothesis: true difference in means is not equal to 0
# 95 percent confidence interval:
# -8.880021 9.155836
# sample estimates:
# mean of x mean of y
# 7.156974 7.019067
其他結果是 results[[2]], . . . ，結果[[20]]。需要該simplify=FALSE引數來保留每次運行的串列結構。
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/424468.html
上一篇：使用文本檔案中的每個串列執行命令