中年知识分子心理健康影响因素的广义相对危险度模型
作者：李新华 邓冰 朱焱
单位：550004 贵阳医学院预防医学系
关键词：心理健康；相对危险度；中年知识分子
中年知识分子心理健康影响因素的 广义相对危险度模型 【摘要】 目的 研究影响中年知识分子心理健康危险因素的作用及模式，方法 对1503名中年知识分子进行SCL-90测试和影响因素的调查，用广义相对危险度模型分析危险因素的作用大小及相互作用形式，结果 约 6.9% 的中年知识分子(男性为 6.1% ，女性为 7.8%)有心理健康问题，按作用大小危险因素依次是离婚或丧偶、接触毒物或噪声、缺少业余爱好、工作时间延长、患慢性病。多因素模型趋向于超相乘结构，危险因素间联合作用趋向于超相乘或中间结构。结论 广义相对危险度模型适用于心理学领域的流行病学研究。
Application of relative risk models in the study of psychological health on middle-age intellectual
, 百拇医药
Li Xinhua, Deng Bing , Zhu Yan. Department of Preventive Medicine , Guiyang Medical College , Guiyang 550004
【Abstract】 Objective To study the effects and models of risk factors of psychological health on middle-age intellectual. Method The psychological health of 1503 middle-age intellectual was measured by SCL-90. The effects and interaction modes of risk factors were analyzed with general relative risk model. Results The rate of psychological disorder was 6.9% in 1503 intellectual. The morbidity was 6.1% in male and 7.8% in female. The results suggested that the major risk factors were, as listed in a order of their effect on psychological disorder , divorcement or lose of spouse,contact of poison and noisiness, lack of avocation, extend of working hours, having chronic diseases. The multivariate model trended towards supermultiplicative structure. The allied effects between risk factors trended towards supermultiplicative or midway structure. Conclusion The general relative risk model is applicable to the epidemiological research in psychology.
, 百拇医药
【Key words】 Psychological health Relative risk Middle-age intellectual
随着医学模式的发展，以特殊人群为研究对象来研究疾病和健康问题，寻求增进人类健康的途径，已逐渐成为现代医学发展的一个方向。中年知识分子处于身体机能和心理状态的特殊阶段，并且担负着工作和家庭的双重重担，各种压力会影响其心理健康，为了解中年知识分子心理健康状况，我们于1998年4月进行了本次研究，并采用广义相对危险度模型来研究危险因素的作用及模式。
对象和方法
一、对象 采用分层整群抽样方法抽取了贵阳地区省市级企业、事业、教育、卫生、科研单位中1503名40～59岁的中年知识分子，文化程度为中专以上，男性占57.6%，女性占42.4%，平均年龄48.74岁，从事专业技术职务者1251人，专任行政职务者252人。
, 百拇医药
二、 方法
(一) 调查方法 全部对象均用Derogatis 氏90项症状自评量表(SCL-90)[1]进行心理状况自评，计算10个因子分和总均分，根据总均分大于等于 3分(即中等严重度)判定对象存在心理问题。同时按统一的项目填写调查表，内容包括生理健康、社会经济、环境因素、行为方式等方面。
(二)信度分析 按SCL-90规定将90个项目计算成10个因子分，用主成分分析法计算出最大的特征值Λ，然后计算可靠性系数
根据Q值大小即可评价本次调查的信度[2]。随机抽取100例调查对象在半年后用SCL-90量表再测，计算两次测定结果总均分及10个因子分的相关系数，以评价再测信度。
, 百拇医药
(三)广义相对危险度模型
流行病学病例对照研究中分析多个危险因素的作用常用对数线性模型、logistic 回归模型和比例危险模型，由这些模型分析的结果获得的一个重要指标就是相对危险度，故称为相对危险度模型。这些模型的假设条件是不同危险因素之间的联合作用呈相乘关系,一些研究发现乘法结构不一定能最正确地说明资料的特征，也不应不经过仔细的检查就例行地接受乘法结构的结果，[3]统计学家以具有某种协变量特征的研究人群的危险率与其处于基准状态时的危险率之比为出发点，研究该比值与个体所具有的协变量特征之间存在的函数关系，提出了一些扩展的相对危险度模型，[3，4，5]Breslow提出的广义相对危险度模型函数是定义R(x)为暴露变量X的相对危险度(Relative Risk)，将其作为一个函数，建立通用模型函数式为:
Ln R(Xβ;λ) = [(1 + Xβ)λ - 1 ]/ λ
, http://www.100md.com
式中x为自变量矢量，β为回归系数矢量，λ为尺度参数，它决定函数式的形状，λ= 0 时
Ln R(Xβ;λ) = Ln (1 + Xβ)
为相加(additive)模型，即两变量之间联合作用的危险度为R(A,B) = R(A) + R(B) - 1 ,当λ= 1 时
Ln R(Xβ;λ) =[(1 + Xβ) - 1]
为相乘模型(logistic模型或对数线性模型)，R(A,B) = R(A)R(B) , λ< 0 为次相加(subadditive)模型， R(A,B) < R(A) + R(B) - 1 ,λ> 1为超相乘(supermultiplicative)模型，R(A,B) > R(A)R(B)，λ介于0和1之间时为中间模型(超相加或次相乘模型)。参数的估计用加权最小二乘法在 GLIM 软件包上进行，通过λ值的改变，用偏差度(Deviance)的大小衡量模型的拟合效果，用卡方检验模型的拟合优度，以 Wald 方法检验回归系数的显著性，选择最佳模型。
, http://www.100md.com
结 果
本次研究分析的可靠性系数Q为0.825，总均分再测相关系数为0.6106，各因子分的再测相关系数见表1，P<0.01，表明信度可靠。按总均分大于等于3分判定，约 6.9% 的中年知识分子(男性为 6.1% ，女性为 7.8%)有心理健康问题。研究因素赋值情况和单因素分析结果如表2。
表1 10个因子分的再测相关系数
因子
1
2
3
4
5
6
, 百拇医药
7
8
9
10
相关系数
0.715
0.619
0.400
0.592
0.493
0.347
0.314
0.426
, http://www.100md.com
0.398
0.679
表2 分析变量赋值及单因素分析结果
暴露变量
赋 值
χ2
P
R(x)
X1
性别
1男 2 女
1.72
>0.05
, 百拇医药
X2
年龄组
140～；2 45～；3 50～；4 55～59
10.71
<0.05
0.81
X3
文化程度
1大学本科以下 2 本科及以上
0.75
>0.05
X4
, http://www.100md.com
家庭经济状况
0能维持开支 1 不能维持开支
4.09
<0.05
1.66
X5
离婚或丧偶
0否 1 是
22.35
<0.01
3.46
X6
吸烟
, 百拇医药
0否 1 是
0.02
>0.05
X7
经常过量饮酒
0否 1 是
11.34
<0.01
2.48
X8
体育锻炼
1每天坚持 2 偶尔或无
3.68
, http://www.100md.com
>0.05
X9
业余爱好
0有 1 无
5.64
<0.05
1.63
X10
重视健康
0不能顾及或无所谓 1 是
5.62
<0.05
0.61
, 百拇医药
x11
慢性病
0无 1 有
5.03
<0.05
1.58
x12
睡眠不足
06小时及以上 1 6小时以下
6.06
<0.05
1.65
x13
, 百拇医药
工作时间延长
0无 1 有
7.95
<0.01
1.77
x14
接触毒物或噪声
0无 1 有
10.12
<0.01
1.93
多因素广义相对危险度模型分析结果如表3。为了探讨危险因素间联合作用的模式，对筛选出的因素进行两两间作用关系的分析，以λ与偏差度之间的动态变化图来判断实际数据趋向于何种模型，结果如表4。
, http://www.100md.com
表3 多因素分析的模型比较
因素
λ=0.0
λ=1.0
λ=1.2
P
R(x)
β
Sβ
β
Sβ
β
Sβ
, 百拇医药
标化β
X5
7.305
4.200
1.183
0.282
1.004
0.230
4.365
<0.01
3.044
X9
1.533
, http://www.100md.com
1.143
0.506
0.211
0.447
0.180
2.483
<0.05
2.082
X11
1.933
1.267
0.413
0.209
, 百拇医药
0.336
0.160
2.288
<0.05
1.941
X13
0.477
0.788
0.456
0.212
0.416
0.181
2.298
, 百拇医药
<0.05
2.041
X14
1.827
1.205
0.587
0.215
0.513
0.185
2.773
<0.01
2.173
偏差度
, 百拇医药
34.80
32.55
32.51
卡方值
31.35
28.44
28.56
表4 危险因素间作用关系的分析(R(x))
暴露状态
X5,X9
X5,X11
X5,X13
, http://www.100md.com
X5,X14
X9,X11
X9,X13
X9,X14
X11,X13
X11,X14
X13,X14
- -
1.00
1.00
1.00
1.00
1.00
, 百拇医药
1.00
1.00
1.00
1.00
1.00
- +
1.55
1.74
1.49
2.51
2.22
2.11
2.13
, http://www.100md.com
2.69
5.40
2.24
+ -
3.01
4.05
2.05
5.14
2.28
2.08
2.08
2.33
3.64
, 百拇医药
2.23
+ +
6.67
5.25
7.37
8.76
3.21
2.24
2.26
4.62
15.89
2.28
偏差度
, 百拇医药
0.0004
0.379
0.0024
0.0053
1.4500
0.0010
0.0020
0.0003
0.0737
0.0034
λ
2.2
0.25
, 百拇医药
2.4
1.4
1.4
2.6
2.6
0.8
0.6
4.2
讨 论
本研究定量分析了生理健康、社会经济、环境因素、行为生活方式对心理健康的影响，单因素的分析显示，家庭经济状况差、经常过量饮酒、离婚或丧偶、缺少业余爱好、患慢性病、睡眠不足、工作时间延长、接触毒物或噪声是影响中年知识分子心理健康的危险因素(相对危险度大于1.0)，年龄组和重视健康则是有利因素(相对危险度小于1.0)，年龄组的这一特点与全国情况相吻合[6]。
, 百拇医药
为了控制混杂因素的作用，应用多因素广义相对危险度模型筛选因素，结果显示以超相乘模型较为适宜，根据标化回归系数判定危险因素的的作用大小依次为离婚或丧偶、接触毒物或噪声、缺少业余爱好、工作时间延长、患慢性病，在固定其他因素作用不变的情况下，其相对危险度依次为3.044、2.173、2.802、2.041和1.941。
危险因素间的相互作用分析提示: 离婚或丧偶与缺少业余爱好、工作时间延长、接触毒物或噪声；缺少业余爱好与患慢性病、工作时间延长、接触毒物或噪声；工作时间延长与接触毒物或噪声之间两两联合作用的形式趋向于超相乘结构。离婚或丧偶与患慢性病；患慢性病与工作时间延长、接触毒物或噪声之间两两联合作用的形式趋向于中间结构模型。
目前国内流行病学调查常用logistic 回归模型分析危险因素的作用，但假设条件是各因素之间的相互作用呈相乘关系，实际工作中的数据不一定都符合此假设，若不符合此假设，模型便会失效，如有的情况下各因素之间的相互作用呈相加关系，应该使用相加模型[7]，或事先不清楚数据适合于何种模型。心理学领域的流行病学研究涉及的因素繁多，因素之间的关系错综复杂，一概采用logistic 回归模型来分析不一定妥当。Breslow提出的广义相对危险度模型函数，通过不同指数λ的幂变换，自变量与相对危险度之间的连接形式范围从次相加结构变化到超相乘结构，并可以从偏差度的图分析方法及拟合优度检验来较为客观地选择最佳模型，因而适用于心理学领域的流行病学研究。
, 百拇医药
*本研究为贵州省科委科学技术基金资助项目
参考文献
1 徐俊冕主编. 医学心理学. 第二版.上海:上海医科大学出版社，1996:183～190.
2 王庆石.统计指标导论.第一版.大连:东北财经大学出版社，1994:54～55.
3 Breslow，N.E.，and Storer，B.E.General relative risk functions for case-control studies. Am J Epidemiol， 1985，122(1)：149～162.
4 项永兵，高玉堂，金凡.相对危险度回归模型.中华流行病学杂志，1994，15(4):244～246.
5 赵宁，俞顺章，孙伟明. 相对危险度函数蔟模型在原发性肝癌危险因素研究中的应用. 中国卫生统计, 1995，12(4):13～17.
6 金华，吴文源，张明圆. 中国正常人SCL-90评定结果的初步分析.中国神经精神疾病杂志，1986，12(5)：260～262.
7 方积乾，陈和年.医学研究中logistic回归模型的正确应用(二).中国卫生统计, 1993，10(5):61～63.
(收稿日期：1999-04-19), 百拇医药
单位：550004 贵阳医学院预防医学系
关键词：心理健康；相对危险度；中年知识分子
中年知识分子心理健康影响因素的 广义相对危险度模型 【摘要】 目的 研究影响中年知识分子心理健康危险因素的作用及模式，方法 对1503名中年知识分子进行SCL-90测试和影响因素的调查，用广义相对危险度模型分析危险因素的作用大小及相互作用形式，结果 约 6.9% 的中年知识分子(男性为 6.1% ，女性为 7.8%)有心理健康问题，按作用大小危险因素依次是离婚或丧偶、接触毒物或噪声、缺少业余爱好、工作时间延长、患慢性病。多因素模型趋向于超相乘结构，危险因素间联合作用趋向于超相乘或中间结构。结论 广义相对危险度模型适用于心理学领域的流行病学研究。
Application of relative risk models in the study of psychological health on middle-age intellectual
, 百拇医药
Li Xinhua, Deng Bing , Zhu Yan. Department of Preventive Medicine , Guiyang Medical College , Guiyang 550004
【Abstract】 Objective To study the effects and models of risk factors of psychological health on middle-age intellectual. Method The psychological health of 1503 middle-age intellectual was measured by SCL-90. The effects and interaction modes of risk factors were analyzed with general relative risk model. Results The rate of psychological disorder was 6.9% in 1503 intellectual. The morbidity was 6.1% in male and 7.8% in female. The results suggested that the major risk factors were, as listed in a order of their effect on psychological disorder , divorcement or lose of spouse,contact of poison and noisiness, lack of avocation, extend of working hours, having chronic diseases. The multivariate model trended towards supermultiplicative structure. The allied effects between risk factors trended towards supermultiplicative or midway structure. Conclusion The general relative risk model is applicable to the epidemiological research in psychology.
, 百拇医药
【Key words】 Psychological health Relative risk Middle-age intellectual
随着医学模式的发展，以特殊人群为研究对象来研究疾病和健康问题，寻求增进人类健康的途径，已逐渐成为现代医学发展的一个方向。中年知识分子处于身体机能和心理状态的特殊阶段，并且担负着工作和家庭的双重重担，各种压力会影响其心理健康，为了解中年知识分子心理健康状况，我们于1998年4月进行了本次研究，并采用广义相对危险度模型来研究危险因素的作用及模式。
对象和方法
一、对象 采用分层整群抽样方法抽取了贵阳地区省市级企业、事业、教育、卫生、科研单位中1503名40～59岁的中年知识分子，文化程度为中专以上，男性占57.6%，女性占42.4%，平均年龄48.74岁，从事专业技术职务者1251人，专任行政职务者252人。
, 百拇医药
二、 方法
(一) 调查方法 全部对象均用Derogatis 氏90项症状自评量表(SCL-90)[1]进行心理状况自评，计算10个因子分和总均分，根据总均分大于等于 3分(即中等严重度)判定对象存在心理问题。同时按统一的项目填写调查表，内容包括生理健康、社会经济、环境因素、行为方式等方面。
(二)信度分析 按SCL-90规定将90个项目计算成10个因子分，用主成分分析法计算出最大的特征值Λ，然后计算可靠性系数
根据Q值大小即可评价本次调查的信度[2]。随机抽取100例调查对象在半年后用SCL-90量表再测，计算两次测定结果总均分及10个因子分的相关系数，以评价再测信度。, 百拇医药
(三)广义相对危险度模型
流行病学病例对照研究中分析多个危险因素的作用常用对数线性模型、logistic 回归模型和比例危险模型，由这些模型分析的结果获得的一个重要指标就是相对危险度，故称为相对危险度模型。这些模型的假设条件是不同危险因素之间的联合作用呈相乘关系,一些研究发现乘法结构不一定能最正确地说明资料的特征，也不应不经过仔细的检查就例行地接受乘法结构的结果，[3]统计学家以具有某种协变量特征的研究人群的危险率与其处于基准状态时的危险率之比为出发点，研究该比值与个体所具有的协变量特征之间存在的函数关系，提出了一些扩展的相对危险度模型，[3，4，5]Breslow提出的广义相对危险度模型函数是定义R(x)为暴露变量X的相对危险度(Relative Risk)，将其作为一个函数，建立通用模型函数式为:
Ln R(Xβ;λ) = [(1 + Xβ)λ - 1 ]/ λ
, http://www.100md.com
式中x为自变量矢量，β为回归系数矢量，λ为尺度参数，它决定函数式的形状，λ= 0 时
Ln R(Xβ;λ) = Ln (1 + Xβ)
为相加(additive)模型，即两变量之间联合作用的危险度为R(A,B) = R(A) + R(B) - 1 ,当λ= 1 时
Ln R(Xβ;λ) =[(1 + Xβ) - 1]
为相乘模型(logistic模型或对数线性模型)，R(A,B) = R(A)R(B) , λ< 0 为次相加(subadditive)模型， R(A,B) < R(A) + R(B) - 1 ,λ> 1为超相乘(supermultiplicative)模型，R(A,B) > R(A)R(B)，λ介于0和1之间时为中间模型(超相加或次相乘模型)。参数的估计用加权最小二乘法在 GLIM 软件包上进行，通过λ值的改变，用偏差度(Deviance)的大小衡量模型的拟合效果，用卡方检验模型的拟合优度，以 Wald 方法检验回归系数的显著性，选择最佳模型。
, http://www.100md.com
结 果
本次研究分析的可靠性系数Q为0.825，总均分再测相关系数为0.6106，各因子分的再测相关系数见表1，P<0.01，表明信度可靠。按总均分大于等于3分判定，约 6.9% 的中年知识分子(男性为 6.1% ，女性为 7.8%)有心理健康问题。研究因素赋值情况和单因素分析结果如表2。
表1 10个因子分的再测相关系数
因子
1
2
3
4
5
6
, 百拇医药
7
8
9
10
相关系数
0.715
0.619
0.400
0.592
0.493
0.347
0.314
0.426
, http://www.100md.com
0.398
0.679
表2 分析变量赋值及单因素分析结果
暴露变量
赋 值
χ2
P
R(x)
X1
性别
1男 2 女
1.72
>0.05
, 百拇医药
X2
年龄组
140～；2 45～；3 50～；4 55～59
10.71
<0.05
0.81
X3
文化程度
1大学本科以下 2 本科及以上
0.75
>0.05
X4
, http://www.100md.com
家庭经济状况
0能维持开支 1 不能维持开支
4.09
<0.05
1.66
X5
离婚或丧偶
0否 1 是
22.35
<0.01
3.46
X6
吸烟
, 百拇医药
0否 1 是
0.02
>0.05
X7
经常过量饮酒
0否 1 是
11.34
<0.01
2.48
X8
体育锻炼
1每天坚持 2 偶尔或无
3.68
, http://www.100md.com
>0.05
X9
业余爱好
0有 1 无
5.64
<0.05
1.63
X10
重视健康
0不能顾及或无所谓 1 是
5.62
<0.05
0.61
, 百拇医药
x11
慢性病
0无 1 有
5.03
<0.05
1.58
x12
睡眠不足
06小时及以上 1 6小时以下
6.06
<0.05
1.65
x13
, 百拇医药
工作时间延长
0无 1 有
7.95
<0.01
1.77
x14
接触毒物或噪声
0无 1 有
10.12
<0.01
1.93
多因素广义相对危险度模型分析结果如表3。为了探讨危险因素间联合作用的模式，对筛选出的因素进行两两间作用关系的分析，以λ与偏差度之间的动态变化图来判断实际数据趋向于何种模型，结果如表4。
, http://www.100md.com
表3 多因素分析的模型比较
因素
λ=0.0
λ=1.0
λ=1.2
P
R(x)
β
Sβ
β
Sβ
β
Sβ
, 百拇医药
标化β
X5
7.305
4.200
1.183
0.282
1.004
0.230
4.365
<0.01
3.044
X9
1.533
, http://www.100md.com
1.143
0.506
0.211
0.447
0.180
2.483
<0.05
2.082
X11
1.933
1.267
0.413
0.209
, 百拇医药
0.336
0.160
2.288
<0.05
1.941
X13
0.477
0.788
0.456
0.212
0.416
0.181
2.298
, 百拇医药
<0.05
2.041
X14
1.827
1.205
0.587
0.215
0.513
0.185
2.773
<0.01
2.173
偏差度
, 百拇医药
34.80
32.55
32.51
卡方值
31.35
28.44
28.56
表4 危险因素间作用关系的分析(R(x))
暴露状态
X5,X9
X5,X11
X5,X13
, http://www.100md.com
X5,X14
X9,X11
X9,X13
X9,X14
X11,X13
X11,X14
X13,X14
- -
1.00
1.00
1.00
1.00
1.00
, 百拇医药
1.00
1.00
1.00
1.00
1.00
- +
1.55
1.74
1.49
2.51
2.22
2.11
2.13
, http://www.100md.com
2.69
5.40
2.24
+ -
3.01
4.05
2.05
5.14
2.28
2.08
2.08
2.33
3.64
, 百拇医药
2.23
+ +
6.67
5.25
7.37
8.76
3.21
2.24
2.26
4.62
15.89
2.28
偏差度
, 百拇医药
0.0004
0.379
0.0024
0.0053
1.4500
0.0010
0.0020
0.0003
0.0737
0.0034
λ
2.2
0.25
, 百拇医药
2.4
1.4
1.4
2.6
2.6
0.8
0.6
4.2
讨 论
本研究定量分析了生理健康、社会经济、环境因素、行为生活方式对心理健康的影响，单因素的分析显示，家庭经济状况差、经常过量饮酒、离婚或丧偶、缺少业余爱好、患慢性病、睡眠不足、工作时间延长、接触毒物或噪声是影响中年知识分子心理健康的危险因素(相对危险度大于1.0)，年龄组和重视健康则是有利因素(相对危险度小于1.0)，年龄组的这一特点与全国情况相吻合[6]。
, 百拇医药
为了控制混杂因素的作用，应用多因素广义相对危险度模型筛选因素，结果显示以超相乘模型较为适宜，根据标化回归系数判定危险因素的的作用大小依次为离婚或丧偶、接触毒物或噪声、缺少业余爱好、工作时间延长、患慢性病，在固定其他因素作用不变的情况下，其相对危险度依次为3.044、2.173、2.802、2.041和1.941。
危险因素间的相互作用分析提示: 离婚或丧偶与缺少业余爱好、工作时间延长、接触毒物或噪声；缺少业余爱好与患慢性病、工作时间延长、接触毒物或噪声；工作时间延长与接触毒物或噪声之间两两联合作用的形式趋向于超相乘结构。离婚或丧偶与患慢性病；患慢性病与工作时间延长、接触毒物或噪声之间两两联合作用的形式趋向于中间结构模型。
目前国内流行病学调查常用logistic 回归模型分析危险因素的作用，但假设条件是各因素之间的相互作用呈相乘关系，实际工作中的数据不一定都符合此假设，若不符合此假设，模型便会失效，如有的情况下各因素之间的相互作用呈相加关系，应该使用相加模型[7]，或事先不清楚数据适合于何种模型。心理学领域的流行病学研究涉及的因素繁多，因素之间的关系错综复杂，一概采用logistic 回归模型来分析不一定妥当。Breslow提出的广义相对危险度模型函数，通过不同指数λ的幂变换，自变量与相对危险度之间的连接形式范围从次相加结构变化到超相乘结构，并可以从偏差度的图分析方法及拟合优度检验来较为客观地选择最佳模型，因而适用于心理学领域的流行病学研究。
, 百拇医药
*本研究为贵州省科委科学技术基金资助项目
参考文献
1 徐俊冕主编. 医学心理学. 第二版.上海:上海医科大学出版社，1996:183～190.
2 王庆石.统计指标导论.第一版.大连:东北财经大学出版社，1994:54～55.
3 Breslow，N.E.，and Storer，B.E.General relative risk functions for case-control studies. Am J Epidemiol， 1985，122(1)：149～162.
4 项永兵，高玉堂，金凡.相对危险度回归模型.中华流行病学杂志，1994，15(4):244～246.
5 赵宁，俞顺章，孙伟明. 相对危险度函数蔟模型在原发性肝癌危险因素研究中的应用. 中国卫生统计, 1995，12(4):13～17.
6 金华，吴文源，张明圆. 中国正常人SCL-90评定结果的初步分析.中国神经精神疾病杂志，1986，12(5)：260～262.
7 方积乾，陈和年.医学研究中logistic回归模型的正确应用(二).中国卫生统计, 1993，10(5):61～63.
(收稿日期：1999-04-19), 百拇医药