数学建模案例分析消费分布规律的分类概率统计方法建模

由天下分享时间：2025/2/23 21:30:20 加入收藏我要投稿点赞

§7 消费分布规律的分类

为研究辽宁、浙江、河南、甘肃、青海5省份在某年城镇居民生活消费的分布规律，需要用调查资料对这5个省分类.数据见下表：指标省份辽宁浙江河南甘肃青海 X1 X2 X3 X4 X5 X6 X7 X8 7.90 39.77 8.49 12.94 19.27 11.05 2.04 13.29 7.68 50.37 11.35 13.30 19.25 14.59 2.75 14.87 9.42 27.93 8.20 8.14 16.17 9.42 1.55 9.76 9.16 27.98 9.01 9.32 15.99 9.10 1.82 11.35 10.06 28.64 10.52 10.05 16.18 8.39 1.96 10.81 其中，X1：人均粮食支出； X2：人均副食品支出； X3：人均烟、酒、茶支出； X4：人均其它副食品支出； X5：人均衣着商品支出； X6：人均日用品支出； X7：人均燃料支出； X8：人均非商品支出.

在科学研究、生产实践、社会生活中，经常会遇到分类的问题.例如，在考古学中，要将某些古生物化石进行科学的分类；在生物学中，要根据各生物体的综合特征进行分类；在经济学中，要考虑哪些经济指标反映的是同一种经济特征；在产品质量管理中，要根据各产品的某些重要指标而将其分为一等品，二等品等等.

这些问题可以用聚类分析方法来解决.

聚类分析的研究内容包括两个方面，一是对样品进行分类，称为Q型聚类法，使用的统计量是样品间的距离；二是对变量进行分类，称为R型聚类法，使用的统计量是变量间的相似系数.

设共有n个样品，每个样品xi有p个变量，它们的观测值可以表示为 xi?(x1i,x2i,?,xpi),i?1,2,?,n 一、样品间的距离

下面介绍在聚类分析中常用的几种定义样品xi与样品xj间的距离. 1、 Minkowski 距离 d(xi,xj)?[?xk?1kipmki?xkj]

1m2、绝对值距离 d(xi,xj)??xk?1p?xkj

2123、欧氏距离 d(xi,xj)?[?xk?1pki?xkj]

二、变量间的相似系数

相似系数越接近1，说明变量间的关联程度越好.常用的变量间的相似系数有 1、夹角余弦

1 / 4

rij??xk?1nikxjkn?(xk?1n

ik)2??(xjk)2k?12、相关系数

rij??(xk?1nik?x(i))(xjk?x(j))n?(xk?1n

ik?x(i))2??(xjk?x(j))2k?1值得注意的是，当指标的测量值相差较大时，直接使用以上各式计算距离或相似系数常使数值较小的变量失去作用，为此需应先对数据进行标准化，然后再用标准化的数据来计算.标准化的具体方法是：

xki?*xki?xk,i?1,2,?,n,k?1,2,?,p sk1n1n其中 xk??xki,sk?(xki?xk)2,k?1,2,?,p ?ni?1n?1i?1三、类与类之间的距离

用Gp和Gq分别代表两个类，它们所包含的样品个数分别记为np和nq，类Gp和Gq之间的距离记为D(Gp,Gq).下面给出三种最常用的定义方法. 1、最短距离 D(Gp,Gq)?min(dijxi?Gp,xj?Gq)

类与类之间的最短距离有如下的递推公式，设Gr为由Gp和Gq合并所得，则Gr与其它类

Gk(k?p,q)的最短距离为

D(Gr,Gk)?min{D(Gp,Gk),D(Gq,Gk)}

2、最长距离 D(Gp,Gq)?max(dijxi?Gp,xj?Gq)

类与类之间的最长距离有如下的递推公式，设Gr为由Gp和Gq合并所得，则Gr与其它类

Gk(k?p,q)的最长距离为

D(Gr,Gk)?max{D(Gp,Gk),D(Gq,Gk)}

3、类平均距离 D(Gp,Gq)?1npnqxi?Gpxj?Gq??dij

2 / 4

类与类之间的类平均距离有如下的递推公式，设Gr为由Gp和Gq合并所得，则Gr与其它类

Gk(k?p,q)的类平均距离

D(Gr,Gk)?npnrD(Gp,Gk)?nqnrD(Gq,Gk)，其中 nr?np?nq.

以上类与类之间的距离，不但适用于Q型聚类，同样也适合于R型聚类，这只要将dij用变量间的相似系数rij代替就行了.为简单起见以下均记成dij.

系统聚类法是目前最流行的方法.

有了样品间的距离（或变量间的相似系数）以及类与类之间的距离后，便可进行系统聚类，基本步骤如下：

1、n个样品（或p个变量）一开始看作n类（p类），计算两两之间的距离（或相似系数），构成一个对称矩阵D0?dij??n?n，此时显然有D(Gp,Gq)?dpq；

2、选择D0中对角线元素以外的下三角部分中的最小元素（相似系数矩阵则选择对角线元素以外的最大者），设其为D(Gp,Gq)，则将Gp和Gq合并为一个新类Gr.在D0中划去Gp和Gq所对应的两行与两列，并加入由新类Gr与剩下的未聚合的各类之间的距离所组成的一行和一列，得到一个新的矩阵D1，它是降低了一阶的对称矩阵；

3、由D1出发，重复步骤2得到对称矩阵D2，依此类推，直到n个样品（或p个变量）聚为一个大类为止；

4、在合并过程中记下两类合并时样品（或变量）的编号以及合并两类时的距离（或相似系数）的大小，并绘成聚类图，然后可根据实际问题的背景和要求选定相应的临界水平以确定类的个数.

上面是一个Q型聚类问题，现在用系统聚类法来解决.将每个省份看成一个样品，并以1，2，3，4，5分别表示辽宁、浙江、河南、甘肃、青海5省，计算样品间的欧氏距离，得到如下的距离矩阵D0

{1} {2} {3} {4} {5}

?0???11.670??? D0??13.8024.630???13.1224.062.200??12.8023.543.512.210???下面给出采用最短距离法的聚类过程：首先将5个省各看成一类，即令Gi?{i},(i?1,2,3,4,5).从D0可以看出，其中最小的元素是D({4},{3})?d43?2.20，故将G3和G4合并成一类G6，然

3 / 4

后利用递推公式计算G6与G1，G2，G5之间的最短距离.

D({3,4},{1})?min{d31,d41}?min{13.80,13.12)?13.12 D({3,4},{2})?min{d32,d42}?min{24.63,24.06)?24.06 D({3,4},{5})?min{d35,d45}?min{3.51,2.21)?2.21

在D0中划去{3},{4}所对应的行和列,并加上新类{3,4}到其它类距离作为新的一行一列,得到 {3,4} {1} {2} {5}

?0???13.120??D1??? 24.0611.670???2.2112.8023.540???重复上面的步骤，依次可得到相应的距离矩阵如下：

{3,4,5} {1} {2}

?0???D2??12.800?,?23.5411.670???{3,4,5} {1,2}

?0D3???12.80??? ?0?最后将5个省合并为一大类,画出聚类图如下：

辽宁

11.67 浙江

12.80 河南

2.20 甘肃