聚类分析定义及分析方法

  • 格式:pdf
  • 大小:257.61 KB
  • 文档页数:6

聚类分析定义及分析⽅法

聚类分析

聚类分析(Cluster Analysis)是根据事物本⾝的特性研究个体分类的⽅法。

聚类分析的原则是同⼀类中的个体有较⼤的相似性,不同类的个体差异很⼤。

根据分类对象不同分为样品聚类和变量聚类。

样品聚类在统计学中⼜称为Q型聚类。⽤SPSS的术语来说就是对事件(cases)进⾏聚类,或是说对观测量进⾏聚类。是根据被观测的对象的各种特征,即反映被观测对象的特征的各变量值进⾏分类。

变量聚类在统计学中有称为R型聚类。反映事物特点的变量有很多,我们往往根据所研究的问题选择部分变量对事物的某⼀⽅⾯进⾏研究。SPSS中进⾏聚类和判别分析的统计过程是由菜单Analyze---Classify导出的

选择Classify 可以显⽰三个过程命令:1 K-Means Cluster进⾏快速聚类过程。

2 Hierarchical Cluster进⾏样本聚类和变量聚类过程。

3 Discriminant进⾏判别分析过程。

通常情况下在聚类进⾏之前 Proximitice 过程先根据反映各类特性的变量对原始数据

进⾏预处理,即利⽤标准化⽅法对原始数据进⾏⼀次转换。并进⾏相似性测度或距离测度。然后 Cluster 过程根据转换后的数据进⾏聚类分析。

在SPSS for Windows 中分层聚类各⽅法都包含了 Proximitice 过程对数据的处理和Cluster 过程。对数据的分析给出的统计量可以帮助⽤户确定最好的分类结果。1.1 主要功能

聚类的⽅法有多种,最常⽤的是分层聚类法。根据聚类过程不同⼜分为凝聚法和分解法。

分解法:聚类开始把所有个体(观测量或变量)都视为属于⼀⼤类,然后根据距离和相似性逐层分解,直到参与聚类的每个个体⾃成⼀类为⽌。

凝聚法:聚类开始把参与聚类的每个个体(观测量或变量)视为⼀类,根据两类之间的距离或相似性逐步合并直到合并为⼀个⼤类为⽌。

⽆论哪种⽅法,其聚类原则都是近似的聚为⼀类,即距离最近或最相似的聚为⼀类。实际上以上两种⽅法是⽅向相反的两种聚类过程。

调⽤此过程可完成系统聚类分析。在系统聚类分析中⽤户事先⽆法确定类别数,系统将所有例数均调⼊内存,且可执⾏不同的聚类算法。系统聚类分析有两种形式,⼀是对研究对象本⾝进⾏分类称为Q聚类。另⼀种是对研究对象的观察指标进⾏分类称为R型聚类。1.2 实例操作

29名⼉童的⾎红蛋⽩(g/100ml)与微量元素(µg/100ml)测定结果如下表。由于微量元素的测定成本⾼、耗时长,故希望通过聚类分析(即R型指标聚类)筛选代表性指标,以便更经济快捷地评价⼉童的营养状态。

编号钙镁铁锰铜⾎红蛋⽩154.8930.86448.70.012 1.0113.5

272.4942.61467.30.008 1.6413

353.8152.86425.610.004 1.2213.75

464.7439.18469.80.005 1.2214

558.837.67456.550.012 1.0114.25643.6726.18395.780.0010.59412.75

754.8930.86448.70.012 1.0112.5

886.1243.79440.130.017 1.7712.25

960.3538.2394.40.001 1.1412

1054.0434.23405.60.008 1.311.75

1161.2337.354460.022 1.3811.5

1260.1733.67383.20.0010.91411.25

1369.6940.01416.70.012 1.3511

1472.2840.12430.80 1.210.75

1555.1333.02445.80.0120.91810.5

1670.0836.81409.80.012 1.1910.25

1763.0535.07384.100.85310

1848.7530.53342.90.0180.9249.75

1952.2827.14326.290.0040.8179.5

2052.2136.18388.540.024 1.029.25

2149.7125.43331.10.0120.8979

2261.0229.27258.940.016 1.198.75

2353.6828.79292.80.048 1.328.5

2450.2229.17292.60.006 1.048.25

2565.3429.99312.80.006 1.038

2656.3929.292830.016 1.357.8

2766.1231.93344.200.6897.5

2873.8932.94312.50.064 1.157.25

2947.3128.55294.70.0050.8387

1.2.1数据准备

激活数据管理窗⼝,定义变量名:钙、镁、铁、锰、铜和⾎红蛋⽩的变量名分别为x1、x2、x3、x4、x5、x6,之后输⼊原始数据。1.2.2 统计分析

激活Statistics菜单选 Classify 中的 Hierarchical Cluster...项,弹出 Hierarchical Cluster Analysis 对话框。从对话框左侧的变量列表中选x1、x2、x3、x4、x5、x6,点击按钮使之进⼊ Variable(s) 框;在 Cluster 处选择聚类类型,其中 Cases 表⽰观察对象聚类,Variables 表⽰变量聚类,本例选择 Variables。

说明:(1) Variable(s) 栏存放分析变量栏。

(2) Label Cases 栏存放标识变量。

(3) Cluster 栏选择聚类类型。

选择 Variable(s) 项要进⾏变量聚类

选择 Cases 项要进⾏观测量聚类

点击Statistics...钮,弹出Hierarchical Cluster Analysis: Statistics对话框,选择Proximity matrix,要求显⽰距离矩阵,点击Continue钮返回Hierarchical Cluster Analysis对话框。

本例要求系统输出聚类结果的树状关系图,故点击 Plots... 钮弹出 Hierarchical Cluster Analysis:Plots 对话框,选择Dendrogram (树形图)项,点击 Continue 钮返回 Hierarchical Cluster Analysis 对话框。

点击Method...钮弹出Hierarchical Cluster Analysis:Method对话框,系统提供7种聚类⽅法供⽤户选择:Between-groups linkage:类间平均链锁法;合并两类的结果使所有的两两项对之间的平均距离最⼩。项对的两个成员分别属于不同的类。该⽅法中使⽤的是各对之间的距离,即⾮最⼤距离也⾮最⼩距离。Within-groups linkage:类内平均链锁法;若当两类合并为⼀类后,合并后的类中的所有项之间的平均距离最⼩。两类间的距离即是合并后的类中所有可能的观测量对之间的距离平⽅。Nearest neighbor:最近邻居法;该⽅法⾸先合并最近的或最相似的两项,⽤两类间最近点间的距离代表两类间的距离。

Furthest neighbor:最远邻居法;⽤两类之间最远点的距离代表两类之间的距离,也称之为完全连接法。

Centroid clustering:重⼼法,应与欧⽒距离平⽅法⼀起使⽤;像计算所有各项均值之间的距离那样计算两类之间的距离,该距离随聚类的进⾏不断减⼩。Median clustering:中间距离法,应与欧⽒距离平⽅法⼀起使⽤;

Ward's method:离差平⽅和法,应与欧⽒距离平⽅法⼀起使⽤。

本例选择类间平均链锁法(系统默认⽅法)。

在选择距离测量技术上,系统提供8种形式供⽤户选择:Euclidean distance:Euclidean距离,即两观察单位间的距离为其值差的平⽅和的平⽅根,该技术⽤于Q型聚类;

Squared Euclidean distance:Euclidean距离平⽅,即两观察单位间的距离为其值差的平⽅和,该技术⽤于Q型聚类;

Cosine:变量⽮量的余弦,这是模型相似性的度量;

Pearson correlation:相关系数距离,适⽤于R型聚类;

Chebychev:Chebychev距离,即两观察单位间的距离为其任意变量的最⼤绝对差值,该技术⽤于Q型聚类;

Block:City-Block或Manhattan距离,即两观察单位间的距离为其值差的绝对值和,适⽤于Q型聚类;

Minkowski:距离是⼀个绝对幂的度量,即变量绝对值的第p次幂之和的平⽅根;p由⽤户指定

Customized:距离是⼀个绝对幂的度量,即变量绝对值的第p次幂之和的第r次根,p 与r由⽤户指定。本例选⽤Pearson correlation,点击Continue钮返回Hierarchical Cluster Analysis对话框,再点击OK钮即完成分析。

1.2.3 结果解释

得到结果如表

①处理数据的基本信息

29例样本进⼊聚类分析,采⽤相关系数测量技术。

②欧⽒不相似系数平⽅矩阵先显⽰各变量间的相关系数,这对于后⾯选择典型变量是⼗分有⽤的。

③聚类的凝聚过程表

显⽰类间平均链锁法的合并进程,即第⼀步,X3与X6被合并,它们之间的相关系数最⼤,为0.863431;第⼆步,X1与X5合并,其间相关系数为0.624839;第三步,X2与第⼀步的合并项被合并,它们之间的相关系数为0.602099;

第四步,它们与第⼆步的合并项再合并,其间相关系数为0.338335;第五步,与最后⼀个变量X4合并,这个相关系数最⼩,为按类间平均链锁法,变量合并过程的冰柱图如下。先是X3与X6合并,接着X1与X5合并,然后X3、X6与X2合并,接着再与X1、X5合并,最后加上X4,六个变量全部合并。

④聚为五类的冰柱图

⽤更为直观的聚类树状关系图表⽰,即X1、X2、X3、X5、X6先聚合后与X4再聚合。这表明,在评价⼉童营养状态时,可在微量元素钙、镁、铁、铜和⾎红蛋⽩5个指标中选择⼀个,再加上微量元素锰即可,其效果与六个指标都⽤是基本等价的,但更经济更迅速。* * * * * * H I E R A R C H I C A L C L U S T E R A N A L Y S I S * * * * * * Dendrogram using Average Linkage (BetweenGroups)

Rescaled Distance Cluster Combine

C A S E 0 5 10 15 20 25

Label Num +---------+---------+---------+---------+---------+

X3 3 ò?òòòòòòòòòòòòò?

X6 6 ò÷ùòòòòòòòòòòòòò?

X2 2 òòòòòòòòòòòòòòò÷ùòòòòòòòòòòòòòòòòòòò?

X1 1 òòòòòòòòòòòòò?òòòòòòòòòòòòòòò÷ó

X5 5 òòòòòòòòòòòòò÷ó

X4 4 òòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòòò÷