R code for data analysis as reported in manuscript.