system: [SYSTEM]

user: 3. 结果分析 采用1:1最近邻匹配 [ 14 ] [ 15 ] (Nearest Neighbor Matching)，匹配容差为0.2。两组的比较运用t检验，若呈现非高斯分布(Non-Gaussian distribution)则一般采用Wilcoxon秩和检验 [ 16 ]。匹配前，ICU患者的生存时间的中位数(即实际入ICU停留时间)大约为40小时(图1)；其中的PTCA手术组中位生存期约为42小时，而非手术组的中位生存期约为25小时。Kaplan-Meier生存曲线分析显示两组之间的差异具有统计学意义，非手术组的ICU患者实际入ICU停留时间短于PTCA手术组的ICU患者(P < 0.001) (图2)。 图1. 匹配前ICU患者生存曲线 图2. 匹配前PTCA手术组与非手术组ICU患者生存曲线 在进行进一步的分析之前，我们先比较两组基线特征，发现PTCA手术组和非手术组两组的专科医生(Physician Speciality)、医师干预类别(Physician Intervention Category)、急性生理学评分(Acute Physiology Score)、阿帕奇分数(Apache Score)、阿帕奇评分算法版本(Predicted ICU Mortality)、ICU预测停留时间(Predicted ICU Los)、ICU实际停留时间(Actual ICU Los)、医院预测死亡率(Predicted Hospital Mortality)、医院实际死亡情况(Actual Hospital Mortality)、住院预测时间(Predicted Hospital Los)、住院实际时间(Actual Hospital Los)、是否有术前性心肌梗塞(preopmi)、是否进行过术前心脏导管检查(preopcardiaccath)、未删减ICU实际停留时间(unabridged unit los)、为删减住院实际时间(unabridged hosp los)这些协变量均有统计学上的差异(均有P < 0.001)；为了匹配两组基线在特征上的差异化，本研究运用倾向性得分匹配(PSM)的匹配方法来平衡两组基线特征偏移，在进行匹配的过程中我们按照1:1的比例匹配后总共筛选出17，686例患者，其中PTCA手术组和非手术组各8843例。我们在此匹配的基础上，比较匹配后的两组患者， χ 2 的检验表明两组的ICU患者在临床病理特征上的差异都明显减小(表1) (图3~4)。 Table 1 Distribution profiles of the clinicopathologic factors of patients in the PTCA operation group and the non-operation group before and after PSM matching (n% 临床病理参数 匹配前 P-Value 匹配后 -Value PTCA手术组 (n = 8862) 非手术组 (n = 139，670) PTCA手术组 (n = 8843) 非手术组 (n = 8843) Physician speciality cardiology 3650 (41.1) 6399 (4.5) <0.001 2586 (29.2) 2121 (23.9) 0.356 Specialty Not Specified 1081 (12.2) 16，661 (11.9) 1178 (13.3) 1706 (19.2) internal medicine 972 (10.9) 20，488 (14.6) 1056 (11.9) 920 (10.4) hospitalist 932 (10.5) 19，581 (14.0) 1178 (13.3) 920 (10.4) others 2227 (25.3) 76，541 (55.0) 2845 (32.3) 2958 (33.7) Physician intervention category I 1661 (18.7) 24，269 (17.3) <0.001 1746 (19.7) 1512 (17.0) 0.942 II 2072 (23.3) 29，916 (21.4) 2445 (27.6) 1520 (17.1) III 1584 (17.8) 28，227 (20.2) 2118 (23.6) 1434 (16.2) IV 190 (2.1) 3654 (2.6) 148 (1.6) 212 (2.3) Unknown 3355 (38.1) 53，559 (38.5) 2386 (27.5) 4165 (47.4) Acute physiology score <38 5338 (60.2) 70，760 (50.7) <0.001 5493 (62.1) 5412 (61.2) 0.153 >38 3524 (39.8) 68，910 (49.3) 3350 (37.9) 3430 (38.8) Apache score <51 5178 (58.4) 71，178 (51.0) <0.001 5350 (60.5) 5272 (59.6) 0.201 >51 3684 (41.6) 68，492 (49.0) 3493 (39.5) 3571 (40.4) Apache version IV 4361 (49.2) 68，581 (49.1) 0.844 4359 (49.2) 4262 (48.1) 0.343 IVa 3684 (41.6) 68，492 (49.0) 4484 (50.8) 4581 (51.9) Predicted ICU mortality <0.023 4656 (52.5) 68，880 (49.3) <0.001 4969 (56.2) 4930 (55.8) 0.084 >0.023 4206 (47.5) 70，790 (50.7) 3874 (43.8) 3913 (44.2) Predicted ICU los <3.2 5848 (66.0) 67，862 (48.6) <0.001 5735 (64.9) 5739 (64.9) 0.005 >3.2 3014 (34.0) 71，808 (51.4) 3108 (35.1) 3104 (35.1) Actual ICU los <1.8 5014 (56.6) 68，775 (49.2) <0.001 5106 (57.7) 5046 (57.1) 0.234 >1.8 3848 (43.4) 70，895 (50.8) 3737 (42.3) 3797 (42.9) Predicted hospital mortality <0.05 4825 (54.4) 71，662 (51.3) <0.001 5019 (56.8) 5043 (57.0) 0.590 >0.05 4037 (45.6) 68，008 (48.7) 3824 (43.2) 3800 (43.0) Actual hospital mortality ALIVE 7410 (83.6) 127，385 (91.2) <0.001 7493 (84.7) 7404 (83.7) 0.095 EXPIRED 1452 (16.4) 12，285 (8.8) 1350 (15.3) 1439 (16.3) Predicted hospital los <8.8 6674 (75.3) 67，023 (48.0) <0.001 6603 (74.7) 3945 (44.6) 0.883 >8.8 2188 (24.7) 72，647 (52.0) 2240 (25.3) 4898 (55.4) Actual hospital los <5.5 6107 (68.9) 68，051 (48.7) <0.001 6071 (68.7) 5974 (67.6) 0.821 >5.5 2755 (31.1) 71，619 (51.3) 2772 (31.3) 2869 (32.4) preopmi 0 8862 (100) 139，246 (99.7) <0.001 8843 (100) 8842 (99.9) 0.317 1 0 (0) 424 (0.3) 0 (0) 1 (0.1) preopcardiaccath 0 8862 (100) 138，533 (99.2) <0.001 8843 (100) 8843 (100) - 1 0 (0) 1137 (0.8) 0 (0) 0 (0) Unabridged unit los <1.8 5014 (56.6) 68，775 (49.2) <0.001 5106 (57.7) 5046 (57.1) 0.234 >1.8 3848 (43.4) 70，895 (50.8) 3737 (42.3) 3797 (42.9) Unabridged hosp los <5.5 6107 (68.9) 68，051 (48.7) <0.001 6071 (68.7) 5974 (67.6) 0.478 >5.5 2755 (31.1) 71，619 (51.3) 2772 (31.3) 2869 (32.4) 表1. PSM匹配前后PTCA手术组和非手术组ICU患者临床病理资料分布特征(n%) 图3. PTCA手术组与非手术组两组的ICU患者PSM前后倾向性评分分布的特征 图4. 经过PSM筛选及去除的ICU患者的分布特征 在经过倾向性得分匹配平衡各协变量之间的差异后，原则上我们就可以对PTCA手术组ICU患者和非手术组ICU患者的生存状况之间是否有显著性差异直接检验了。常规的方法两组的比较运用t检验，但是这种方法在数据呈非高斯分布(即非正态分布)的情况下检验效果较差，且由检验结果(表2)我们可以看出匹配后我们的PTCA手术组和非手术组两组数据都是呈现的非高斯分布。故我们采用大样本性质较好且不要求检验数据呈高斯分布前提的Wilcoxon秩和检验来判断两组之间是否有显著性差异。显然，Wilcoxon秩和检验的结果(表3)告诉我们PTCA手术组ICU患者和非手术组ICU患者的生存状况之间是明显的差异，PTCA手术组ICU患者生存率高于非手术组，即非手术组预后的状况差于PTCA手术组；其次，我们很清晰的看到ICU患者的生存期中位数(实际入ICU停留时间)为19小时左右(图5)；其中PTCA手术组的ICU患者生存期中位数为22小时左右，而非手术组的ICU患者生存期中位数为16小时左右(图6)。结合Kamplan-Meier生存曲线分析(图6)、加速失效模型 [ 17 ] [ 18 ] [ 19 ] (accelerated failure time， AFT) (表4)显示两组生存差异具有统计学意义，PTCA手术组实际入ICU停留时间长于非手术组。 Table 2 Normal test of PTCA operation group and non operation group after matchin 检验方法：One-sample Kolmogorov-Smirnov test Data：PTCA手术组 Data：非手术组 检验统计量D = 0.5 P-Value < 0.001 检验统计量D = 0.5 P-Value < 0.001 Alternative hypothesis: two-sided Alternative hypothesis: two-sided 表2. 匹配后PTCA手术组与非手术组各自的正态性检验 Table 3 Wilcoxon rank sum test of PTCA operation group and non operation group after matchin 检验方法：Wilcoxon rank sum test with continuity correction Data：PTCA手术组与非手术组 检验统计量W = 38，126，595 P-Value < 0.001 alternative hypothesis: true location shift is not equal to 0 表3. 匹配后PTCA手术组与非手术组的Wilcoxon秩和检验 图5. 匹配后ICU患者生存曲线 图6. 匹配后PTCA手术组与非手术组ICU患者生存曲线 Table 4 Survival time of AFT mode 变量 回归系数 标准误 P-Value 备注 截距项 3.227 0.069 P < 0.001 PTCA with 24 h (非手术组) −0.189 0.015 P < 0.001 基准组：PTCA手术组 表4. 生存时长AFT模型