Add R-codes and plots