我有一个 Stan 模型(称为fit
),我得到了一些奇怪的结果。我知道我可能遗漏了一些东西,但我看不出在哪里和什么。在这种情况下,我想查看参数sigma_h
。
如果我用来as.matrix
从中获取样本,fit
我会得到以下结果:
posterior <- as.matrix(fit)
trace_sigma_h <- posterior[,colnames(posterior) %in% "sigma_h"]
ts.plot(trace_sigma_h)
这与我使用该rstan::traceplot
函数时得到的图相同。
如果我从对象中提取样本fit
并绘制它,我会得到
samples <- extract(fit,inc_warmup = FALSE)
ts.plot(samples$sigma_h)
这些地块并不相同。我错过了什么?fit
仅将对象设为矩阵和提取样本之间有区别吗?