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

Question

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

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

Answer 1

你可以执行

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

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