--- title: "HW 10" author: "Name here" output: html_document --- Please use D2L to turn in both the HTML output and your R Markdown file in. We will revisit the Swiss birds dataset for this lab to construct a Poisson regression model for abundance of the Willow Tit. ```{r} library(knitr) swiss.birds <- read.csv('http://math.montana.edu/ahoegh/teaching/stat491/data/willowtit2013_count.csv') kable(head(swiss.birds)) ``` This dataset contains 242 sites and 6 variables. The only difference from the lab is: - bird.count, count variable for abndance of birds observed, rather than a binary outcome #### 1. (5 points) Model Specification Clearly write out the model, using proper notation for the variables in bird dataset. You don't need use all variables, but you should state which are included. \vfill #### 2. (5 points) Priors Describe and justify the necessary priors for this model. \vfill #### 3. (5 points) Fit MCMC Fit the JAGS code for this model. You will have to put this together following the specification in the previous examples, but the following statement can be used for the sampling model portion. ```{r, eval = F} model { for (i in 1:Ntotal) { y[i] ~ dpois(mu[i]) mu[i] <- exp(beta0 + sum( beta[1:Nx] * x[i,1:Nx] )) } # priors inserted here } ``` #### 4. (5 points) Summarize inferences from model Talk about the model and discuss which and how predictor variables influence the observation of a bird. Furthermore, discuss the differences between this model and the logistic regression setting.