多元线性回归分析ppt课件.ppt

《多元线性回归分析ppt课件.ppt》由会员分享，可在线阅读，更多相关《多元线性回归分析ppt课件.ppt（68页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第十五章第十五章 多元线性回归分析多元线性回归分析(Multiple Linear Regression) 讲述内容讲述内容第一节第一节 多元线性回归多元线性回归第二节第二节 自变量选择方法自变量选择方法第三节第三节 多元线性回归的应用多元线性回归的应用 及其注意事项及其注意事项 人的人的体重体重与与身高、胸围身高、胸围血压值血压值与与年龄、性别、劳动强度、饮年龄、性别、劳动强度、饮食习惯、吸烟状况、家族史食习惯、吸烟状况、家族史糖尿病人的糖尿病人的血糖血糖与与胰岛素、糖化血红胰岛素、糖化血红蛋白、血清总胆固醇、甘油三脂蛋白、血清总胆固醇、甘油三脂射频治疗仪定向治疗脑肿瘤过程中，射频治疗仪定向

2、治疗脑肿瘤过程中，脑皮质的脑皮质的毁损半径毁损半径与辐射的与辐射的温度温度、与照、与照射的射的时间时间目的目的：作出以多个自变量估计应变量的多：作出以多个自变量估计应变量的多元线性回归方程。元线性回归方程。资料资料：应变量为定量指标：应变量为定量指标, ,一定是随机的；一定是随机的；自变量全部或大部分为定量指标，可以是自变量全部或大部分为定量指标，可以是随机变动的随机变动的, ,也可以是人为选定的也可以是人为选定的, ,若有少若有少量定性或等级指标需作转换。量定性或等级指标需作转换。用途用途：解释和预报。更精确：解释和预报。更精确第一节第一节 多元线性回归多元线性回归 变量：变量：应变量应变量

3、 1 个，自变量个，自变量m 个，共个，共 m+1 个。个。 样本含量：样本含量：n 数据格式见表数据格式见表15-1 回归模型一般形式：回归模型一般形式：eXXXYmm22110一、多元线性回归模型一、多元线性回归模型表表15-1 多元回归分析数据格式多元回归分析数据格式 条件条件(1)Y与mXXX,21之间具有线性关系。 (2)各例观测值), 2 , 1(niYi相互独立。 (3)残差 e 服从均数为 0、 方差为2的正态分布， 它等价于对任意一组自变量mXXX,21值，应变量 Y 具有相同方差，并且服从正态分布。 两自变量与应变量的散点图一般步骤一般步骤mmXbXbXbbY22110建立

4、回归方程建立回归方程(2)检验并评价回归方程检验并评价回归方程及各自变量的作用大小及各自变量的作用大小 二、多元线性回归方程的建立 例例15-1 27名糖尿病人的血清总胆固名糖尿病人的血清总胆固醇、甘油三脂、空腹胰岛素、糖化血红醇、甘油三脂、空腹胰岛素、糖化血红蛋白、空腹血糖的测量值列于表蛋白、空腹血糖的测量值列于表15-2中，中，试建立血糖与其它几项指标关系的多元试建立血糖与其它几项指标关系的多元线性回归方程。线性回归方程。表15-2 27名糖尿病人的血糖及有关变量的测量结果 用用 最小二乘法最小二乘法y 2201122()()kkQYYYbb Xb Xb X11 11221121 1222

5、221 122kkYkkYkkkkkkYl bl bl bll bl bl bll bl bl bl01122()kkbYb Xb Xb X求偏导数（一阶）()() , , j=1,2,k()(), 1,2,ijijiijjijjjYjjjXXlXXXXX XinXYlXXYYX Yjkn原 理最小二乘法统计软件包统计软件包最小最小43216382027060351501424094335X.X.X.X.Y三、假设检验及其评价 012:0mH ， 1:jH各 （j=1,2, ,m)不全为 0， 0.05 1. 方差分析法方差分析法：残回总SSSSSS残回残回（MSMSmnSSmSSF) 1/（

6、一）对回归方程) 1(mn ,mFF方差分析法可将回归方程中所方差分析法可将回归方程中所有自变量有自变量X1,X2,Xm,作为一个作为一个整体来检验他们与应变量整体来检验他们与应变量Y之之间是否具有线性关系，并对回间是否具有线性关系，并对回归方程的预测或解释能力做出归方程的预测或解释能力做出综合评价。综合评价。变异来源 自由度 SS MS F P 总变异 n-1 SS总 回 归 m SS回 SS回 /m MS回/MS残 残 差 n-m-1 SS残 SS残 /(n-m-1) 表15-4 例15-1的方差分析表 表15-3 多元线性回归方差分析表(0.05)(0.05)SPSS软件输出结果软件输出

7、结果 Model Sum of Squaresdf Mean SquareFSig.Regression ResidualTotal133.7107 88.8412 222.5519 4222633.42774.0382 8.28 0.000102 R，说明自变量mXXX,21能够解释Y变化的百分比，其值愈接近于 1，说明模型对数据的拟合程度愈好。本例 6008. 05519.2227107.1332R 表明血糖含量变异的 60%可由总胆固醇、 甘油三脂、胰岛素和糖化血红蛋白的变化来解释。 21SSSSRSSSS回残总总2. 决定系数决定系数R 2：3.复相关系数复相关系数 可用来度量应变量

8、与多个自变量间的线性相与多个自变量间的线性相关程度，亦即观察值 与估计值 之间的相关程度。Y YY 计算公式：2RR，本例7751060080.R 若m=1自变量，则有| r |R ， r 为简单相关系数。 （二）对各自变量的假设检验与评价 指明方程中的每一个自变量对Y的影响含义 回归方程中某一自变量jX的偏回归平方和表示模型中含有其它 m-1 个自变量的条件下该自变量对 Y的回归贡献， 相当于从回归方程中剔除jX后所引起的回归平方和的减少量， 或在 m-1个自变量的基础上新增加jX引起的回归平方和的增加量。 1. 偏回归平方和 在其它自变量存在于回归方在其它自变量存在于回归方程中的程中的下，

9、考察某一自下，考察某一自变量变量Xj对应变量对应变量Y的回归效的回归效应应 )(jXSS回表示偏回归平方和，其值愈大说明相应的自变量愈重要。 ()/1/(1)jjSSXFSSnm回残一般情况下，m-1 个自变量对 Y 的回归平方和由重新建立的新方程得到，而不是简单地把jjXb从有 m 个自变量的方程中剔出后算得。 12 1, 1nm平方和（变异） 回归方程中 包含的自变量 SS回 SS残 4321X,X,X,X 133.7107 88.8412 432X,X,X 133.0978 89.4540 431XX,X 121.7480 100.8038 421XX,X 113.6472 108.90

10、47 321XX,X 105.9168 116.6351 表表15-5 对例对例15-1数据作回归分析的部分中间结果数据作回归分析的部分中间结果 各自变量的偏回归平方和可以通过拟合包含不同各自变量的偏回归平方和可以通过拟合包含不同自变量的回归方程计算得到，表自变量的回归方程计算得到，表15-515-5给出了例给出了例15-115-1数数据分析的部分中间结果。据分析的部分中间结果。 11234234()(,)(,) 133.7107-133.0978=0.6129SSXSSXXXXSSX XX回回回 21234134()(,)(,) 133.7107-121.748011.9627SSXSSXX

11、XXSSX XX回回回 31234124()(,)(,) 133.7107-113.647220.0635SSXSSXXXXSSX XX回回回 41234123()(,)(,) 133.7107-105.916827.7939SSXSSXXXXSSX XX回回回 152. 0) 1427( /8412.881 /6129. 01F， 962. 2) 1427/(8412.881 /9627.112F 968. 4) 1427/(8412.881 /0635.203F ， 883. 6) 1427/(8412.881 /7939.274F 结结果果2. t 检验法 是一种与偏回归平方和检验完全等

12、价的一种方法。计算公式为jbjjSbt检验假设： H0：0j，jt服从自由度为1mn的 t 分布。如果12mn ,/jt|t |，则在（0.05）水平上拒绝 H0，接受 H1，说明jX与Y有线性回归关系。 jb为偏回归系数的估计值，jbS是jb的标准误。 390036560142401.t 721120420351502.t 229212140270603.t 623224330638204.t 结结 果果结结 论论0742222050.t,/.，074.2|34 tt， P值 均 小 于0.05， 说 明3b和4b有 统 计 学 意义 ， 而1b和2b则 没 有 统 计 学 意 义 。 对于

13、同一资料，不同自变量的对于同一资料，不同自变量的t值间可值间可以相互比较，以相互比较，t的绝对值越大，说明该的绝对值越大，说明该自变量对自变量对Y的回归所起的作用越大。的回归所起的作用越大。3标准化回归系数标准化回归系数 自变量取值的单位及其离散程度是不同的，因自变量取值的单位及其离散程度是不同的，因此量纲不同的各回归系数之间不能直接比较大小，此量纲不同的各回归系数之间不能直接比较大小，可对变量进行标准化变换，然后拟合回归方程，这可对变量进行标准化变换，然后拟合回归方程，这样获得的回归系数称为标准化回归系数样获得的回归系数称为标准化回归系数。 变量标准化是变量标准化是将原始数据减去相应变量的均

14、数，将原始数据减去相应变量的均数，然后再除以该变量的标准差。然后再除以该变量的标准差。()jjjjXXXS 计算得到的回归方程称作标准化回归方程，相应的回归系数即为标准化回归系数。 YjjYYjjjjSSbllbb注意：注意： 一般回归系数一般回归系数有单位，用来解释各自有单位，用来解释各自变量对应变量的影响，表示在其它自变量保变量对应变量的影响，表示在其它自变量保持不变时，持不变时， 增加或减少一个单位时增加或减少一个单位时Y Y的平的平均变化量均变化量 。不能用各。不能用各 来比较各来比较各 对对 的影响大小。的影响大小。 标准化回归系数标准化回归系数无单位，用来比较各无单位，用来比较各自

15、变量对应变量的影响大小，自变量对应变量的影响大小， 越大，越大， 对对 的的 影响越大。影响越大。jXjbjXYjb jXY11.5934S ，22.5748S ，33.6706S ，41.8234S ，2.9257YS 0776. 09257. 25934. 11424. 01b 309309257257482351502.b 339509257267063270603.b 397709257282341638204.b 结结论论结果显示， 对血糖影响大小的顺序依次为糖化血红蛋白)(4X、 胰岛素)(3X、 甘油三脂)(2X和总胆固醇)(1X。 第四节 自变量筛选目的目的：使得预报和（或）解

16、释效果好 多元线性回归方程中所包含的自变量是根据专多元线性回归方程中所包含的自变量是根据专业知识和经验事先选择好的，但在许多回归分析的、业知识和经验事先选择好的，但在许多回归分析的、应用中，由于没有清晰的理论依据，回归模型所包应用中，由于没有清晰的理论依据，回归模型所包含的自变量难以预先确定，如果将一些不重要的自含的自变量难以预先确定，如果将一些不重要的自变量也引入方程，会降低模型的精度，因此选择有变量也引入方程，会降低模型的精度，因此选择有意义的自变量是回归分析的第一步。意义的自变量是回归分析的第一步。 选择自变量的基本思路是：尽可能将回归效果选择自变量的基本思路是：尽可能将回归效果显著的自

17、变量选入回归方程中，将作用不显著的特显著的自变量选入回归方程中，将作用不显著的特别是与自变量有密切线性关系的自变量排除在外。别是与自变量有密切线性关系的自变量排除在外。一、全局择优法目的：根据一些准则建立根据一些准则建立 “最优最优”回归模型回归模型意义：对自变量各种不同的组合所建立 的回归方程进行比较 ， 择优选择方法：决定系数决定系数R2的缺点：的缺点： 当回归方程中包含有很多自变量，即当回归方程中包含有很多自变量，即使其中有一些自变量对解释应变量变使其中有一些自变量对解释应变量变异的贡献极小，随着回归方程的自变异的贡献极小，随着回归方程的自变量的增加，量的增加，R2 值表现为只增不减。值

18、表现为只增不减。1校正决定系数2cR选择法，其计算公式为 总残MSMSpnnRRc111)1 (122 n 为样本含量，2R为包含)(mpp个自变量的回归方程的决定系数。2cR的变化规律是：当2R相同时，自变量个数越多2cR越小。所谓“最优”回归方程是指2cR最大者。 2. pC选择法 )1(2)()(pnMSSSCmpp残残 pSS)(残是由)(mpp个自变量作回归 的误差平方和，mMS)(残是从全部 m 个自 变量的回归模型中得到的残差均方。 例15-2 用全局择优法对例15-1数据的自变量进行选择。 方 程 中 的 自 变 量 2cR pC 方 程 中 的 自 变 量 2cR pC X2

19、,X3,X4 0.546 3.15 X2,X3 0.408 9.14 X1,X2,X3,X4 0.528 5.00 X1,X3 0.375 10.78 X1,X3,X4 0.488 5.96 X4 0.347 11.63 X1,X2,X4 0.447 7.97 X1 0.284 14.92 X1,X4 0.441 7.42 X1,X2 0.275 15.89 X2,X4 0.440 7.51 X3 0.231 17.77 X3,X4 0.435 7.72 X2 0.179 20.53 X1,X2,X3 0.408 9.88 m=4, 故回归方程拟合数为4212115m 。 最优组合均为X2，X

20、3，X4，即由甘油三脂、胰岛素和糖化血红蛋白与血糖建立的回归方程最优。 如果自变量个数为4，则所有的回归模型有241 15个；当自变量数个数为10时，所有可能的回归为 2101 1023个；当自变量数个数为50时，所有可能的回归为25011015个。二、逐步选择法1. 1.前进法，回归方程中的自变量从无到有、从少到多逐个引入回归方程。此法已基本淘汰。 2. 后退法，先将全部自变量选入方程，然后逐步剔除无统计学意义的自变量。 剔除自变量的方法是在方程中选一个偏回归平方和最小的变量，作F检验决定它是否剔除，若无统计学意义则将其剔除，然后对剩余的自变量建立新的回归方程。重复这一过程，直至方程中所有的

21、自变量都不能剔除为止。理论上最好，建议使用采用此法。 3.逐步回归法，逐步回归法是在前述两种方法的基础上，进行双向筛选的一种方法。该方法本质上是前进法。 1.从不包含任何自变量的模型开始（截距）从不包含任何自变量的模型开始（截距）2.添加具有最大统计学意义的变量进入模添加具有最大统计学意义的变量进入模型，如型，如p值最小，小于预先给定的进入标值最小，小于预先给定的进入标准的变量。准的变量。3.重复第二步，直到没有重复第二步，直到没有p值小于预先给定值小于预先给定的进入标准的变量。的进入标准的变量。1.从整个模型开始（包含研究所涉及所有自变量及从整个模型开始（包含研究所涉及所有自变量及其交互项其

22、交互项2.剔除对参数进行剔除对参数进行t检验或方差分析的结果中最无统检验或方差分析的结果中最无统计学意义的变量。如计学意义的变量。如p值最大，大于预先给定的剔值最大，大于预先给定的剔除标准的变量。除标准的变量。3.重复第二步，直到没有重复第二步，直到没有p值大于预先给定的剔除值大于预先给定的剔除标准的变量。标准的变量。1.从不包含任何自变量的模型开始。从不包含任何自变量的模型开始。2.添加具有最大统计学意义的变量进入模型。添加具有最大统计学意义的变量进入模型。3.剔除对参数进行剔除对参数进行t检验或方差分析的结果检验或方差分析的结果中最无统计学意义的变量。中最无统计学意义的变量。3.重复第二步

23、和第三步，直到模型外没有可重复第二步和第三步，直到模型外没有可进入的变量或模型内没有可剔除的变量为止。进入的变量或模型内没有可剔除的变量为止。例 15-3 试用逐步回归方法分析例 15-1 数据（100.入，15. 0出） 。 步 骤 (l) 引 入 变 量 剔 除 变 量 变 量 个 数p 2R ( )()lSSXj回 ( ) lSS残 F值 P值 1 X4 1 0.372 82.714 139.837 14.788 0.0007 2 X1 2 0.484 25.076 114.762 5.244 0.0311 3 X3 3 0.547 13.958 100.804 3.185 0.0875

24、 4 X2 4 0.601 11.963 88.841 2.962 0.0993 5 X1 3 0.598 0.613 88.841 0.152 0.7006 表15-7 逐步回归过程 变异来源 自由度 SS MS F P 总变异 26 222.5519 回 归 3 133.098 44.366 11.41 0.0001 残 差 23 89.454 3.889 表15-8 例15-3方差分析表 “最优”回归方程为4326632. 02871. 04023. 04996. 6XXXY 结果表明：血糖的变化与甘油三脂、胰岛素和糖化血红蛋白有线性回归关系，其中与胰岛素负相关。由标准化回归系数看出，糖

25、化血红蛋白对空腹血糖的影响最大。 表15-9 例15-3的回归系数的估计及检验结果 变 量 回归 系数b 标准误 bS 标准回 归系数b t 值 P值 常数项 6.4996 2.3962 0 2.713 0.0124 X2 0.4023 0.1540 0.3541 2.612 0.0156 X3 0.2870 0.1117 0.3601 2.570 0.0171 X4 0.6632 0.2303 0.4133 2.880 0.0084 一、多元线性回归的应用影响因素分析，控制混杂因素预测：由自变量值推出应变量Y的值控制：指定应变量Y的值查看自变量的改变量 （一）变量的数量化（一）变量的数量化

26、（二）样本含量（二）样本含量 （三）统计（三）统计“最优最优”与专业的与专业的“最优最优” （四）多重共线性（四）多重共线性 （五）交互作用（五）交互作用 （六）残差图（六）残差图1. 影响因素分析 例如影响例如影响高血压的因素高血压的因素可能有可能有年龄年龄、饮饮食习惯食习惯、吸烟状况、吸烟状况、工作紧张度工作紧张度和和家族史家族史等，在影响高血压的众多可疑因素中，需等，在影响高血压的众多可疑因素中，需要研究哪些因素有影响，哪些因素影响较要研究哪些因素有影响，哪些因素影响较大。大。 在临床试验中，则可能由于种种原因难以在临床试验中，则可能由于种种原因难以保证各组的指标基线相同，如在保证各组的

27、指标基线相同，如在年龄年龄、病病情情等指标不一致出现混杂的情况下，如何等指标不一致出现混杂的情况下，如何对不同的治疗方法进行比较等。对不同的治疗方法进行比较等。 这些问题都可以利用这些问题都可以利用回归分析回归分析来处理。来处理。控制混杂因素控制混杂因素( (confounding factor) )的一个的一个简单办法就是将其引入回归方程中，简单办法就是将其引入回归方程中，与其与其他主要变量一起进行分析他主要变量一起进行分析 2. 估计与预测 如由儿童的如由儿童的心脏横径心脏横径、心脏纵径心脏纵径和和心脏宽心脏宽径径估计估计心脏的表面积心脏的表面积；由胎儿的；由胎儿的孕龄孕龄、头头颈颈、胸径

28、胸径和和腹径腹径预测预测出生儿出生儿体重体重等。等。 3. 统计控制 逆估计。 例如采用射频治疗仪治疗脑肿瘤，脑例如采用射频治疗仪治疗脑肿瘤，脑皮质的皮质的毁损半径毁损半径与与射频温度射频温度及及照射时间照射时间有有线性回归关系，建立回归方程后可以按预线性回归关系，建立回归方程后可以按预先给定的脑皮质毁损半径，确定最佳控制先给定的脑皮质毁损半径，确定最佳控制射频温度和照射时间。射频温度和照射时间。 二、多元线性回归应用的注意事项二、多元线性回归应用的注意事项（1）定量，非线性 线性 （2）定性，转为定量： （0，1）变量，哑变量，伪变量，指示变量。 1指标的数量化指标的数量化 分2类，可用一个

29、（0，1）变量。如性别0 1 X男女 分k类，k-1个（0，1）变量，如血型。 血型 X1 X2 X3 O 0 0 0 A 1 0 0 B 0 1 0 AB 0 0 1 编号 X1 X2 X3 Y 1 1 0 0 2 0 0 0 3 0 1 0 n 0 0 1 数据格式回归方程 建立回归方程建立回归方程 0112233Ybb Xb Xb Xb1 ：相当相当A 型相对于型相对于O 型的差别型的差别b2 ：相当：相当B 型相对于型相对于O 型的差别型的差别b3 ：相当：相当AB 型相对于型相对于O 型的差别型的差别 （3）等级 定量。 一般是将等级从弱到强转换为 （或 ）如文化程度分为小学、中学、

30、大学、大学以上四个等级。Y为经济收入。1 1 2 3 X小学中学大学 4 大学以上 011Ybb X解释：b(b1)反映X（X1） 增加1个单位， 增加b个单位（如：500元）。 表示中学文化较小学收入多500， 大学较中学多500，余类推。, 3 , 2 , 1X, 2 , 1 , 0XYb1,b2,b3分别反映中学、大学、大学以上相对于小分别反映中学、大学、大学以上相对于小学文化程度者经济收入差别的大小学文化程度者经济收入差别的大小也可将K个等级转换为K-1个（0，1）变量 2样本含量: n =(510)m。 3关于逐步回归: 对逐步回归得到的结果不要盲目的信任，所谓的“最优”回归方程并不

31、一定是最好的，没有选入方程的变量也未必没有统计学意义。例如，例15-3中若将选入标准和剔除标准定为 和 选入的变量是 ， 而不是 ， 结果发生了改变。 不同回归方程适应于不同用途，依专业知识定。05. 0入10. 0出41, XX432,XXX不同准则、方法得出的“最优”方程不同不同的引入、剔除标准获得的“最优”方程不同；方程还受数据的正确性、共线性影响研究者应结合问题本身和专业知识及经验来决定 自变量间存在着相关关系，使一个或几个自变量可以由另外的自变量线性表示时，称为该变量与另外的自变量间存在有共线性(collinearity)。xxxxxxxxxxxxxX1X2Y 4. 多重共线性 即指

32、一些自变量之间存在较强的线性关系。如高血压与年龄、吸烟年限、饮白酒年限等，这些自变量通常是高度相关的，有可能使通过最小二乘法建立回归方程失效，引起下列一些不良后果： （1）参数估计值的标准误变得很大，从而t值变得很小。 （2）回归方程不稳定，增加或减少某几个观察值，估计值可能会发生很大的变化。 （3）t检验不准确，误将应保留在模型中的重要变量舍弃。 （4）估计值的正负符号与客观实际不一致。 消除多重共线性：剔除某个造成共线性的自变量，重建回归方程；合并自变量；采用逐步回归方法。 筛选自变量用主成分回归岭回归回归系数的符号与专业知识不符变量的重要性与专业不符R2高，但各自变量对应的回归系数均不显

33、著。方差膨胀因子（Variance Inflation Factors ，VIF) 10 当某一自变量对应变量的作用大小与另一个自变量的取值有关时，则表示两个变量有交互作用（interaction）。 检验两变量间有无交互作用，普遍的做法是在方程中加入它们的乘积项再做检验。如考察X1、X2间的交互作用，可在模型中加入X1X2项。如对表15-2数据的分析， 已经选出甘油三酯)(2X、 胰岛素)(3X和糖化血红蛋白)(4X三个变量，在方程中加入43XX项，如果这一乘积项显著， 则说明胰岛素和糖化血红蛋白之间有交互作用。 为此需要定义一个新的变量43XXZ ，按方程 ZbXbXbXbbYZ44332

34、20 估计参数。若经检验， 拒绝0H：0z，便可认为除3X和4X的主效应外还有交互作用。本例结果为 ZXXXY1785. 05097. 12267. 13690. 07898. 0432 经检验Z的作用显著(P 0.01)，说明糖尿病患者体内胰岛素对血糖的影响依赖于糖化血红蛋白的含量。 （六）（六）、残差图（、残差图（residual plot） x or y 0ResidualsHomoscedasticity: Residuals appear completely random. No indication of model inadequacy.0ResidualsCurved pat

35、tern in residuals resulting from underlying nonlinear relationship.0ResidualsResiduals exhibit a linear trend with time.Time0ResidualsHeteroscedasticity: Variance of residuals changes when x changes.x or y x or y 从菜单选择 Regression LinearMethod中的选项： Enter 所选自变量强行全部选入方程 Stepwise 逐步回归法重要的结果R、R Square、 A

36、djust R Square反映模型对应变量变异的解释程度。Adjust R Square尤其用于自变量个数不等的模型之间的比较。ANOVA方差分析表中的Sig.反映模型是否有统计学意义Coefficients表中B为各自变量系数（回归系数）。Beta为标准化回归系数，其绝对值用于说明自变量的重要性。t、Sig.检验各系数是否有统计学意义。逐步回归的重要选项Method要选为StepwiseOptions中要设定合适的选入选出标准（注：Removal中设定的p值必须大于Entry中的设定。） 参考书目参考书目1.孙振球主编孙振球主编 . 医学统计学医学统计学 . 第第 2 版版 . 北京：人民

37、北京：人民卫生出版社，卫生出版社， 2005 2. 方积乾主编方积乾主编 , 孙振球副主编孙振球副主编 . 卫生统计学卫生统计学 . 第第 5 版版 . 北京：人民卫生出版社，北京：人民卫生出版社， 2003 3. 徐勇勇主编徐勇勇主编 , 孙振球副主编孙振球副主编 . 医学统计学医学统计学 . 第第 2 版版 . 北京：高等教育出版社，北京：高等教育出版社， 2004 4. 方积乾主编方积乾主编 . 医学统计学与电脑实验医学统计学与电脑实验 . 第第 2 版版 . 上海：上海科学技术出版社，上海：上海科学技术出版社， 2001 5. 赵耐青主编赵耐青主编 . 医学统计学医学统计学 . 高等教育出版社，高等教育出版社， 2004