MCMC链收敛,但是对数后验拓扑不适合收敛

时间:2019-04-10 16:10:13

标签: mcmc stan

我在STAN中使用贝叶斯模型来估计测量误差回归模型。诊断说这条链融合了。但是,当我查看参数与对数后验图时。我发现链没有探索对数后验最大的空间。该图显示了参数与后验图。我们可以看到对数后部仍在向右增加。
Imgur

我试图了解是什么导致了此行为,因为链条正在收敛到这个意外点,这不是此边际分布的后验模式。为什么HMC / NUTS不探索该参数对数后验较高的右侧空间?

我适合的型号是:

data {
  int<lower=0> N;   // Sample size
  vector[N] X;      // Predictor
  vector[N] Y;      // Response
  real<lower=0> sigmax; // Measurement error in X
}
parameters {
  real beta0;       // intercept
  real beta1;       // X coefficient
  real<lower=0> sigma;  // error in Y
  vector[N] Xu;         // latent unmeasured X
}
model {
  X ~ normal(Xu,sigmax);
  Y ~ normal(beta0+beta1*Xu,sigma);  // likelihood
}

0 个答案:

没有答案