Add R code for Q1-4