Simulazione RJAGS per la regressione di Poisson

Nel video precedente abbiamo costruito un modello di regressione di Poisson del volume \(Y\)_i in funzione dello stato di giorno feriale \(X\)_i e della temperatura \(Z\)_i:

• verosimiglianza: \(Y\)_i \(\sim Pois(l\)_i) dove \(log(l\)_i\() = a + b X\)_i \(+ c Z\)_i
• prior: \(a \sim N(0, 200^2)\), \(b \sim N(0, 2^2)\) e \(c \sim N(0, 2^2)\)

Combinando gli spunti ricavati dai dati osservati RailTrail e i prior specificati qui, definirai, compilerai e simulerai un modello posterior di questa relazione usando RJAGS. Per metterti alla prova in quest’ultima simulazione RJAGS del corso, ti forniremo meno codice del solito!

I dati RailTrail sono nel tuo workspace.