我正在使用stancmd,这是进行贝叶斯推理的好工具。就我而言,我确实需要推断参数的均值和方差信息。我了解这对于MCMC是可能的,只需从out.csv中读取stddev。例如,
# Output
# mean se_mean sd 2.5% 25% 50% 75% 97.5%
n_eff Rhat
# alpha 3.33 0.01 0.40 2.57 3.06 3.33 3.59 4.13
1229 1
# beta 2.14 0.01 0.40 1.37 1.89 2.14 2.40 2.98
1470 1
# sigma 1.92 0.01 0.27 1.45 1.71 1.90 2.09 2.51
1211 1
# lp__ -31.92 0.05 1.30 -35.27 -32.50 -31.63 -30.96 -30.43
769 1
如何使用变分推理来做同样的事情?我只能从标准输出文件中获取均值。
非常感谢