subclinical TB project R script