data analysis using r code