Eu usei o rjags para executar o MCMC em um modelo, especificado na linguagem JAGS. Existe uma boa maneira de extrair esse modelo e executar previsões com ele (usando as distribuições posteriores dos meus parâmetros)? Posso re-especificar o modelo em R e conectar os modos dos meus parâmetros posteriores; Só estou me perguntando se existe uma maneira menos redundante de fazer isso.
Eu acredito que http://sourceforge.net/p/mcmc-jags/discussion/610037/thread/0ecab41c está fazendo a mesma pergunta.
Respostas:
Geralmente você pode fazer as previsões no JAGS. Abaixo está um exemplo de regressão com VEF (algo relacionado à capacidade pulmonar) como variável dependente e indicador de idade e tabagismo como preditores.
VEF20s e VEF20ns são os valores previstos para um fumante de 20 anos e um não fumante de 20 anos.
Exemplo de: Ideias Bayesianas e Análise de Dados
fonte