R scripts for analyzing TB case control genetic data