R - BayesFactor：从 for 循环创建图答案

【问题标题】：R - BayesFactor: Creating a plot from a for loopR - BayesFactor：从 for 循环创建图
【发布时间】：2015-09-28 13:58:06
【问题描述】：

我是 R 新手，我正在尝试根据 for 循环的结果创建一个绘图。我正在使用 BayesFactor 包来查找几个贝叶斯因子

N <- seq (10, 500, by = 5)

for (i in 1:length(N)) {ttest.tstat(2.24, N[i], simple = TRUE}

现在我想用 N（在 x 轴上）和 BayesFactors（在 y 轴上）创建一个基本图。如何将 for 循环中生成的贝叶斯因子转换为向量？

谢谢！

【问题讨论】：

【解决方案1】：

类似：

N <- seq (10, 500, by = 5)
tstat_vec = sapply(N, ttest.stat, t = 2.24, simple = TRUE)
plot(N, tstat_vec)

应该可以解决问题。

或者使用ggplot2:

df = data.frame(N = N, tstat = tstat_vec)
ggplot(df) + geom_point(aes(x = N, y = tstat))

【讨论】：

【解决方案2】：

我发现的另一种方法（坚持 for 循环）是

N <- seq(10, 500, by=5)
BF <- rep(0, 99)

for (i in 1:length(N)) {BF[i]=ttest.tstat(2.24, N[i], simple=TRUE)[[1]]}

plot.default(N,BF)

这也可以解决问题。

【讨论】：