2015-09-28 4 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 (по оси х) и BayesFactors (по оси ординат). Как перевести полученные байесовские факторы из цикла 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

Альтернативный способ, который я нашел (настаивания на цикл) является

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) 

Это также делает трюк.

Смежные вопросы