我在STAN中使用贝叶斯模型来估计测量误差回归模型。诊断说这条链融合了。但是,当我查看参数与对数后验图时。我发现链没有探索对数后验最大的空间。该图显示了参数与后验图。我们可以看到对数后部仍在向右增加。
我试图了解是什么导致了此行为,因为链条正在收敛到这个意外点,这不是此边际分布的后验模式。为什么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
}