Updated all R codes and graphs