Gibbs Sampler 输出:多少个马尔可夫链?

机器算法验证 贝叶斯 条件概率 马尔可夫链蒙特卡罗 吉布斯 马尔可夫链
2022-04-05 09:48:03

当运行带有两个完整条件的 Gibbs 采样器(对于迭代)时,我得到输出 .n=200x=(x1(n),x2(n))n=1,...,200

所以是 Gibbs Markov 链的实现,即所谓的 Gibbs 序列。但是是链的两个实现?x(x1(n))n[1,...,200],(x2(n))n[1,...,100]

1个回答

当子链本身仍然是马尔可夫链时,这个两块吉布斯采样器是唯一的通用情况,因为是通过内核生成的 有关更多详细信息, 请参阅我们的 MCMC 书,但这是一种交错属性的情况,它也保证 ' s 与随时间差而减小的相关性呈正相关,并且 Rao-Blackwellisation 总是会降低结果估计的方差。(X1(n))

K(x1,x1)=f2(x2|x1)f1(x1|x2)dx2
X1(n)

两条马尔可夫链被称为 具有交错性质(或 交错性质)彼此共轭,如果(X(t))(Y(t))

  1. X(t)有条件地独立于 X(t+1)Y(t)
  2. Y(t1)有条件地独立于 Y(t)X(t)
  3. (X(t),Y(t1))在平稳性下同分布。(X(t),Y(t))