Simulación con RJAGS para regresión multivariante
Considera el siguiente modelo bayesiano del volumen \(Y\)i según el estado de día laborable \(X\)i y la temperatura \(Z\)i:
- verosimilitud: \(Y\)i \(\sim N(m\)i, \(s^2)\) donde \(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)\)
Tu exploración previa de la relación entre volume, weekday y hightemp en los datos de RailTrail te dio algunas pistas sobre esta relación. Ahora combinarás esto con la información de las priors para construir un modelo posterior de esta relación usando RJAGS. Los datos de RailTrail están en tu espacio de trabajo.
Este ejercicio forma parte del curso
Modelado bayesiano con RJAGS
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# DEFINE the model
rail_model_2 <-