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 mytest_relativeDiff_stat_indp(x, y):
# define
stat = Stat.RelativeDifferenceIndp(x, y, name='Test RelativeDifferenceIndp')
print('Testing RelativeDifferenceIndp:')
print_results(stat)
