R code for creating maps using the census of religion dataset