python生存率列线图

以下是使用Python绘制生存率列线图的示例代码：

import matplotlib.pyplot as plt

def survival_rate(ls_hr, ls_xvar, basic_rate):
    ls_beta = [np.log(x) for x in ls_hr]
    ls_weight = [a * b for a, b in zip(ls_beta, ls_xvar)]
    pi = sum(ls_weight)
    survival_rate = basic_rate ** np.exp(pi)
    return survival_rate

# 示例数据
ls_hr = [0.5, 0.6, 0.7, 0.8, 0.9]
ls_xvar = [1, 2, 3, 4, 5]
basic_rate = 0.95

# 计算生存率
survival_rates = [survival_rate(ls_hr, ls_xvar, basic_rate) for _ in range(10)]

# 绘制列线图
plt.bar(range(len(survival_rates)), survival_rates)
plt.xlabel('Sample')
plt.ylabel('Survival Rate')
plt.title('Survival Rate Distribution')
plt.show()

这段代码首先定义了一个计算生存率的函数`survival_rate`，然后使用示例数据计算了10个样本的生存率。最后，使用Matplotlib库绘制了一个列线图，横轴表示样本，纵轴表示生存率。