我有一个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
对象设为矩阵并提取样本之间是否存在差异?