本文整理汇总了Python中lifelines.estimation.CoxPHFitter._get_efron_values方法的典型用法代码示例。如果您正苦于以下问题:Python CoxPHFitter._get_efron_values方法的具体用法?Python CoxPHFitter._get_efron_values怎么用?Python CoxPHFitter._get_efron_values使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类lifelines.estimation.CoxPHFitter
的用法示例。
在下文中一共展示了CoxPHFitter._get_efron_values方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: test_efron_computed_by_hand_examples
# 需要导入模块: from lifelines.estimation import CoxPHFitter [as 别名]
# 或者: from lifelines.estimation.CoxPHFitter import _get_efron_values [as 别名]
def test_efron_computed_by_hand_examples(self, data_nus):
cox = CoxPHFitter()
X = data_nus['x'][:, None]
T = data_nus['t']
E = data_nus['E']
# Enforce numpy arrays
X = np.array(X)
T = np.array(T)
E = np.array(E)
# Want as bools
E = E.astype(bool)
# tests from http://courses.nus.edu.sg/course/stacar/internet/st3242/handouts/notes3.pdf
beta = np.array([[0]])
l, u = cox._get_efron_values(X, beta, T, E)
l = -l
assert np.abs(l[0][0] - 77.13) < 0.05
assert np.abs(u[0] - -2.51) < 0.05
beta = beta + u / l
assert np.abs(beta - -0.0326) < 0.05
l, u = cox._get_efron_values(X, beta, T, E)
l = -l
assert np.abs(l[0][0] - 72.83) < 0.05
assert np.abs(u[0] - -0.069) < 0.05
beta = beta + u / l
assert np.abs(beta - -0.0325) < 0.01
l, u = cox._get_efron_values(X, beta, T, E)
l = -l
assert np.abs(l[0][0] - 72.70) < 0.01
assert np.abs(u[0] - -0.000061) < 0.01
beta = beta + u / l
assert np.abs(beta - -0.0335) < 0.01