Simulação com RJAGS para regressão multivariada

Considere o seguinte modelo Bayesiano do volume \(Y\)_i por status de dia da semana \(X\)_i e temperatura \(Z\)_i:

• verossimilhança: \(Y\)_i \(\sim N(m\)_i, \(s^2)\) em que \(m\)_i \(= a + b X\)_i \(+ c Z\)_i .
• priors: \(a \sim N(0, 200^2)\), \(b \sim N(0, 200^2)\), \(c \sim N(0, 20^2)\), \(s \sim Unif(0, 200)\)

Sua exploração anterior da relação entre volume, weekday e hightemp nos dados RailTrail trouxe algumas pistas sobre esse relacionamento. Você vai combinar isso com o que os priors indicam para desenvolver um modelo a posteriori dessa relação usando RJAGS. Os dados RailTrail estão no seu ambiente de trabalho.