《判别分析实例.doc

例：人文与发展指数是联合国开发计划署于1990年5月发表的第一份《人类发展报告》中公布的。该报告建议，目前对人文发展的衡量指标应当以人生的三大要素为重点。衡量人生的三大要素的指标分别为：实际人均GDP指数、出生时的预期寿命指数、受教育程度指数（由成人识字率指数和综合总人学率指数按2/3、1/3的权重加权而得），将一生三个指数合成为一个指数就是人文发展指数。今从2007年世界各国人文发展指数（2005年）的排序中，选取高发展水平、中等发展水平和低发展水平国家各6个作为三组样品，另选四个国家作为待判样品，资料如下表所示。试用判别分析过程对以下数据资料进行判别分析，并据此对待选的四个国家进行判别归类。 国家 人均GDP（美元） 出生时的预期寿命（岁） 成人识字率（%） 初等、中等和高等教育入学率（%） 第一类：高发展水平国家 美国 41890 77.9 99.5 93.3 德国 29461 79.1 99.2 88 希腊 23381 78.9 96 99 新加坡 29663 79.4 92.5 87.3 意大利 28529 80.3 98.4 90.6 韩国 22029 77.9 99 96 第二类：中等发展水平国家 古巴 6000 77.7 99.8 87.6 罗马尼亚 9060 71.9 97.3 76.8 巴西 8402 71.7 88.6 87.5 泰国 8677 69.6 92.6 71.2 菲律宾 5137 71 92.6 81.1 土耳其 8407 71.4 87.4 68.7 第三类：低发展水平国家 尼泊尔 1550 62.6 48.6 58.1 尼日利亚 1128 46.5 69.1 56.2 喀麦隆 2299 49.8 67.9 62.3 巴基斯坦 2370 64.6 49.9 40 越南 3071 73.7 90.3 63.9 印度尼西亚 3843 69.7 90.4 68.2 待判组 日本 31267 82.3 99 85.9 印度 3452 63.7 61 63.8 中国 6757 72.5 90.9 69.1 南非 11110 50.8 82.4 77 data develop; input type gdp life rate zhrate@@; cards; 1 41890 77.9 99.5 93.3 1 29461 79.1 99.2 88 1 23381 78.9 96 99 1 29663 79.4 92.5 87.3 1 28529 80.3 98.4 90.6 1 22029 77.9 99 96 2 6000 77.7 99.8 87.6 2 9060 71.9 97.3 76.8 2 8402 71.7 88.6 87.5 2 8677 69.6 92.6 71.2 2 5137 71 92.6 81.1 2 8407 71.4 87.4 68.7 3 1550 62.6 48.6 58.1 3 1128 46.5 69.1 56.2 3 2299 49.8 67.9 62.3 3 2370 64.6 49.9 40 3 3071 73.7 90.3 63.9 3 3843 69.7 90.4 68.2 . 31267 82.3 99 85.9 . 3452 63.7 61 63.8 . 6757 72.5 90.9 69.1 . 11110 50.8 82.4 77 ; proc discrim simple wcov distance list;/*simple:要求技术各类样品的简单描述统计量；选项WCOV要求计算类内协方差阵；选项DISTANCE要求计算马氏距离；选项LIST要求输出重复替换归类结果。由于没有给出方法选项，所以系统按缺省时的正态分布进行有关参数的估计和归类。*/ class type; var gdp life rate zhrate; run; proc discrim pool=test slpool=0.05 list; /*simple: */ class type; priors 1=0.3 2=0.4 3=0.3 ; run; proc discrim method=npar k=2 list; /*simple: */ class type; run; proc candisc out=result ncan=2; /*simple: */ class type; var gdp life rate zhrate; run; proc gplot data=reult; plot can1*can2=type; run; proc