def print_comparative_outcomes(sim_outcomes_none, sim_outcomes_anti): """ prints average increase in survival time, discounted cost, and discounted utility under combination therapy compared to mono therapy :param sim_outcomes_mono: outcomes of a cohort simulated under mono therapy :param sim_outcomes_combo: outcomes of a cohort simulated under combination therapy """ # increase in mean survival time under combination therapy with respect to mono therapy increase_survival_time = Stat.DifferenceStatIndp( name='Increase in mean survival time', x=sim_outcomes_anti.survivalTimes, y_ref=sim_outcomes_none.survivalTimes) # estimate and CI estimate_CI = increase_survival_time.get_formatted_mean_and_interval( interval_type='c', alpha=D.ALPHA, deci=2) print( "Increase in mean survival time and {:.{prec}%} confidence interval:". format(1 - D.ALPHA, prec=0), estimate_CI) # increase in mean discounted cost under combination therapy with respect to mono therapy increase_discounted_cost = Stat.DifferenceStatIndp( name='Increase in mean discounted cost', x=sim_outcomes_anti.costs, y_ref=sim_outcomes_none.costs) # estimate and CI estimate_CI = increase_discounted_cost.get_formatted_mean_and_interval( interval_type='c', alpha=D.ALPHA, deci=2, form=',') print( "Increase in mean discounted cost and {:.{prec}%} confidence interval:" .format(1 - D.ALPHA, prec=0), estimate_CI) # increase in mean discounted utility under combination therapy with respect to mono therapy increase_discounted_utility = Stat.DifferenceStatIndp( name='Increase in mean discounted utility', x=sim_outcomes_anti.utilities, y_ref=sim_outcomes_none.utilities) # estimate and CI estimate_CI = increase_discounted_utility.get_formatted_mean_and_interval( interval_type='c', alpha=D.ALPHA, deci=2) print( "Increase in mean discounted utility and {:.{prec}%} confidence interval:" .format(1 - D.ALPHA, prec=0), estimate_CI) # increase in mean discounted utility under combination therapy with respect to mono therapy increase_num_strokes = Stat.DifferenceStatIndp( name='Increase in mean discounted utility', x=sim_outcomes_anti.nTotalStrokes, y_ref=sim_outcomes_none.nTotalStrokes) # estimate and CI estimate_CI = increase_num_strokes.get_formatted_mean_and_interval( interval_type='c', alpha=D.ALPHA, deci=2) print( "Increase in number of strokes and {:.{prec}%} confidence interval:". format(1 - D.ALPHA, prec=0), estimate_CI)
def print_comparative_outcomes(multi_cohort_no_drug, multi_cohort_with_drug): """ prints expected and percentage increase in average survival time when drug is available :param multi_cohort_no_drug: multiple cohorts simulated when drug is not available :param multi_cohort_with_drug: multiple cohorts simulated when drug is available """ # increase in survival time increase_stat = Stat.DifferenceStatIndp( name='Increase in mean survival time', x=multi_cohort_with_drug.multiCohortOutcomes.meanSurvivalTimes, y_ref=multi_cohort_no_drug.multiCohortOutcomes.meanSurvivalTimes) # mean and prediction interval mean = increase_stat.get_mean() pred_int = increase_stat.get_PI(alpha=D.ALPHA) print( "Expected increase in mean survival time (years) and {:.{prec}%} prediction interval:" .format(1 - D.ALPHA, prec=0), mean, pred_int) # % increase in mean survival time relative_diff_stat = Stat.RelativeDifferenceIndp( name='% increase in mean survival time', x=multi_cohort_with_drug.multiCohortOutcomes.meanSurvivalTimes, y_ref=multi_cohort_no_drug.multiCohortOutcomes.meanSurvivalTimes) # estimate and prediction interval mean = relative_diff_stat.get_mean() pred_int = relative_diff_stat.get_PI(alpha=D.ALPHA) print( "Expected percentage increase in mean survival time and {:.{prec}%} confidence interval:" .format(1 - D.ALPHA, prec=0), mean, pred_int)
def print_comparative_outcomes(cohort_no_drug, cohort_with_drug): """ prints expected and percentage increase in survival time when drug is available :param cohort_no_drug: a cohort simulated when drug is not available :param cohort_with_drug: a cohort simulated when drug is available """ # create a difference statistics for the increase in survival time increase_stat = Stat.DifferenceStatIndp( name='Increase in survival time', x=cohort_with_drug.cohortOutcomes.survivalTimes, y_ref=cohort_no_drug.cohortOutcomes.survivalTimes) # get mean and t-based confidence interval for the increase in survival time mean = increase_stat.get_mean() conf_int = increase_stat.get_t_CI(alpha=D.ALPHA) print( "Average increase in survival time (years) and {:.{prec}%} confidence interval:" .format(1 - D.ALPHA, prec=0), mean, conf_int)
def mytest_diff_stat_indp(x, y): # define stat = Stat.DifferenceStatIndp(x, y, name='Test DifferenceStatIndp') print('Testing DifferenceStatIndp:') print_results(stat)