my R codes and figures