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)
