4

我有一个 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仅将对象设为矩阵和提取样本之间有区别吗?

4

0 回答 0