def print_comparative_outcomes(sim_outcomes_none, sim_outcomes_treat):
""" 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 number of strokes
change_n_cured = Stat.DifferenceStatIndp(
name='Change in number of patients cured',
x=sim_outcomes_treat.nTotalCured,
y_ref=sim_outcomes_none.nTotalCured)
estimate_CI = change_n_cured.get_formatted_mean_and_interval(
interval_type='c', alpha=D.ALPHA, deci=2)
print(
"Change in the number of patients cured:".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_treat.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)
def print_comparative_outcomes(sim_outcomes_amino, sim_outcomes_immuno):
""" prints average increase in survival time, discounted cost, and discounted utility
under immunosuppresive therapy compared to aminosalicylate therapy
:param sim_outcomes_amino: outcomes of a cohort simulated under mono therapy
:param sim_outcomes_immuno: outcomes of a cohort simulated under combination therapy
"""
# increase in mean survival time under immuno therapy with respect to amino therapy
increase_survival_time = Stat.DifferenceStatIndp(
name='Increase in mean survival time',
x=sim_outcomes_immuno.survivalTimes,
y_ref=sim_outcomes_amino.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 immuno therapy with respect to amino therapy
increase_discounted_cost = Stat.DifferenceStatIndp(
name='Increase in mean discounted cost',
x=sim_outcomes_immuno.costs,
y_ref=sim_outcomes_amino.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 number of patients alive under immuno therapy with respect to amino therapy
increase_patients_alive = Stat.DifferenceStatIndp(
name='Increase in mean number of patients alive',
x=sim_outcomes_immuno.numPatientsAlive,
y_ref=sim_outcomes_amino.numPatientsAlive)
# estimate and CI
estimate_CI = increase_patients_alive.get_formatted_mean_and_interval(
interval_type='c', alpha=D.ALPHA, deci=2)
print(
"Increase in mean number of patients alive and {:.{prec}%} confidence interval:"
.format(1 - D.ALPHA, prec=0), estimate_CI)
def print_comparative_outcomes(sim_outcomes_SOC, sim_outcomes_NSB):
"""
:param sim_outcomes_SOC: outcomes of a cohort simulated under SOC diagnostics
:param sim_outcomes_NSB: outcomes of a cohort simulated under NSB diagnostics
"""
# 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_combo.survivalTimes,
y_ref=sim_outcomes_mono.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_combo.costs,
y_ref=sim_outcomes_mono.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_combo.utilities,
y_ref=sim_outcomes_mono.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)
def print_comparative_outcomes(sim_outcomes_NO, sim_outcomes_ASP):
# increase in mean survival time under therapy with respect to no therapy
decrease_survival_time = Stat.DifferenceStatIndp(
name='Increase in mean survival time',
x=sim_outcomes_ASP.survivalTimes,
y_ref=sim_outcomes_NO.survivalTimes)
# estimate and CI
estimate_CI = decrease_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 therapy with respect to no therapy
increase_discounted_cost = Stat.DifferenceStatIndp(
name='Increase in mean discounted cost',
x=sim_outcomes_ASP.costs,
y_ref=sim_outcomes_NO.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 no therapy
increase_discounted_utility = Stat.DifferenceStatIndp(
name='Increase in mean discounted utility',
x=sim_outcomes_ASP.utilities,
y_ref=sim_outcomes_NO.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)
def mytest_diff_stat_indp(x, y):
# define
stat = Stat.DifferenceStatIndp('Test DifferenceStatIndp', x, y)
print('Testing DifferenceStatIndp:')
print_results(stat)