def print_comparative_outcomes(multi_trees_P, multi_trees_N): """ prints average increase in discounted cost, and discounted utility under combination therapy compared to mono therapy :param multi_trees_P: outcomes of a multi-trees simulated under Palivizumab :param multi_trees_N: outcomes of a multi-trees simulated under no prophylaxis """ # increase in mean discounted cost under combination therapy with respect to mono therapy increase_mean_discounted_cost = Stat.DifferenceStatPaired( name='Increase in mean discounted cost', x=multi_trees_P.PCosts, y_ref=multi_trees_N.NCosts) # estimate and PI estimate_PI = increase_mean_discounted_cost.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2, form=',') print("Increase in mean discounted cost and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean discounted QALY under combination therapy with respect to mono therapy increase_mean_discounted_qaly = Stat.DifferenceStatPaired( name='Increase in mean discounted QALY', x=multi_trees_P.PQALYs, y_ref=multi_trees_N.NQALYs) # estimate and PI estimate_PI = increase_mean_discounted_qaly.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean discounted utility and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI)
def print_comparative_outcomes(multi_cohort_outcomes_no, multi_cohort_outcomes_asp): """ prints average increase in survival time, discounted cost, and discounted utility under combination therapy compared to mono therapy :param multi_cohort_outcomes_mono: outcomes of a multi-cohort simulated under mono therapy :param multi_cohort_outcomes_combo: outcomes of a multi-cohort simulated under combination therapy """ # increase in mean survival time under combination therapy with respect to mono therapy increase_mean_survival_time = Stat.DifferenceStatPaired( name='Increase in mean survival time', x=multi_cohort_outcomes_asp.meanSurvivalTimes, y_ref=multi_cohort_outcomes_no.meanSurvivalTimes) # estimate and PI estimate_PI = increase_mean_survival_time.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean survival time and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean discounted cost under combination therapy with respect to mono therapy increase_mean_discounted_cost = Stat.DifferenceStatPaired( name='Increase in mean discounted cost', x=multi_cohort_outcomes_asp.meanCosts, y_ref=multi_cohort_outcomes_no.meanCosts) # estimate and PI estimate_PI = increase_mean_discounted_cost.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2, form=',') print("Increase in mean discounted cost and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean discounted QALY under combination therapy with respect to mono therapy increase_mean_discounted_qaly = Stat.DifferenceStatPaired( name='Increase in mean discounted QALY', x=multi_cohort_outcomes_asp.meanQALYs, y_ref=multi_cohort_outcomes_no.meanQALYs) # estimate and PI estimate_PI = increase_mean_discounted_qaly.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean discounted utility and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI)
def mytest_diff_stat_paired(x, y): # define stat = Stat.DifferenceStatPaired('Test DifferenceStatPaired', x, y) print('Testing DifferenceStatPaired:') print_results(stat)
def print_comparative_outcomes(multi_cohort_outcomes_SOC, multi_cohort_outcomes_NSB): """ prints average increase in survival time, discounted cost, and discounted utility under combination therapy compared to mono therapy :param multi_cohort_outcomes_mono: outcomes of a multi-cohort simulated under mono therapy :param multi_cohort_outcomes_combo: outcomes of a multi-cohort simulated under combination therapy """ # increase in mean survival time under combination therapy with respect to mono therapy increase_mean_survival_time = Stat.DifferenceStatPaired( name='Increase in mean survival time', x=multi_cohort_outcomes_NSB.meanSurvivalTimes, y_ref=multi_cohort_outcomes_SOC.meanSurvivalTimes) # estimate and PI estimate_PI = increase_mean_survival_time.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean survival time and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean number of hospitalizations increase_mean_nHosp = Stat.DifferenceStatPaired( name='Increase in mean number of kids hospitalized', x=multi_cohort_outcomes_NSB.nHospitalized, y_ref=multi_cohort_outcomes_SOC.nHospitalized) # estimate and PI estimate_PI = increase_mean_nHosp.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean number of hospitalizations and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean number of YLLs increase_mean_YLLs = Stat.DifferenceStatPaired( name='Increase in mean number of Years of Life Lost', x=multi_cohort_outcomes_NSB.totalYLL, y_ref=multi_cohort_outcomes_SOC.totalYLL) # estimate and PI estimate_PI = increase_mean_YLLs.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean number of Years of Life Lost and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean number of patients w/ TBM increase_mean_TBM = Stat.DifferenceStatPaired( name='Increase in mean number of cases of TB meningitis', x=multi_cohort_outcomes_NSB.nTBM, y_ref=multi_cohort_outcomes_SOC.nTBM) # estimate and PI estimate_PI = increase_mean_TBM.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean number of cases of TB meningitis and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean 2-month mortality increase_mean_mortality2Month = Stat.DifferenceStatPaired( name='Increase in mean 2-month mortality', x=multi_cohort_outcomes_NSB.mortality2Month, y_ref=multi_cohort_outcomes_SOC.mortality2Month) # estimate and PI estimate_PI = increase_mean_mortality2Month.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean 2-month mortality and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean 2-year mortality increase_mean_mortality2Year = Stat.DifferenceStatPaired( name='Increase in mean 2-year mortality', x=multi_cohort_outcomes_NSB.mortality2Years, y_ref=multi_cohort_outcomes_SOC.mortality2Years) # estimate and PI estimate_PI = increase_mean_mortality2Year.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean 2-year mortality and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean 5-year mortality increase_mean_mortality5Year = Stat.DifferenceStatPaired( name='Increase in mean 5-year mortality', x=multi_cohort_outcomes_NSB.mortality5Years, y_ref=multi_cohort_outcomes_SOC.mortality5Years) # estimate and PI estimate_PI = increase_mean_mortality5Year.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2) print("Increase in mean 5-year mortality and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI) # increase in mean discounted cost under combination therapy with respect to mono therapy increase_mean_discounted_cost = Stat.DifferenceStatPaired( name='Increase in mean discounted cost', x=multi_cohort_outcomes_NSB.meanCosts, y_ref=multi_cohort_outcomes_SOC.meanCosts) # estimate and PI estimate_PI = increase_mean_discounted_cost.get_formatted_mean_and_interval(interval_type='p', alpha=D.ALPHA, deci=2, form=',') print("Increase in mean discounted cost and {:.{prec}%} uncertainty interval:" .format(1 - D.ALPHA, prec=0), estimate_PI)