def test_se_against_Survival_Analysis_by_John_Klein_and_Melvin_Moeschberger(self):
# see table 8.1 in Survival Analysis by John P. Klein and Melvin L. Moeschberger, Second Edition
df = load_larynx()
cf = CoxPHFitter()
cf.fit(df, duration_col='time', event_col='death')
# standard errors
actual_se = cf._compute_standard_errors().values
expected_se = np.array([[0.0143, 0.4623, 0.3561, 0.4222]])
npt.assert_array_almost_equal(actual_se, expected_se, decimal=2)
def test_se_against_Survival_Analysis_by_John_Klein_and_Melvin_Moeschberger(self):
# see table 8.1 in Survival Analysis by John P. Klein and Melvin L. Moeschberger, Second Edition
df = load_larynx()
cf = CoxPHFitter(normalize=False)
cf.fit(df, duration_col='time', event_col='death')
# standard errors
actual_se = cf._compute_standard_errors().values
expected_se = np.array([[0.0143, 0.4623, 0.3561, 0.4222]])
npt.assert_array_almost_equal(actual_se, expected_se, decimal=2)
