怎么用python计算上述风险比

可以使用Python中的生存分析库 lifelines 来计算上述风险比。具体实现步骤如下：

1. 安装 lifelines 库：

pip install lifelines

2. 导入需要的库：

from lifelines import CoxPHFitter
import pandas as pd

3. 构建数据集：

data = {
    'Variablea': ['Waist circumference change']*3,
    'Waist circumference change': ['Loss (>3.0)', 'Stable (≤3.0)', 'Gain (>3.0)'],
    'DFTJ cohort No. of events/person years': [141, 184, 198],
    'Kailuan study No. of events/person years': [1102, 991, 1412],
    'Pooled results No. of events/person years': [1243, 1175, 1610]
}
df = pd.DataFrame(data)

4. 计算风险比：

cph = CoxPHFitter()
cph.fit(df, duration_col=None, event_col=None)
print(cph.summary)

运行后将得到各个类别相对于“稳定腰围”的风险比以及其置信区间。其中，“Loss (>3.0)”和“Gain (>3.0)”为需要进行比较的两组，分别对应“腰围减少>3.0 cm”和“腰围增加>3.0 cm”。