能告诉你什么
Crosstabs 在你的设备上运行与 SPSS COXREG 对应的过程,采用 Breslow 法处理结值,支持指示编码的分类协变量和可选分层,并给出风险比 Exp(B) 及其 95% 置信区间和协变量均值处的生存函数。
- 每个协变量的风险比 Exp(B):在其他协变量不变时,该协变量每增加一个单位,事件风险乘以多少。
- 通过得分卡方与似然比卡方检验,判断全部协变量合在一起是否优于不含协变量的模型。
- 每个系数的 Wald 检验;对于多水平分类协变量,还给出整个项的联合检验。
- 各协变量均处于平均值的个案的估计生存曲线和累积风险。
适用场景
适合使用
- 含删失的时间—事件结局,需要在比较组别时调整年龄、严重程度等协变量。
- 估计剂量、月消费等连续预测变量对事件风险的影响。
- 基线风险因研究中心或队列而异:将该变量设为分层,使每层拥有各自的基线风险。
改用其他方法
- 风险随时间明显交叉(某组早期较差、后期较好):比例风险假设不成立,可比较 Kaplan-Meier 曲线或将随访分段分析。
- 需要普通协变量之间的交互项:请先计算乘积项,再作为协变量纳入。
- 没有时间维度,只关心事件是否发生:使用二元 logistic 回归。
- 含抽样权重、分层和整群的调查数据:使用复杂抽样 Cox 回归。
在 Crosstabs 中如何操作
- 01
打开数据
打开每行一个个案的文件,并切换到“统计”模式。
- 02
选择 Cox 回归
在“生存分析”组中选择“Cox 回归”(COXREG)。
- 03
设置时间、状态和事件值
选择“时间”“状态”变量以及“表示事件已发生的值”。
- 04
选择协变量
勾选最多 12 个“协变量”。在“分类协变量(指示编码,最后一个类别为参照)”中勾选需要虚拟编码的变量;文本变量和名义变量默认已勾选。
- 05
按需添加分层
选择“分层(可选)”变量,使每个水平拥有独立的基线风险,同时共享回归系数。
- 06
选择图形并运行
将“图形”设为协变量均值处的生存或风险(累积);如有权重,将其声明为频数权重,然后点击“运行”。
计算示例
以下结果由 Crosstabs 引擎在构建本页时根据模拟数据实际计算得出,并非手工填写。
数据
120 名模拟患者交替分配到药物组或安慰剂组,年龄来自正态分布(均值 60,标准差 10)。生存时间服从指数模型:年龄每增加一岁风险上升 3%,药物组风险低 40%;随访在 36 个月截止。固定随机种子。
| 项目 | N | 百分比 |
|---|---|---|
| 事件 | 78 | 65.0% |
| 已删失 | 42 | 35.0% |
| 总计 | 120 | 100.0% |
| -2 对数似然 | 总体(得分)卡方 | df | 显著性 | 相对上一步的变化卡方 | df | 显著性 |
|---|---|---|---|---|---|---|
| 651.272 | 29.068 | 2 | <.001 | 29.027 | 2 | <.001 |
起始块 -2 对数似然:680.300。结值使用 Breslow 方法处理。
| 项目 | B | SE | Wald | df | 显著性 | Exp(B) | Exp(B) 的 95% 置信区间 |
|---|---|---|---|---|---|---|---|
| age | 0.028 | 0.013 | 4.466 | 1 | .035 | 1.028 | 1.002 – 1.055 |
| treatment(1) | −1.231 | 0.248 | 24.620 | 1 | <.001 | 0.292 | 0.180 – 0.475 |
treatment:指示编码;(1) = Drug;参照 = Placebo。
结果解读
模型总体显著(似然比 χ²(2) = 29.03,p < .001)。在年龄不变的条件下,药物使风险乘以 Exp(B) = 0.292(95% 置信区间 0.180 – 0.475);年龄每增加一岁,风险乘以 1.028(95% 置信区间 1.002 – 1.055)。风险比小于 1 表示事件的瞬时风险较低,而不是生存时间相差某个固定值。
输出表格说明
- 个案处理摘要
- 事件数、删失数及总数和百分比。经验上,每个系数至少需要约十个事件。
- 步骤摘要
- 使用逐步方法时:先列出起始模型,再列出每一步输入(得分统计量)或剔除(条件、LR 或 Wald 统计量)的协变量,以及模型的 -2 对数似然及其相对上一步的变化。
- 模型系数的综合检验
- 拟合模型的 -2 对数似然、总体得分卡方以及相对空模型的似然比变化,均附自由度和显著性。起始块的 -2 对数似然见脚注。
- 方程中的变量
- 每个系数的 B、标准误、Wald、自由度、显著性、Exp(B) 及其 95% 置信区间;多水平分类项另有一行联合检验。脚注列出指示编码及参照类别。
- 协变量均值
- 每个设计列(包括虚拟变量列)的(加权)均值;生存表即在这些取值处计算。
- 生存表
- 每个事件时间的基线累积风险(所有协变量为 0)、协变量均值处的生存率及其标准误,以及均值处的累积风险;每层各一组。
前提假设与注意事项
- 比例风险:每个协变量在所有时间点都使风险乘以相同的倍数。
- 连续协变量对对数风险呈线性作用;若不合理,应进行变换或分类。
- 在给定协变量的条件下删失为非信息性删失,且观测相互独立。
- 事件数足够,且没有协变量能完全分隔事件时间;否则系数会发散,面板会给出提示。
结果如何验证
在含一个连续协变量和一个三水平协变量的 150 例参考数据上,无论是否加权,系数与标准误均与 statsmodels 0.15 PHReg 及 R survival 3.8.12 coxph(Breslow 结处理)一致到 1e-7;-2 对数似然与得分卡方一致到 1e-9;协变量均值处的 Kalbfleisch-Prentice 生存率及其标准误与 R 一致到 1e-7;分层模型与使用 strata() 的 coxph 一致到 1e-7。采用条件、LR 与 Wald 剔除的向前、向后逐步选择与基于 R survival 独立编写的逐步实现逐步一致(统计量误差在 1e-6 以内),时依协变量 x·ln(T_) 与使用 tt() 的 coxph 一致到 1e-6。 验证证据(英文)
常见问题
如何解读 Exp(B)?
Exp(B) 即风险比。处理指示变量的 Exp(B) = 0.70 表示在其他协变量不变时,处理组在任一时点的瞬时风险比参照类别低 30%。对于连续协变量,它是每增加一个单位的风险比。
哪个类别是参照?
按排序后的最后一个类别,与 SPSS 默认的指示编码一致。如需其他参照类别,请先重新编码变量;“方程中的变量”表下的脚注列出了编码方式。
应看得分检验还是似然比检验?
两者都检验全部系数,通常结论一致。小样本时一般优先采用似然比检验(相对上一步的变化);得分检验不需要迭代。
如何处理结值事件时间?
采用 SPSS 所用的 Breslow 近似。结值很多时,Efron 法的估计会略有不同,因此结果可能与默认使用 Efron 的软件略有差异。
数据会被发送到别处吗?
不会。模型在浏览器中拟合,适合必须保留在本机上的临床或客户数据。
相关过程
- Kaplan-Meier 生存分析KM
- 寿命表SURVIVAL
- 复杂抽样 Cox 回归CSCOXREG
- 二元 Logistic 回归LOGISTIC REGRESSION
- 广义线性模型GENLIN
- 免费的 SPSS 替代工具:支持范围
用你自己的数据运行 Cox 回归
打开 .sav、CSV 或 Excel 文件。所有计算都在浏览器中完成,文件不会离开你的设备。
SPSS 是 IBM Corp. 的商标。Crosstabs 与 IBM 无关联,也未获得 IBM 认可。