判别分析学习教案

1、会计学1判别分析判别分析第一页，共75页。2*重点：判别分析概念、应用、结果(ji gu)解释、注意的问题。第1页/共74页第二页，共75页。3,后果。后果。第2页/共74页第三页，共75页。4判别指标判别指标: 生后生后1 分钟阿氏评分分钟阿氏评分(X1) 、窒息复苏时间、窒息复苏时间(X2 ) 、惊厥持续天数、惊厥持续天数(X3 ) 、急性期、急性期CT 改改变变(X4) 及治疗及治疗3 天后原始反射情况天后原始反射情况(X5).求得判别值求得判别值Z = 77, 以以Z 77 作为预后不良作为预后不良(bling)的标准的标准.不良不良(bling)预后预后:脑性瘫痪脑性瘫痪, 严重智能

2、低严重智能低下下第3页/共74页第四页，共75页。5例例: 根据经济指标根据经济指标,人均收入、人人均收入、人均工农产值、人均消费水平等判断均工农产值、人均消费水平等判断不同地区经济发展程度不同地区经济发展程度(chngd)类型。类型。第4页/共74页第五页，共75页。6第5页/共74页第六页，共75页。7第6页/共74页第七页，共75页。8第7页/共74页第八页，共75页。9第8页/共74页第九页，共75页。101 122.mmzc xc xc x(18-1)第9页/共74页第十页，共75页。11正常人冠心病人22ABABZZSSz1z2Z第10页/共74页第十一页，共75页。12 准则(z

3、hnz)：寻找组间变异(类间均数)/ 组内变异的比值最大化. 1 122.mmzc xc xc x22ABABZZSS22BABAsszz英国(yn u)统计学家FisherR A 爵士(18901962)第11页/共74页第十二页，共75页。13mmmmmmmmmmDCSCSCSDCSCSCSDCSCSCS.22112222212111212111(18-3)Sij为第i指标(zhbio)和第j个指标(zhbio)的合并协方差1 122.mmzc xc xc x类间均数差值第12页/共74页第十三页，共75页。142BAczzz(18-5)cicicizzzzzz判为A类判为B类判为任意(r

4、ny)一类第13页/共74页第十四页，共75页。15第14页/共74页第十五页，共75页。167 . 2 , 8 . 5, 3 . 28 . 5, 2 .38, 3 .203 . 2, 3 .20, 3 .175,33,32,31232221131211ssssssssss计算步骤(bzhu):1.计算各类均数和合并(A、B)的协方差距阵(S)第15页/共74页第十六页，共75页。171)(1)(2211221112BBBBBAAAAAnxxxxnxxxxs1)(1)(21121111BBBAAAnxxnxxs2)(xxlxx)(yyxxlxy第16页/共74页第十七页，共75页。1827 .

5、 28 . 53 . 298 . 52 .383 .2073 . 23 .203 .175321321321CCCCCCCCC318. 0,225. 0,07. 0321CCC1230.070.2250.318zxxx 类间均数差值第17页/共74页第十八页，共75页。19722. 1,428. 1BAzz1230.070.2250.318zxxx 将各类每个个体(gt)的变量值代入判别方程，得到zi，得到zA和zB的均数.004. 02)722. 1(428. 12BAczzz预测: 某病人测定了x1、x2、x3值，代入方程(fngchng)z，计算的z-0.004，为A类。第18页/共74

6、页第十九页，共75页。20z-0.004，为A类第19页/共74页第二十页，共75页。21类的例数类的例数判为BABAp)/(类的例数类的例数判为ABABp)/(BAAB误判为 类类的例数总误判率类类的例数第20页/共74页第二十一页，共75页。22第21页/共74页第二十二页，共75页。23 SPSS统计软件统计软件(run jin)的判别分析的判别分析与结果与结果 (Fisher判别或典则判别判别或典则判别)第22页/共74页第二十三页，共75页。24第23页/共74页第二十四页，共75页。25第24页/共74页第二十五页，共75页。26表1第25页/共74页第二十六页，共75页。27表2

7、第26页/共74页第二十七页，共75页。28本例p0.05,满足(mnz)齐性条件.表3第27页/共74页第二十八页，共75页。29z1z2第28页/共74页第二十九页，共75页。302BAczzz(18-5)cicicizzzzzz判为A类判为B类判为任意(rny)一类第29页/共74页第三十页，共75页。31表4和表5反映建立判别方程提取信息量.和有无(yu w)统计意义表4表5第30页/共74页第三十一页，共75页。320.020.4 1 0.127 20.179 3zxxx 根据表6系数(xsh)建立判别方程表6第31页/共74页第三十二页，共75页。33表表7提供了各变量提供了各变量

8、(binling)对判别分类的重要对判别分类的重要性。性。表7系数(xsh)的绝对值反映重要性第32页/共74页第三十三页，共75页。34表8灵敏度特异(ty)度第33页/共74页第三十四页，共75页。35第34页/共74页第三十五页，共75页。36总误判率低于总误判率低于0.2,认为判别认为判别(pnbi)函数可用函数可用.第35页/共74页第三十六页，共75页。37国家国家类别类别期望寿命期望寿命识字率识字率GDP美国美国1.0076.0099.005374.00日本日本1.0079.5099.005359.00训练样本训练样本瑞士瑞士1.0078.0099.005372.00阿根廷阿根廷

9、1.0072.1095.905242.00阿联酋阿联酋1.0073.8077.705370.00保加利亚保加利亚2.0071.2093.004250.00古巴古巴2.0075.3094.903412.00巴拉圭巴拉圭2.0070.0091.203390.00格鲁吉亚格鲁吉亚2.0072.8099.002300.00南非南非2.0062.9080.603799.00中国中国2.0068.5079.301950.00待判样本待判样本罗马尼亚罗马尼亚2.0069.9096.902840.00希腊希腊1.0077.6093.805233.00哥伦比亚哥伦比亚1.0069.3090.305158.00第

10、36页/共74页第三十七页，共75页。38)/ )()./ )()/ )(21klkklklkysxpysxpysxpp (18-7)第37页/共74页第三十八页，共75页。39第38页/共74页第三十九页，共75页。40第39页/共74页第四十页，共75页。4136.6白细胞白细胞 x7 23.7(单位单位)第40页/共74页第四十一页，共75页。4200015. 057. 010. 0.65. 021. 00047. 061. 018. 0.55. 034. 00018. 028. 032. 0.37. 034. 00017. 008. 061. 008. 095. 072. 011. 0

11、57. 04321pppp该病例预测(yc)为第3类 - 坏疽型该法主要得到 表18-3条件(tiojin)概率,要求例数较多.第41页/共74页第四十二页，共75页。431l12l2lmg1l12l2lmk=1() ()/) ()/)()/)(/ )() ()/) ()/)()/)KkkmkkKkkmkP YP X SY P XSYP XSYP YaP YP X SY P XSYP XSY()kp y: 第K类的先验概率,在事件(shjin)中的比例.第42页/共74页第四十三页，共75页。44()kp y第43页/共74页第四十四页，共75页。4510.2 0.57 0.11 0.72 0

12、.95 0.08 0.61 0.080.000033p 10.000030(/ )0.0150.000033 0.00090.0011750.000075p ya 2(/ )0.412p ya 3(/ )0.538p ya 4(/ )0.034p ya 结果(ji gu)与前法相同第44页/共74页第四十五页，共75页。46第45页/共74页第四十六页，共75页。47mmGGGmmmmxcxcxccpyxcxcxccpyxcxcxccpy.ln.ln.ln111101111022211110111已知各类的先验概率Gp第46页/共74页第四十七页，共75页。48mmGGmmmmxcxcxccy

13、xcxcxccyxcxcxccy.1111011110211111011（18-9）gPK1其中C为判别(pnbi)系数。系数c的与fisher计算相同.第47页/共74页第四十八页，共75页。49gtclckkyyyyp1)exp()exp(预测方法(fngf):1.将个体判为YG值中最大的类.2. 计算个体各类后验概率,判为概率最大类.两方法(fngf)结果等价(公式18-13)第48页/共74页第四十九页，共75页。50第49页/共74页第五十页，共75页。51第50页/共74页第五十一页，共75页。5249.149252.10206.14118.79152.102662.26836.1

14、395.2906.141836.13388.15356.64188.79195.29356.6446.2074S X1 X2 X3 X42. 按讲义公式(gngsh)18-11计算出各系数C第51页/共74页第五十二页，共75页。53311)()()221gypypyP（mjkjJKkkXCypC1021)(lg第52页/共74页第五十三页，共75页。5443366. 137202. 424568.1710396. 0094.19044140. 1366. 42097.1810448. 053.1994582. 135492. 4241.1910739. 051.223321xxxxyxxxx

15、yxxxxy应用(yngyng): 将个体的m个变量值代入上面3个判别式,个体归为最大y值类.见表18-4.第53页/共74页第五十四页，共75页。55后验概率后验概率X1X2X3X4原分类原分类1类类2类类3类类判别结果判别结果6.0-11.519901.981.018.0001-11.0-18.525-363.000.140.860390.2-17.01732.002.547.4512-4.0-15.013541.969.030.0011第54页/共74页第五十五页，共75页。56第55页/共74页第五十六页，共75页。57第56页/共74页第五十七页，共75页。58rrrTw(18-14

16、)11grgNFFF r指变量(binling)个数第57页/共74页第五十八页，共75页。59第58页/共74页第五十九页，共75页。6011129042.470.97296520.28wT 222215.43800.99215.4612wT 333373.26190.66558.2353wT 44420894.860.744328070.23wT 见395页的W和T距阵第59页/共74页第六十页，共75页。61717. 32031766. 066. 0111grgNF第二步，X4选入，F=5.714第三步，X2选入，F=2.192第四步，X3剔除，F=0.1174.方程内选入X4、X2，建

17、立(jinl)方程结果见396页.74. 314, 2 ,05. 0F第60页/共74页第六十一页，共75页。62SPSS统计软件的逐步统计软件的逐步(zhb)判别分析与判别分析与结果结果 (Bayes判别法判别法)第61页/共74页第六十二页，共75页。63逐步(zhb)判别分析第62页/共74页第六十三页，共75页。64选入和剔除(tch)方程的概率用0.2和0.3第63页/共74页第六十四页，共75页。65变量(binling)在方程内的情况第64页/共74页第六十五页，共75页。66123101.489.8652 20.9533 474.928.4737 20.8009 462.767.7397 20.7215 4yxxyxxyxx Bayes判别函数方程(fngchng)第65页/共74页第六十六页，共75页。67第66页/共74页第六十七页，共75页。68第67页/共74页第六十八页，共75页。69给出判别后个体(gt)是各类的概率和判别结果。判为第1、第2、第3的概率(gil)。第68页/共74页第六十九页，共75页。70对判别的指标要求高对判别的指标要求高.第69页/共74页第七十页，共75页。71第70页/共74页第七十一页，共75页。72第71页/共74页第七十二页，共75页。73给出预测后的后验概率根据方程平均分界点(Z)根据G类方程得到后验