R code for data analysis of Beijing intestinal disease dataset