Как извлечь стандартный код из объекта rstanarm

Question

Есть ли возможность извлечь стандартный код, используемый для выборки MCMC в rstanarm?

Я хотел бы сравнить свою собственную параметризацию модели и предыдущих вариантов с той, что использовалась в rstanarm.

Answer 1

Вы можете выполнить

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

чтобы увидеть код Стэна или посмотреть его на GitHub. Тем не менее, он не предназначен для чтения людьми и не особенно читаем из-за #include заявления и большое количество опций поддерживается.