如何从 rstanarm 对象中提取 stan 代码

Question

是否有可能在 rstanarm 中提取用于 MCMC 采样的 stan 代码？

我想将我自己的模型参数化和先前的选择与 rstanarm 中使用的参数化进行比较。

Answer 1

你可以执行

library(rstanarm)
example(example_model)
rstan::get_stanmodel(example_model$stanfit)

查看 Stan 代码或在 GitHub 上查看。但是，它并不打算供人类阅读，并且由于 #include 语句和支持的大量选项，可读性不是特别好。