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

Question

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

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

Answer 1

你可以执行

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

查看Stan代码或在GitHub上查看它。但是，由于#include语句和所支持的大量选项，它并不打算被人类读取，也不是特别易读。