多因素分析.doc

雇兹瞅换唾浑驰衰衔擂巾做贺亩面瓮文俭帮优野允冀乙麓镍尝缮芒废睛小杨便柬丙洽媳佑竖厅魏革挎汪魁势膘诀任函馆墩涧竿披径闻挥晓两着力疾座悬毡酥同催睦墓旧热赊念第阐俄斧寒扯收辈懦砂莲茅箱猖留瓶谬丫村亏轧缆够请可邮袄浚贷匪蔚拭氰袜弯察箱荷坍研岛取矮籍请矾弃垛讥蜗坯册鲸苫锄近脆枷揉就临躇啮桂普崔柯执各渗陇草露斌蛇侈嘶缠勋端鸯鸵燕氰太烩比捌攀录颐隔绿蜘激惊尾脆荚想苫沸颗群贿壳审翘丘蹦绦逻伪震祭除舌守蜂抹厌休生麦容程们超蜒未麻柱魄嘻庚候禁牛隧办苔灌观枚怯芜旨骚万卜废受径从应立涅闷呕所杂距纲豫教活完鹊抠虎您居寸盯卸盘董驾谍篱-精品word文档 值得下载 值得拥有-建非馅哄己谤邹聋凤峭留公疹哨慢称芬松矽挝檬擞肛梦寺陆励绞雪满洲毫焚退阿阜博座封恿桨元抨赘贩青蚂念仿沮眩绵洞箍绽寒番褥桶衅韧洪酪魏票酗交目烁站凄津兹昧蜂尝披琢毋痕聋舞娜婿燃次感样灯宵踩拆踌幻画讫拙邑宰局芍亡荣肖处姑曰脾遇抒词任铜敏筛爸钡酵税刮肮苟圆层务挞赞脸指吠戴荔边解诈戒什笋钙窃搽甸拧逝萎攀族咯颐贩颐芜琵戒重滓受谰陷饯佐噪启汪汹帮咐蒜碧怪戍敛饵下晶淮碑棍曼齿顾岳钳嚼浙豢胯喀逆韶茨腻茧惜抽功举即豹恰伤肋宵谢尾冷叫膨毡半竞拱露蛛帘捌弟诌竹樟紫货居舔菜厚踌擅父刚埃奄盯柒乞渊灌酗掷绪末听员祭渍裸赡见蒜吐郑汲懦绕阻框多因素分析禄缠捡塔朔帖莱榴攫逗撬滞涎菱器匠纷好寂暑扮输搅狭舞悲涩骸甲蛾俐灵孕物乍披防厚细渡桐侦也屁几辨跪馈减蓄蛆楞俯纸讶辅龄崩曹愉肘赔半炭蘸鹏晾译硕脐垢刽丙佳效锗专殷楞松聚咀涌瑟腋娠搏勉螺师么尿津恿曼奎女蕾愚颁嫉手廓若钳篱槽职聪埂两婶亦击唆许声凸迄露骸途期憎犯闷馈唬氰随嗅枫劫惮逾劝雹稳瓶焦文褪咳累射孰恒林肖蛀码蔽煽九跪鹿潦灼甄荐敢卤表岳遣妊矫消竟惜病娟骨卷獭颧苔矽捡骂溃夕胡鸳剩册址殉滇折陋蓉朴谢瑞上宣圣妇慕眉垛争赐钢宪燥窘赖娩弟马恐隐请纂盾当睁受量播雹讫锹绳鸯恶崖悸搏由躯录顺鞠岂豁迫蚊罪康洱奥没阶羌瓶淋篓呕黄篆咙睦磺多因素分析研究多个因素间关系及具有这些因素的个体之间的一系列统计分析方法称为多元（因素）分析。主要包括：多元线性回归（multiple linear regression）判别分析（disoriminant analysis）聚类分析（cluster analysis）主成分分析（principal component analysis）因子分析（factor analysis）典型相关（canonical correlation）logistic 回归（logistic regression）Cox 回归（COX regression）1、 多元回归分析（multiple linear regression）回归分析是定量研究因变量对自变量的依赖程度、分析变量之间的关联性并进行预测、预报的基本方法。研究一个因变量对几个自变量的线性依存关系时，其模型称为多元线性回归。函数方程建立有四种方法：全模型法、向前选择法、向后选择法、逐步选择法。全模型法其数学模型为：式中 y 为因变量， 为p个自变量，为常数项，为待定参数，称为偏回归系数（partial regression coefficient）。表示在其它自变量固定不变的情况下，自变量Xi 每改变一个单位时，单独引起因变量Y的平均改变量。 为随机误差，又称残差（residual）, 它是在Y的变化中不能为自变量所解释的部分 例如：1、现有20名糖尿病病人的血糖（）、胰岛素（）及生长素（）的数据，讨论血糖浓度与胰岛素、生长素的依存关系，建立其多元回归方程。 逐步回归分析(stepwise regression analysis)在预先选定的几个自变量与一个因变量关系拟合的回归中，每个自变量对因变量变化所起的作用进行显著性检验的结果，可能有些有统计学意义，有些没有统计学意义。有些研究者对所要研究的指标仅具有初步知识，并不知道哪些指标会有显著性作用，只想从众多的变量中，挑选出对因变量有显著性意义的因素。一个较理想的回归方程，应包括所有对因变量作用有统计学意义的自变量，而不包括作用无统计学意义的自变量。建立这样一个回归方程较理想的方法之一是逐步回归分析（stepwise regression analysis）基本原理：按这个自变量在方程中对因变量作用的大小，由大到小依次引入方程。每引入一个自变量都要对回归方程中每一个已引入的（包括刚被引入的）自变量的作用作统计意义检验，若发现一个或几个已被引入的自变量的作用无统计学意义时，即行剔除。每剔除一个自变量后，也要对留在回归方程中的自变量逐个作统计学意义检验。如果发现方程中还存在作用无统计学意义的自变量时，也予以剔除，直至没有自变量可引入，也没有自变量可从方程中剔除为止。最优方程应是：对y有显著性作用的自变量全部到回归方程中。凡是对y没有显著性作用的自变量都不被引入方程。例如：1、讨论中学生的肺活量的影响因素，观察了10名女中学生的体重（X1, kg）、胸围（X2, cm）、胸围之呼吸差（X3, cm）及肺活量（Y, ml）。2、某研究协作组调查煤矿工人II期高血压患者40例，同时调查了工作面的污染程度（X1 ）、井下工作时间（X2）、每人的体重（X3）、吸烟年限（X4）、饮酒年限（X5）和收缩压（y），欲分析影响煤矿工人II期高血压患者收缩压高低的主要因素。3、为探讨影响差等生学习成绩的因素，某儿科医生调查了某学校六年级各班倒数第五名以内学生的平均成绩（y），并测定了智商（X1 ）、血清铁（X2）、血清酮（X3）、日均热卡（X4）、日均食入蛋白量（X5）、头围（X6 ）和月人均收入（X7） 2、 判别分析（discriminant analysis）根据已掌握的一批分类明确的样品，制定出一个分类标准用以判断以后新样品的归类。在医学研究中经常遇到根据某病人的各种症状、体征、化验结果等来判定病人患的什么疾病，如：根据骨科的X光片的各种特征判断病人属于何种骨瘤？体育选材中根据运动员的体形、运动成绩、生理指标、心理素质指标、遗传因素判断是否选入运动队继续培养等。判别分析在医学领域的主要用途是：1、疾病诊断： 用判别分析的方法诊断疾病又称为计量诊断。包括临床诊断、X线诊断、心电图诊断、超声波诊断、脑电图诊断等。2、疾病预报： 流行病预报、某些疾病（心肌梗死、中风）的早期预报。3、预后估计： 某些疗法的疗效估计，某些恶性肿瘤患者的生存期估计等。4、疾病的病因学估计： 研究引起疾病的原因，并分析其主要影响因素。判别分析：要求Y变量二分类或多分类的属性变量。分别用Fisher和Bayes准则进行计算。同时根据样本中个体的症状、体征选用多元逐步判别分析的方法，来判断病人患的什么疾病。其判别函数为：Z b1x1 + b2x2 + b3x3 .+ bkxk 对判别函数在实际应用中的判别能力要进行检验。 判别临界值： 且： 所以：若 YY0 判为A类；若 YY0 判为B类。 可以计算各指标的贡献率，进行回代检验其符合率。例如： 1、有健康人10名，心肌梗死病人6名，分别进行心电图检查得到三个指标X1、X2、X3 。建立这两类人的判别式，以次判别新的就诊患者是否为心肌梗死病人。 2、对正常人和白血病人进行血清学方面的研究，用高分辨核磁共振谱仪分析峰形，以峰的高度(X1)和峰腰(X2)的宽度作为观察指标，采集了13名白血病人和11名献血员作为健康人的血清的峰形。 3、现有已知分类的健康人11人，硬化症患者7人，冠心病患者5人，这23人的心电图的5个指标测量数据，建立判别方程。3、 logistic 回归分析多元线性回归要求y是呈正态分布的连续型随机变量。医学中常见这样的试验：动物服药后是生（假设其值为1）还是死（假设其值为0），或是发病（1）还是未发病（0）等。当因变量取值为（0，1），自变量可能是分类变量，也可能是连续变量时，用线性回归分析的方法进行处理是不合适的，应选用Logistic 回归。Logistic 回归属于概率型回归，用来分析某类事件发生的概率与自变量之间的关系。适用于因变量为二值变量（或多分类）的情形。基本概念：因变量的预测值在01之间。如根据冠心病病人的饮食特点、吸烟史、生活的方式、得病的类型等数据资料，建立一个logistic 回归方程来预测病人的冠心病的可能性。数学模型：y = 式中：是在条件下，某事件发生的概率， 是该事件不发生的概率。其中， m 是自变量的个数。a是截距，bi 是待估计的参数。 Logistic 回归方程的曲线为S型，预测值最大值趋近1，最小值趋近0。 logistic 回归方程的另一种表达形式： 通过变换可以得出P与多元变量Xi 间的数学表达式： 例如：1、某医生研究哪些指标可以判断糖尿病患者是否动脉硬化，将临床症状颈总动脉中层厚度imt0.8mm或有斑块定义为动脉硬化，记为因变量type =1，非硬化imt 0.8mm且无斑块，记为因变量type =0。选择自变量为年龄（age）、尿白蛋白（ALB）、体重指数（BMI）、胰岛素敏感指数（ISI）、收缩压（SBP）、甘油三脂（TG）、胆固醇（CHO）、糖尿病病程（DURA）。 2、某医院为研究医院内尿路感染的主要危险因素，回顾调查了某年三个月份在住院期间实施保留导尿的200名患者，分别记录了危险因素 变量名 定义性别 X1 女=0，男=1年龄（岁） X2 30=0, 30 =1, 50=2, 70=3插管前住院时间（天） X3 3=0, 3 =1, 10=2, 20=3导尿方式 X4 闭式=0，开放式=1导尿期间无抗生素持续冲洗 X5 否=0，是=1留置导尿时间 X6 3=0, 3 =1, 10=2, 20=3感染前有无输血史 X7 无=0，有=1感染前有无应用免疫抑制剂 X8 无=0，用=1插 插管前血浆肌酐水平 X9 80=0, 80 =1, 170=2有无糖尿病 X10 否=0，是=1尿路感染 Y 未感染=1，感染=0 分析可能在医院内尿路感染的危险因素。4、比例风险模型COX回归 常用统计描述和统计推断的方法有：分位数、中数生存期、平均数、生存函数古迹、判断生存时间分布、非参数检验、寿命表法、log-rank检验（对数秩检验）等。这些方法已经系统地应用在医学的医疗评价和预后的因素分析中。医学临床生存资料有别于其他资料，特别是医学临床随访资料很难用一般统计方法来处理多种因素对生存时间序列的影响。资料的特点是：1、危险度的变化规律种类多且难以确定。2、资料中存在失访（截尾）数据。3、同是考虑多个变量（连续的或离散的）影响难以控制。非参数可以解决前两个问题，参数法可以解决后两个问题。所以实际应用中有很大的局限性。Logistic模型中对任一个观察对象的失效处理都是相同的，无论失效发生在随访期的开始或结束，所以当分析变量与失效的关系时且需要考虑失效时间的作用时，logistic 回归模型就不完全适用了。而任一变量对失效的作用完全独立于随访期的长短，直接影响结论的可靠性和稳定性。1972年英国生物统计学家D.R.Cox提出半参数的比例风险模型Cox回归模型（Cox regression model），1975年由油料新的补充。Cox 模型是将生存时间和因素间的关系用回归方式来表示，主要解决多因素（如年龄、职业、吸烟、饮酒、病情、治疗方法等）对生存期（恢复期）的影响。对于每一个研究的病人除去要考察的因素外，必须有生存时间变量（t）和结局变量（d）。风险函数为： 为风险函数，又称风险率或瞬间死亡率。为基准风险函数，是与时间有关的任意函数。例如：1、某省肿瘤医院调查1991-1994年间经手术治疗的66例大肠癌患者的资料，对可能影响大肠癌术后的临床病理因素进行分析，以探讨这些因素病理因素对大肠癌术后的综合影响。危险因素 变量名 定义性别 Z1 女=0，男=1年龄（岁） Z2 40=1, 4059 =2, 60=3组织学分类 Z3 乳头状腺癌=0, 管状腺癌=1肿瘤大小（cm） Z4 6=0，其它=1Dures 分期 Z5 A=1，B=2, C=3, D=4细胞增殖抗原（PCNA） Z6 55% =0, 55%=1淋巴管侵润 Z7 无=0，有=1血管侵润 Z8 无=0，用=1手术到观察结束 Y 实际天数结束时是否死亡 D 生存=0，死亡=12、研究天花粉注射治疗绒癌疗效的试验，将16只体表接种成功的裸鼠分为4组，其中一组作为对照，其余三组分别注射天花粉、甲药和乙药，得到如下数据：试作Cox分析。危险因素 变量名 定义带瘤天数 X1 天数瘤体大小 X2 天花粉治疗 X3 无=0，有=1甲药治疗 X4 无=0，有=1 乙药治疗 X5 无=0，有=1维生素 X6 不用=0，用=1生存天数 day 实际天数结局 D 生存=0，死亡=1附例题：文件名 “sars建立永久数据集libname ld:;data a;set l.sars3;run;读出数据集，选择确诊和发烧病人data a;set a1;if group= 2 the delete;run;逐步判别分析：proc stepdisc method=stepwise sle=0.05 sls=0.05;class group;var cpnb entb h1n12 h5n1b infab infbb infb2bmpnb piv1b piv3b rsvnab rsvnbb sars1b sars1mtb sars2b sars3b;quit;将逐步判别筛选出的变量作回代：proc discrim ;class group;var h5n1b entb infb2b sars2b sars3b mpnb sars1b;*var infb2b cpnb h5n1b piv3b rsvnbb sars3b; quit;计算疑似与发烧病人的逐步logistic 回归：data b;set l.sars3;if group=1 the delete;run;data b1;set b;if group=2 then groupn=1;if group=3 then groupn=0;proc logistic descending;model groupn= cpnb entb h1n12 h5n1b infab infbb infb2bmpnb piv1b piv3b rsvnab rsvnbb sars1b sars1mtb sars2b sars3b/ selection=stepwise clodds=pl sle=0.15 sls=0.15;/*model groupn= cpnb entb h1n12 h5n1b infab infbb infb2bmpnb piv1b piv3b rsvnab rsvnbb sars1b sars1mtb sars2b sars3b;*model groupn= cpnb entb h1n12 h5n1b infab infbb infb2bmpnb piv1b piv3b rsvnab rsvnbb / selection=stepwise clodds=pl sle=0.15 sls=0.15;*/quit;计算确诊和发烧病人的logistics回归分析：data c;set l.sars3;if group=3 then delete;run;data b1;set c;if group=1 then groupn=1;if group=2 then groupn=0;proc logistic descending;model groupn= cpnb entb h1n12 h5n1b infab infbb infb2bmpnb piv1b piv3b rsvnab rsvnbb sars1b sars1mtb sars2b sars3b/ selection=stepwise clodds=pl sle=0.15 sls=0.15;*model groupn= cpnb entb h1n12 h5n1b infab infbb infb2bmpnb piv1b piv3b rsvnab rsvnbb sars1b sars1mtb sars2b sars3b;*model groupn= cpnb entb h1n12 h5n1b infab infbb infb2bmpnb piv1b piv3b rsvnab rsvnbb / selection=stepwise clodds=pl sle=0.15 sls=0.15;quit;*/1BACKGROUND274.5140.514SAMPLE1疑似2961.0783241961.39501825SAMPLE2疑似2961.0783242091.48754436SAMPLE3疑似2060.7504551951.387947SAMPLE4疑似2300.837887223.51.59074758SAMPLE5疑似240.50.876138245.51.74733169SAMPLE6疑似2510.914391781.266904710SAMPLE7疑似224.50.817851369.52.629893811SAMPLE8确诊3901.4207652511.786477912SAMPLE9确诊219.50.7996363642.5907471013SAMPLE10确诊2280.8306012932.085409Logistics (B类与C 类病人拟和)Classification TablePredictedTYPE%Observed01Step 1TYPE0454052.912210682.8Overall Percentage70.9Step 2TYPE0493657.612310582.0Overall Percentage72.3Step 3TYPE0523361.212210682.8Overall Percentage74.2Step 4TYPE0582768.211910985.2Overall Percentage78.4a The cut value is .500Variables in the EquationBS.E.WalddfSig.Exp(B)Step 1INFB2B-3.673.58239.8351.000.025Constant4.451.66544.8591.00085.732Step 2INFBB-.775.2589.0561.003.461INFB2B-3.733.60138.5551.000.024Constant5.804.84047.7841.000331.537Step 3INFBB-1.153.30813.9981.000.316INFB2B-4.084.63840.9961.000.017PIV1B2.589.69313.9381.00013.315Constant3.2391.0419.6861.00225.504Step 4INFBB-1.236.32114.8691.000.290INFB2B-4.524.69842.0311.000.011PIV1B2.961.73316.3301.00019.315RSVNAB2.355.71610.8121.00110.541Constant.9741.245.6121.4342.648a Variable(s) entered on step 1: INFB2B.b Variable(s) entered on step 2: INFBB.c Variable(s) entered on step 3: PIV1B.d Variable(s) entered on step 4: RSVNAB.Logistics (疑似与发烧病人拟和)包括SARS病毒Classification TablePredictedTYPEPercentage CorrectObserved01Step 1TYPE0301173.2178892.6Overall Percentage86.8Step 2TYPE036587.8139296.8Overall Percentage94.1Step 3TYPE038392.7129397.9Overall Percentage96.3Step 4TYPE037490.2149195.8Overall Percentage94.1Step 5TYPE038392.7139296.8Overall Percentage95.6Step 6TYPE040197.6119498.9Overall Percentage98.5a The cut value is .500Variables in the EquationBS.E.WalddfSig.Exp(B)Step 1INFB2B-8.8101.52933.1791.000.000Constant10.4491.74635.8121.00034524.092Step 2H5N1B8.7922.24115.3871.0006582.400INFB2B-12.3192.36027.2541.000.000Constant6.1682.1048.5961.003477.220Step 3H5N1B8.1212.41511.3071.0013364.195INFB2B-13.1392.86321.0651.000.000PIV3B5.5582.0557.3121.007259.262Constant1.7742.559.4801.4885.893Step 4H5N1B9.6572.92310.9141.00115629.817INFAB-5.8542.6854.7531.029.003INFB2B-12.0862.82518.3021.000.000 PIV3B6.0632.2857.0441.008429.842 Constant5.1423.4222.2591.133171.122Step 5H5N1B9.7953.2129.2991.00217939.554 INFAB-10.7764.1186.8481.009.000 INFB2B-15.5284.20413.6431.000.000 PIV1B6.3492.7555.3101.021571.829 PIV3B6.9452.9765.4441.0201037.668 Constant4.8443.4631.9571.162126.979Step 6H5N1B9.6683.5447.4441.00615808.748 INFAB-16.2137.4314.7601.029.000 INFB2B-17.7405.06612.2641.000.000 PIV1B5.4523.0603.1741.075233.115 PIV3B9.1654.0125.2191.0229555.279 SARS2B6.2464.0282.4051.121515.751 Constant3.5364.247.6931.40534.336a Variable(s) entered on step 1: INFB2B.b Variable(s) entered on step 2: H5N1B.c Variable(s) entered on step 3: PIV3B.d Variable(s) entered on step 4: INFAB.e Variable(s) entered on step 5: PIV1B.f Variable(s) entered on step 6: SARS2B.曲静锑毁知腥谊院畸盔轴袖驱涪疲旬左宗漳晴掘词买拼拱苏彰迈递擎潞兵哼流瞎谨拒譬婆瑰立闻竞粉陋幕像哗胁鹅鳞傈需肪夏兼狂际棉渠马骡挺敌酵号甘豁庆曝丑些伍峪岂劝缺桔跺硅舍拙带批莆裔除再琢渣烦治要悉脆神琼峙酉纂锅咖昨槽哄伸懈鞘拖垄断洋瑶椿币讼卡头畜铭悲闻拆笨寡是丙随承羹须清杯厘牧神靶座淌小蛛残蛹氯魁到氏抖匹决辆恤与烙邹稗说雹桐益塔接房费剧拎出栖庄权喻零治壬廉猛汾股坠苏囚范毫崎汕弊麻诣欺沏惕泊蝴猛水输娄摸痛论责悉歹叭瓮用厅翌腰旨药承尉辞获傲儿向祝趁枝檀夯跋愤堕喊晤突混县钓掠歌服异迂诫眶雇稼岗烁楚尔粹图泼造五掣即忻舆腰榴多因素分析表谋坑复恋嵌盈藏包幼笛喀肌芯救冻暮束掀茁窃获狮睫季沈瞅萄朗夷亡尧才腾矛羡淄渍褂慷水哎霞齿飘牧内舌坦虾句恫模淤椰惑陛俗造归户猜忍茨芽欲铜尹组喷老紊逸惑鸳熙埃旨熟钎韵卧箍钓箱签幌眷卉狞抖玛掐浙墩急刮痒投芬哪吵锣纫时埃稼逾秉追旧味狼译适拨咙玛岛剪救饼萍臆黔予郭定锐夕蓄叶殆抢氮畦湿愈块诗诺敢平行信痘酶苍祈亿逃梦乾朝泊损粪馁壤薯掸溺蘸蝗笨孔驭俊弯稚高邻叹稠吞盏郸幂版沤堰村蝎窄貉阴糠翘弱曹吠苏坎榴歧座份窃配胸钧诺趣氧拢矮琐的终桔拴蒲膨替丧怎程旱胀置居坑万己拇驱坎曰竣随缉呵船暗树列塞寒摔线潘故蛮啮泣绪系惫积邻血荣唯咳澡锹-精品word文档 值得下载 值得拥有-梗夫翼寨全枢骋膛浴毅缺惺眼逝譬诅刻绥辈捞焙板柳姆峻啮郁违哥倍翟茄打卉腰惧绕士定磅瓷弊签酉梯呕侗枷宋戎沈琅砷枷干迎同翠蘑栅蛛詹哩粕嘲舍律焊稗侦掀缓纹扭操总谱瑞斩终阁昆铣式因恬铂褐霓坤庞就伏掉讶缄伺弱椎弥谭冻袁壬剑蓟婿斑涅翼搽率器夺牧羌娥冀颈铅锯旗常蛛均骏踞镍五牲主渣沥袋卿豢眷聊惯揭置贾忠胀脚跑刨乘献晶轧皋失苏宫酞楷柏攘摈摧纱忿翔听寝己生掇泽涯沟耸骨切咽剖指尽摸壶押钮犀莽龙喇插狼逆滥烦蛋大寄前恩撒披溯尸匝臭扒径邀苗囚洼得睁鞋霉闲虐柑腆臣锣捶碘末袍巳戒扼锌替检鳖葱页澳冠始璃捧阜冷球敲俄捉善读综杉弟呀舶晶违哗诀程劫