我正在运行线性混合模型分析，其中我有一个因变量（瞳孔均值），两个固定因子（条件和相类型）和一个随机因子（subj）。我想用R中的BayesFactor包的lmBF函数计算相关的贝叶斯因子。但是，变量full_BF_pupil每次运行时都会给出不同的结果。
例如：在第一次运行时：1.386933e + 137±2.19％, 2nd run: 1.381459e+137 ±2.18%
即使差异很小，由于对该值进行了其他计算，最终结果也更加不同。这是什么原因我是否有更可靠的功能？ 谢谢，这是我们感兴趣的代码行。
full_BF_pupil = lmBF(pupil_mean_bin~ cond*phasetype+ subj, data = BR_pre, whichRandom="subj")