判别分析实例汇总

例：人文与发展指数是联合国开发计划署于1990年5月发表的第一份《人类发展报告》中

公布的。该报告建议，目前对人文发展的衡量指标应当以人生的三大要素为重点。衡量人生

的三大要素的指标分别为：实际人均GDP指数、出生时的预期寿命指数、受教育程度指数

（由成人识字率指数和综合总人学率指数按2/3、1/3的权重加权而得），将一生三个指数合成为一个指数就是人文发展指数。今从2007年世界各国人文发展指数（2005年）的排序中, 选取高发展水平、中等发展水平和低发展水平国家各6个作为三组样品，另选四个国家作为待判样品，资料如下表所示。试用判别分析过程对以下数据资料进行判别分析，并据此对待

选的四个国家进行判别归类。

data develop;

in put type gdp life rate zhrate@@;

cards ;

1 4189077.999.593.3

1 2946179.199.288

1 2338178.99699

1 2966379.492.587.3

1 2852980.398.490.6

1 2202977.99996

2 6000 77.799.887.6

2 9060 71.997.376.8

2 8402 71.788.687.5

2 8677 69.692.671.2

2 5137 7192.681.1

2 8407 71.487.468.7

3 1550 62.648.658.1

3 1128 46.569.156.2

3 2299 49.867.962.3

3 2370 64.649.940

3 3071 73.790.363.9

3 3843 69.790.468.2

.3126782.39985.9

.3452 63.76163.8

.6757 72.590.969.1

.1111050.882.477

proc discrim simple wcov dista nee list ; /*simple: 要求技术各类样品的简单描述

统计量；选项WCOV要求计算类内协方差阵；选项DISTANCE要求计算马氏距离；选项LIST要求输出重复替换归类结果。由于没有给出方法选项,所以系统按缺省时的正态分布进行有关参数的

估计和归类。*/

class type;

run ;

proc discrim poo匸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 can disc out =result ncan =2; /*simple: */

class type;

var gdp life rate zhrate;

run ;

proc gplot data =reult; ―|

plot can 1*ca n2=type;

run ;

proc discrim data =result dista nee list ;

class type;

var canl can2;

run ;

表1已知样本分类水平信息

The DISCRIM Procedure

Observations Variables Classes

18

4

3

DF To怙1

DF Within Classes

DF Between Classes

17

15

2 Class Level Informal ion

Variable Prior type Name Frequency Weight Proportion Probabi1 ity 1_16G.00000.$888380.333333 2J6 6.00000.3333330.333333

36G.00000.3333330.333333表2样本统计量信息

type = 1, DF = 5

Variable gdp life rate zhrate gdp49408532.97-1234.124172.07-11022.03 1 ife-1234.120.85-0.88-2.09

rate4172.07-0.887.43 2.74 zhrate-11022.03-2.09 2.7421.19

type = 2, DF = 5

Variable gdp life rate zhrate gdp2642240.567-2026.117•2419.950-6404.957 1 ife-2026.1177.8868.86113.946 rate-2419.9508.86123.15114.327 zhrate-6404.95713.94614.32764.438

type = 3, DF = 5

Variable gdp life rate zhrate gdp976170.96677840.770012624.07334200.8033 1 ife7840.7700117.611073.166015.3730 rate12624.073373.1660338.6067136.1087 2hrate4200.803315.3730136.108796.9017