R code and output goes here