输卵管注药绝育术效果影响因素的logistic回归分析
输卵管注药绝育术效果影响因素的logistic回归分析
郭秀花 刘云嵘 李东
关键词：输卵管注药绝育术 logistic回归分析 一、资料与方法
1.调查对象：输卵管注药绝育术多中心临床研究在1986年4月至1991年6月开展。被接纳进入该研究的共有1 705名受术合格对象。其纳入条件为术时年龄20～35岁，已有两个或两个以上孩子的健康妇女，自愿接受本术，无禁忌者。术时按随机原则将复方苯酚糊剂(PAP)和显影苯酚胶浆(PM)分配给受术对象，PAP组871名，PM组834名。施术采用经宫颈的输卵管间质部插管定量注药法，往受术者两侧输卵管内分别注入绝育药物0.08 ml。手术后7～14 d、3个月、1年、2年及5年进行随访调查。最后共失访27例，更换其他避孕方法6例，因病死亡2例，因子宫肌瘤行全摘术2例。对于以上37例受术者最后终止日期不明。其余1 668例受术者或有明确诊断的妊娠(施术失败)时间或有5年随访时未发生妊娠(施术成功)的完整资料，应用logistic回归模型探讨施术效果的影响因素[1]。
, 百拇医药
2.方法：将1 705例随访资料用Visual Foxpro 5.0建立数据库，并对全部资料进行了逻辑审校与修正。对其中1 668例符合条件的数据采用SAS 6.04版统计软件进行统计分析。初选施术可能的影响因素12个，先进行单因素logistic回归分析，再用多因素非条件logistic回归来提取主要因素，并计算比数比。
二、结果
设应变量(Y)为注药绝育术结局(成功＝0，失败＝1)；自变量的设置和数量化值方法见表1。
1.单因素分析：首先对初选的12项指标按等级顺序赋予初值，然后以注药绝育术结局变量Y分别与12个因素进行单因素的logistic回归分析。初选的12个因素的赋值方法及单因素logistic回归分析的主要结果见表1。
2.多因素非条件logistic回归：从表1中可知，P<0.10的因素有5个：绝育药物种类X1、年龄X3、受术时机X6、术时有否取环X11和术后盆腔X线片结果X12；从专业上选取：术时宫腔深变X8。将这6个因素与Y一起，采用多因素非条件logistic回归逐步向前剔除法和用逐步向后剔除法。
, 百拇医药
用逐步向前剔除法按P＝0.05为入选界值时，选入模型的因素有3个，按入选的先后顺序依次是：术后盆腔X线片X12结果(OR＝1.910 8)、绝育药物种类X1(OR＝0.301 4)和年龄X3(OR＝1.669 1)，详细结果见表2。
3个因素联合作用的影响显著(似然比χ2＝129.648，P＝0.000 1)，模型的预测能力是：一致性(concordant)＝73.2%，不和谐(不一致)(discordant)＝22.7%。
用逐步向后剔除法，得到回代结果的分类表，当割点概率为0.5时，灵敏度为91.5%、专一性为5.7%、假正比率为8.0%、假负比率52.9%。
三、讨论
1.输卵管注药绝育术是经宫颈注入药物闭塞输卵管导致女性绝育的一种新技术,其先进性已有报道。本文的目的是利用logistic回归模型分析探讨影响手术成功的主要因素。结果表明：在影响该手术成功方面，术后盆腔X线片结果，即绝育药物在输卵管内的充盈长度与本术效果有关，充盈长度每增加1个等级，成功可能性增加1.910 8倍；显影苯酚胶浆是复方苯酚糊剂节育效果的0.301 4倍，说明复方苯酚糊剂是一种节育效果较好的输卵管绝育粘堵剂；在20～35岁年龄段内施行本术时，年龄每增加1岁，成功可能性增加1.669 1倍。
, http://www.100md.com
表1 自变量各因素的赋值方法及单因素统计分析结果
变量名称
因素名称(赋值方法)
χ2值
P值
X1
绝育药物种类(复方苯酚糊剂0，显影苯酚胶浆1)
26.913 6
0.000 1
X2
文化程度(文盲 0，小学 1，初中 2，高中及以上3)
, 百拇医药
0.023 0
0.879 4
X3
年龄(20岁～ 0，25岁～ 1，≥30岁 2)
14.658 0
0.000 1
X4
分娩次数(≤2次 0，3次 1，≥4次 2)
1.236 9
0.266 1
X5
人工流产次数(0次 0，1次 1，2次 2，≥3次 3)
, 百拇医药
0.034 5
0.852 7
X6
受术时机(月经前半周期 0，月经后半周期 1，哺乳期月经不规则 2，哺乳期4～6个月 3，哺乳6个月以上 4)
4.979 5
0.025 6
X7
术时子宫位置(前倾 0，中位 1，后倾 2)
0.065 1
0.798 6
X8
, http://www.100md.com
术时宫腔深变(50 cm～ 0，60 cm～ 1，70 cm～ 2，≥80 cm 3)
2.037 0
0.153 5
X9
扩宫(无 0，有 1)
0.857 3
0.354 5
X10
术时插管有无困难(顺利 0，一般 1，困难 2)
0.687 8
0.406 9
, http://www.100md.com
X11
术时取环(无 0，有 1)
3.861 7
0.049 4
X12
术后盆腔X线片结果(未达峡部 0，峡部0.5～1.4 cm 1，峡部1.5～1.9 cm 2，峡部2.0～3.9 cm 3，峡部≥4.0 cm 4，壶腹部 5)
79.687 8
0.000 1
表2 3个入选因素logistic回归分析结果
名称
, http://www.100md.com
β
标准误
χ2值
P值
标准值
OR值
OR的95%CI
截距
0.638 4
0.266 9
5.722 1
0.016 8
-
, 百拇医药
-
-
X1
-1.199 3
0.203 2
34.842 0
0.000 1
-0.330 638
0.301 4
0.202 4～0.448 9
X3
0.512 3
, 百拇医药
0.125 4
16.693 7
0.000 1
0.205 713
1.669 1
1.305 4～2.134 2
X12
0.647 5
0.070 3
84.914 1
0.000 1
0.416 360
, 百拇医药
1.910 8
1.664 8～2.193 0
2.在用SAS软件进行logistic回归分析时值得注意的问题是：关于响应变量Y的赋值方法与其他许多教材和统计软件不同。SAS软件将响应变量Y某事件发生赋值为0、某事件不发生赋值为1，这是因为SAS软件的logistic回归分析的模型与其他软件(如STATA)不同，致使输出结果中偏回归系数估计的符号相反。
(本文所用数据来自国家计划生育委员会项目“输卵管注药绝育术有效性安全性五年随访研究”；完成本项目和数据收集的主要研究人员有：吴裕浩、陈学煌、黄真嘉、徐铎、蔡光宗等，在此一并致谢)基金项目：国家计划生育委员会资助项目(91025)
作者单位：郭秀花(解放军北京军医学院数学教研室 100071)
李东(计划生育办公室)
刘云嵘(国家计划生育委员会科学技术研究所流行病学与社会医学研究室)
参考文献
1，郭秀花.未婚青年手淫情况的Logistic回归分析.中国公共卫生，1999，15∶243-244., 百拇医药
郭秀花 刘云嵘 李东
关键词：输卵管注药绝育术 logistic回归分析 一、资料与方法
1.调查对象：输卵管注药绝育术多中心临床研究在1986年4月至1991年6月开展。被接纳进入该研究的共有1 705名受术合格对象。其纳入条件为术时年龄20～35岁，已有两个或两个以上孩子的健康妇女，自愿接受本术，无禁忌者。术时按随机原则将复方苯酚糊剂(PAP)和显影苯酚胶浆(PM)分配给受术对象，PAP组871名，PM组834名。施术采用经宫颈的输卵管间质部插管定量注药法，往受术者两侧输卵管内分别注入绝育药物0.08 ml。手术后7～14 d、3个月、1年、2年及5年进行随访调查。最后共失访27例，更换其他避孕方法6例，因病死亡2例，因子宫肌瘤行全摘术2例。对于以上37例受术者最后终止日期不明。其余1 668例受术者或有明确诊断的妊娠(施术失败)时间或有5年随访时未发生妊娠(施术成功)的完整资料，应用logistic回归模型探讨施术效果的影响因素[1]。
, 百拇医药
2.方法：将1 705例随访资料用Visual Foxpro 5.0建立数据库，并对全部资料进行了逻辑审校与修正。对其中1 668例符合条件的数据采用SAS 6.04版统计软件进行统计分析。初选施术可能的影响因素12个，先进行单因素logistic回归分析，再用多因素非条件logistic回归来提取主要因素，并计算比数比。
二、结果
设应变量(Y)为注药绝育术结局(成功＝0，失败＝1)；自变量的设置和数量化值方法见表1。
1.单因素分析：首先对初选的12项指标按等级顺序赋予初值，然后以注药绝育术结局变量Y分别与12个因素进行单因素的logistic回归分析。初选的12个因素的赋值方法及单因素logistic回归分析的主要结果见表1。
2.多因素非条件logistic回归：从表1中可知，P<0.10的因素有5个：绝育药物种类X1、年龄X3、受术时机X6、术时有否取环X11和术后盆腔X线片结果X12；从专业上选取：术时宫腔深变X8。将这6个因素与Y一起，采用多因素非条件logistic回归逐步向前剔除法和用逐步向后剔除法。
, 百拇医药
用逐步向前剔除法按P＝0.05为入选界值时，选入模型的因素有3个，按入选的先后顺序依次是：术后盆腔X线片X12结果(OR＝1.910 8)、绝育药物种类X1(OR＝0.301 4)和年龄X3(OR＝1.669 1)，详细结果见表2。
3个因素联合作用的影响显著(似然比χ2＝129.648，P＝0.000 1)，模型的预测能力是：一致性(concordant)＝73.2%，不和谐(不一致)(discordant)＝22.7%。
用逐步向后剔除法，得到回代结果的分类表，当割点概率为0.5时，灵敏度为91.5%、专一性为5.7%、假正比率为8.0%、假负比率52.9%。
三、讨论
1.输卵管注药绝育术是经宫颈注入药物闭塞输卵管导致女性绝育的一种新技术,其先进性已有报道。本文的目的是利用logistic回归模型分析探讨影响手术成功的主要因素。结果表明：在影响该手术成功方面，术后盆腔X线片结果，即绝育药物在输卵管内的充盈长度与本术效果有关，充盈长度每增加1个等级，成功可能性增加1.910 8倍；显影苯酚胶浆是复方苯酚糊剂节育效果的0.301 4倍，说明复方苯酚糊剂是一种节育效果较好的输卵管绝育粘堵剂；在20～35岁年龄段内施行本术时，年龄每增加1岁，成功可能性增加1.669 1倍。
, http://www.100md.com
表1 自变量各因素的赋值方法及单因素统计分析结果
变量名称
因素名称(赋值方法)
χ2值
P值
X1
绝育药物种类(复方苯酚糊剂0，显影苯酚胶浆1)
26.913 6
0.000 1
X2
文化程度(文盲 0，小学 1，初中 2，高中及以上3)
, 百拇医药
0.023 0
0.879 4
X3
年龄(20岁～ 0，25岁～ 1，≥30岁 2)
14.658 0
0.000 1
X4
分娩次数(≤2次 0，3次 1，≥4次 2)
1.236 9
0.266 1
X5
人工流产次数(0次 0，1次 1，2次 2，≥3次 3)
, 百拇医药
0.034 5
0.852 7
X6
受术时机(月经前半周期 0，月经后半周期 1，哺乳期月经不规则 2，哺乳期4～6个月 3，哺乳6个月以上 4)
4.979 5
0.025 6
X7
术时子宫位置(前倾 0，中位 1，后倾 2)
0.065 1
0.798 6
X8
, http://www.100md.com
术时宫腔深变(50 cm～ 0，60 cm～ 1，70 cm～ 2，≥80 cm 3)
2.037 0
0.153 5
X9
扩宫(无 0，有 1)
0.857 3
0.354 5
X10
术时插管有无困难(顺利 0，一般 1，困难 2)
0.687 8
0.406 9
, http://www.100md.com
X11
术时取环(无 0，有 1)
3.861 7
0.049 4
X12
术后盆腔X线片结果(未达峡部 0，峡部0.5～1.4 cm 1，峡部1.5～1.9 cm 2，峡部2.0～3.9 cm 3，峡部≥4.0 cm 4，壶腹部 5)
79.687 8
0.000 1
表2 3个入选因素logistic回归分析结果
名称
, http://www.100md.com
β
标准误
χ2值
P值
标准值
OR值
OR的95%CI
截距
0.638 4
0.266 9
5.722 1
0.016 8
-
, 百拇医药
-
-
X1
-1.199 3
0.203 2
34.842 0
0.000 1
-0.330 638
0.301 4
0.202 4～0.448 9
X3
0.512 3
, 百拇医药
0.125 4
16.693 7
0.000 1
0.205 713
1.669 1
1.305 4～2.134 2
X12
0.647 5
0.070 3
84.914 1
0.000 1
0.416 360
, 百拇医药
1.910 8
1.664 8～2.193 0
2.在用SAS软件进行logistic回归分析时值得注意的问题是：关于响应变量Y的赋值方法与其他许多教材和统计软件不同。SAS软件将响应变量Y某事件发生赋值为0、某事件不发生赋值为1，这是因为SAS软件的logistic回归分析的模型与其他软件(如STATA)不同，致使输出结果中偏回归系数估计的符号相反。
(本文所用数据来自国家计划生育委员会项目“输卵管注药绝育术有效性安全性五年随访研究”；完成本项目和数据收集的主要研究人员有：吴裕浩、陈学煌、黄真嘉、徐铎、蔡光宗等，在此一并致谢)基金项目：国家计划生育委员会资助项目(91025)
作者单位：郭秀花(解放军北京军医学院数学教研室 100071)
李东(计划生育办公室)
刘云嵘(国家计划生育委员会科学技术研究所流行病学与社会医学研究室)
参考文献
1，郭秀花.未婚青年手淫情况的Logistic回归分析.中国公共卫生，1999，15∶243-244., 百拇医药