如何在添加参与者时绘制顺序贝叶斯因子

Question

我目前正在使用顺序贝叶斯因子方法分析眼球追踪数据，我想绘制生成的贝叶斯因子（BF；根据平均观看时间计算）如何随着参与者的添加而变化。

我希望 x 轴代表参与计算的人数，y 轴代表生成的贝叶斯因子。

例如，当包括参与者 1-10 时，BF = [y 值]，即图表上的一个绘图点。当参与者 1-11 被包括在内时，BF = [y 值 ]，这是图表上的第二个绘图点。

有没有办法在 R 中做到这一点？

比如我有这个数据集：

ID      avg_PTL
 <chr>     <dbl>
 1 D07   -0.0609  
 2 D08    0.0427  
 3 D12    0.112   
 4 D15   -0.106   
 5 D16    0.199   
 6 D19    0.0677  
 7 D20    0.0459  
 8 d21   -0.158   
 9 D23    0.0650  
10 D25    0.0579  
11 D27    0.0463  
12 D29    0.00822 
13 D30    0.00613 
14 D36   -0.0484  
15 D37    0.0312  
16 D39    0.000547
17 D44    0.0336  
18 D46    0.0514  
19 D48    0.236   
20 D51   -0.000487
21 D60    0.0410  
22 D61    0.0622  
23 D62    0.0337  
24 D64   -0.125   
25 D65    0.215   
26 D66    0.200

我用以下方法计算 BF：

bf.mono.correct = ttestBF(x = avg_PTL_mono_correct$avg_PTL)

非常感谢任何提示！

Answer 1

您可以多次使用 sapply 到运行测试，每次只对观察向量进行子集化。例如

srange <- 10:nrow(avg_PTL_mono_correct)
BF <- sapply(srange, function(i) {
  extractBF(ttestBF(x = avg_PTL_mono_correct$avg_PTL[1:i]), onlybf=TRUE)
})

plot(srange, BF)

将导致

如何在添加参与者时绘制顺序贝叶斯因子

How to plot a Sequential Bayes Factor as participants are added

plot

r

bayesian