This work presents the methodologies for integrating environmental variables into Bayesian GAM prevalence disease modeling. I includes scripts for initial data exploration, variable selection, model fitting and plots. Data files used for the analysis are also included.