Add exponentiate mortality model for mortality model (in time_frame_main.R) and add event per patient-year (incidence rate) for both group in time_to_event.R.