Enhancements
parallel_backend argument to adjustedsurv() and adjustedcif() to allow usage of FORK based parallel backends on Linux and MacOS, as suggested by @HUI950319Bug Fixes
adjustedsurv() with method="iptw_km", when one or more rows of the supplied data contained an event exactly at time 0.Enhancements
as_ggsurvplot_df(), even when method = "km" was used. Additionally, with multiple imputation, the correctly pooled values are now included, instead of nothing.ggplot2 based functions, so that the weird messages about unused labels no longer appear with newer versions of ggplot2.Bug Fixes
plot.adjustedsurv() with the times argument in case a dataset with multiple imputation was supplied.cuminc() function when using adjustedcif() with method="aalen_johansen". Essentially, the step functions were evaluated at the time before the actual step, so that any results obtained by using the previous versions were off by one time step. With many points in time the changes should be small, but with only a few events there might be substantial differences.adjustedcif() was used with non-competing risks data and method="iptw_pseudo", method="aiptw_pseudo" or method="direct_pseudo".New features
extend_to_last argument to method="iptw_km" in adjustedsurv(), which allows users to draw the IPTW survival curves up to the last observed point in time per group whether that time was censored or not. Unpublished simulation studies have shown that estimates beyond the last event time are very unstable, which is why in previous versions of this package (<= 0.11.2), this argument did not exist, but was essentially set to FALSE. To get the same results as with old versions, set this argument to FALSE.Refactor:
difference and ratio argument from adjusted_rmst() and adjusted_rmtl()Bug Fixes
censoring_model argument being ignored when using method="aiptw" in adjustedcif().risk_table=TRUE with risk_table_type="n_events" or risk_table_type="n_cens", events or censored observations at the same point in time were mistakenly excluded (a < was used internally instead of <=). This has been fixed, potentially leading to different results than earlier with discrete values of time.Bug Fixes
...) syntax not being correctly evaluated. For example, supplying estimand="ATT" to adjustedsurv() with method="iptw_km" and a formula in the treatment_model argument would result in the estimand argument not being passed to weightit().plot.adjustedsurv() when the adjustedsurv object was created in a function or loopEnhancements
New features
plot.adjustedsurv() (all arguments starting with risk_table)plot_rmst_curve() functionplot_rmtl_curve() functionadjusted_surv_quantile() functionRefactored
adjsurv and adjcif output objects of adjustedsurv() and adjsutedcif() respectively to adjdifference and ratio arguments into one contrast argument in adjusted_rmst(), adjusted_rmtl(), adjusted_surv_quantile(), plot_rmst_curve() and plot_rmtl_curve() functionstmle in adjustedsurv() and adjustedcif() due to concrete being removed from CRANDocumentation
Enhancements
iso_reg and force_bounds to adjustedsurv() and adjustedcif() functions to allow applying correction techniques outside plotting as wellvariable, ev_time or event contain missings, includes the new mi_extrapolation argument in adjustedsurv() and adjustedcif()ratio argument to adjusted_rmst() and adjusted_rmtl() functionsto values in adjusted_rmst() and adjusted_rmtl(), which speeds up plot_rmst_curve() and plot_rmtl_curve() considerablyBug Fixes
adjusted_rmst(), adjusted_rmtl() and output of adjusted_curve_test() and functions that rely on those functions. This is fixed now, but may lead to slightly different estimates than in previous versions.New Features
surv_tmle and cif_tmle, based on the concrete packagesurv_iv_2SRIFsurv_prox_iptw, surv_prox_aiptw based on code from Andrew Yingadjusted_curve_ratio() functionplot_curve_ratio() functionDocumentation
WeightIt as suggested by Noah GreiferWeightIt packageprint method to be equal to summary method