R code for paper

R code for analysis and graphs