def test_p_value_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') # p-values actual_p = cf._compute_p_values() expected_p = np.array([0.1847, 0.7644, 0.0730, 0.00]) npt.assert_array_almost_equal(actual_p, expected_p, decimal=2)