应用于stan对象的extract()和as.matrix()之间的差异

时间:2017-10-27 08:16:05

标签: r rstan

我有一个Stan模型(称为fit),我得到了一些奇怪的结果。我知道我可能遗漏了一些东西,但我无法看到哪里和那些东西。在这种情况下,我想查看参数sigma_h

如果我使用as.matrixfit获取样本,我会得到以下结果:

posterior <- as.matrix(fit) 
trace_sigma_h <- posterior[,colnames(posterior) %in% "sigma_h"]
ts.plot(trace_sigma_h)

enter image description here

这与我使用rstan::traceplot函数时获得的情节相同。

如果我从fit对象中提取样本并绘制它,我会得到

samples <- extract(fit,inc_warmup = FALSE)
ts.plot(samples$sigma_h) 

enter image description here

这些情节不一样。我错过了什么?仅仅将fit对象设为矩阵并提取样本之间是否存在差异?

0 个答案:

没有答案