本例对样本数据 VA Lung Cancer.jmp 使用比例风险模型。该数据是从随机化的临床试验收集的,在这些试验中对无法手术的男性肺癌患者实施了标准或新的(检验)化学疗法(治疗)。该试验的主要目的是评估治疗类型是否对生存时间有影响,特别关注肿瘤类型(细胞类型)的影响。有关该数据集的其他详细信息,请参见 Prentice (1973) 和 Kalbfleisch and Prentice (2002)。
年龄诊断时间KPS 是连续测量值,细胞类型治疗之前治疗是分类(名义型)变量。细胞类型的四个名义型水平包括“Adeno”、“Large”、“Small”和“Squamous”。
1.
选择帮助 > 样本数据库,然后打开 VA Lung Cancer.jmp
2.
选择分析 > 可靠性和生存 > 拟合比例风险
3.
选择时间作为事件时间
4.
选择删失作为删失
5.
选择细胞类型治疗之前治疗年龄诊断时间KPS,然后点击添加
6.
点击运行
图 15.6 具有多个效应和水平的比例风险模型的报表窗口
在“效应似然比检验”报表中,“概率>卡方”值指示 KPS细胞类型的至少一个水平是显著的,而治疗之前治疗年龄诊断时间等效应不显著。
图 15.7 “风险比”报表
图 15.7显示了连续效应(年龄诊断时间KPS)以及名义型效应(细胞类型治疗之前治疗)的风险比。为了便于说明,我们侧重 VA Lung Cancer.jmp 样本数据的连续效应年龄和具有四个水平的名义型效应(细胞类型)。
对于 VA Lung Cancer.jmp 样本数据中的连续效应年龄,按以下方式计算风险比:
exp(β) = exp(-0.0085494) = 0.991487
exp[β(xmax - xmin)] = exp(-0.0085494 * 47) = 0.669099
对于名义型效应细胞类型,计算所有水平对并将结果显示在“‘细胞类型’的风险比”表中。请注意,对于包含 k 个水平的分类变量,仅使用 k -1 个设计变量或水平。在“参数估计值”表中,只显示了细胞类型四个水平中的三个水平(“Adeno”、“Large”和“Small”)的参数估计值。未显示“Squamous”水平,但是它是其他估计值之和的负数。此处是“细胞类型的风险比”计算的两个示例:
Large/Adeno = exp(βLarge)/exp(βAdeno) = exp(-0.2114757)/exp(0.57719588) = 0.4544481
Squamous/Adeno = exp[-(βAdeno + βLarge + βSmall)]/exp(βAdeno)