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

Question

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

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

Answer 1

你可以执行

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

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