2015-09-28 54 views
0

我是R新手,正在嘗試基於for循環的結果創建一個繪圖。我現在用的是BayesFactor包找幾個貝葉斯因子R - BayesFactor:從for循環創建一個圖

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循環轉換爲矢量?

謝謝!

回答

0

喜歡的東西:

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)) 
0

,我已經找到一個替代方法(堅持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) 

這也做的伎倆。