add a function similar to [summary_psm](https://github.com/Roche/flexsurvPlus/blob/main/R/summaryPSM.R) which will return a tibble of: - KM estimates of survival - Hazard estimates (pwe) - Hazard estimates (smoothed) - Cumulative hazard That could then be used for diagnostic plots