聚类分析及判别分析实验报告范例.docx

1、聚类分析及判别分析实验报告范例 上海电力学院应用多元统计分析判别分析与聚类分析学院： 姓名： 学号： 2016年4月 我国部分城市经济发展水平的聚类分析和判别分析摘要：本文基于中国统计年鉴（2012年版）统计数据，寻找评价城市经济发展水平的指标，包括第二三产业发展水平、固定投资额、社会消费零售总额和进出口贸易交流五个指标，利用统计软件SPSS综合考虑各指标，对所选城市进行K-Means 聚类分析，利用Fisher 线性判别待判城市类型，进一步验证所建模型的有效性。关键字：聚类分析，判别分析，SPSS ，城市经济发展水平1，引言经过改革开放后三十多年的长足进展，中国城市化已步入中期阶段，步伐加快

2、，质量显著提高。同时，中国城市化又处于期转折点上，上一期行将结束，下一期将要开始。2011年中国城市化率首次突破50，意味着中国城镇人口首次超过农村人口，中国城市化进入关键发展阶段，这必将引起深刻的社会变革。根据2011年4月公布的第六次人口普查数据，2010年中国居住城镇的人口接近6.6亿人，城镇化率达到49.68%，全国已有近一半的人口居住在城镇，这意味着中国将进入城镇时代。在过去30多年中，中国的城市化发展取得了很大成绩。然而，总体上中国的城市化道路是城市化滞后于工业化的非均衡道路；是土地城市化快于人口城市化的非规整道路；是以抑制农村、农业、农民的经济利益来支持城市发展，导致不能兼顾效率

3、和公平的非协调道路；是片面追求城市发展的数量和规模，而以生态环境损失为代价的非持续道路；是以生产要素的高投入，而不是投入少、产值高、依靠科技拉动经济增长的非集约道路。传统的城市化存在着诸多弊端，中国未来的城市化必须走出一条具有自身特色的新型城市化道路。具体而言，中国城市经济发展水平受限于地理、环境、资源以及政策等因素的影响，我国不同区域的城市化进程尚存在很大差异。2012年中国城市发展报告中指出，从区域角度看，目前沿海一带城市发展起步早，与国际贸易交流往来频率高，经济发展水平较高，西部地区受到政策的大力扶持，表现出了强劲的增长势头，西部主要城市经济发展水平仅次于沿海发达地区，而中部地区城市发展

4、的水平已经落到了最后。显然，通过研究不同城市的经济发展状况和经济类型，指出其发展差异所在，可以为政府在出台相关政策来平衡区域经济发展，缩小不同地区人民生活水平的差异提供一定的指导意见，也为我们深刻理解相关政策提供了扎实的基础。 2，相关统计基础理论2.1，聚类分析聚类分析指将物理或抽象对象的集合分组成为由类似的对象组成的多个类的分析过程。聚类分析的目标就是在相似的基础上收集数据来分类。从统计学的观点看，聚类分析是通过数据建模简化数据的一种法。传统的统计聚类分析法包括系统聚类法、分解法、加入法、动态聚类法、有序样品聚类、有重叠聚类和模糊聚类等。采用k-均值、k-中心点等算法的聚类分析工具已被加入

5、到多著名的统计分析软件包中，如SPSS、SAS等。本文使用统计软件SPSS对所收集的数据进行快速聚类，其特点是：在确定类别数量基础上，先给定一个粗糙的初始分类，然后按照某种原则进行反复修改，直至分类较为合理。在选定类中心作为凝聚点的基础上进行分类和修正的法有很多，本文使用的是K-Means 算法。K-Means 算法接受输入量k ；然后将n个数据对象划分为 k个聚类以便使得所获得的聚类满足：同一聚类中的对象相似度较高；而不同聚类中的对象相似度较小。聚类相似度是利用各聚类中对象的均值所获得一个“中心对象”（引力中心）来进行计算的。K-Means 算法的工作过程说明如下：首先从n个数据对象任意选择

6、 k 个对象作为初始聚类中心；而对于所剩下其它对象，则根据它们与这些聚类中心的相似度（距离），分别将它们分配给与其最相似的（聚类中心所代表的）聚类；然后再计算每个所获新聚类的聚类中心（该聚类中所有对象的均值）；不断重复这一过程直到标准测度函数开始收敛为止。一般都采用均差作为标准测度函数。 一般而言，k个聚类具有以下特点：各聚类本身尽可能的紧凑，而各聚类之间尽可能的分开。2.2，判别分析判别分析是市场研究的重要分析技术，也是多变量分析技术。它可以就一定数量的个体的一个分类变量和相应的其它多元变量的已知信息，确定分类变量与其它多元变量之间的数量关系，建立判别函数，并利用判别函数构建Biplot二元

7、判别图（概念图）。同时，利用这一数量关系对其他已知多元变量的信息、但未知分组的子类型的个体进行判别分组。判别分析属于监督类分析法，例如：市场细分研究中，常涉及判别个体所属类型的问题，也常涉及不同品牌在一组产品属性之间的消费者偏好和认知概念，判别分析可以很好地对这种差异进行鉴别。并在低维度空间表现这种差异。判别分析主要有距离判别、贝叶斯（Bayes）判别、费舍尔（Fisher）判别等几种常用法。距离判别的基本原理是：首先对样本到总体G之间的距离进行合理规定，然后依照“就近”原则判定样本的归属，常用马氏距离（Mahalanobis）规定为：式中为p元总体G的协差阵，x是取自G的样品，则该式即为样品

8、x到总体G的马氏距离。贝叶斯判别既考虑了先验分布产生的影响，也考虑到误判损失产生的影响，是衡量一个判别优劣的比较合理的准则。费舍尔判别的基本思想与主成分分析十分相似，当总体是高维向量时，先把其综合成一个一维变量，然后在对一维变量进行距离判别，费舍尔判别实际上是一种降维处理，降维压缩后，样品y到各个总体的距离可以用欧式距离度量，即：由此导出Fisher判别规则为：， 则本文及使用Fisher判别建立线性判别函数进行距离判别。3，模型建立3.1 设置变量本文综合考虑了评价城市发展指数衡量因素，选取衡量一个城市经济发展水平的主要因素，城市化进程总是伴随着工业化发展，发达的服务业水平是衡量现代新兴城市

9、的主要指标，此外，综合考虑了固定资产投资总额与社会消费品零售总额以及货物进出口总额作为类别分析的主要经济指标：X1：城市第二产业产值（亿元）X2：城市第三产业产值（亿元）X3：城市固定资产投资总额（亿元）X4：城市社会消费品零售总额（亿元）X5：货物进出口总额（亿元）从区域发展角度从上面5个经济指标将城市经济发展水平划分为三大类：G1：发达城市G2：中度发达城市G3：欠发达城市3.2 数据收集和整理本文所有数据来源于中国统计年鉴（2012），选取2011年度36个城市主要经济发展水平做模型建立及分析。其中前32个城市相关经济指标水平作为初始样本用于划分类别，建立类别总体G；最后四个城市（、）及

10、其相关经济发展水平用作待判样品，利用判别函数进行判别分析。所有相关数据经过量纲统一规则化处理见表1所示。表 1 我国部分城市相关经济发展水平（2011年）序号城市第二产业（亿元）第三产业（亿元）固定资产投资总额（亿元）社会消费品零售总额（亿元）货物进出口总额（亿元）1 北 京 3752.512363.25851.52016900.324623374.98842 天 津 5928.35219.27483.69733395.066203.46423 家 庄 2031.91635.83026.97781662.9864850.11124 太 原 949.21097.11024.1444973.293

11、7513.63065 呼和浩特 790.01277.81031.6781890.0478121.47366 阳 3026.92609.84577.0942426.8655637.2157 大 连 3204.22550.74580.05851924.7943630.58748 长 春 2092.71620.22356.61891515.85371040.93229 哈 尔 滨 1647.22147.83011.9712070.4129307.054810 上 海 7927.911142.95064.26246814.826246.15111 南 京 2760.83220.43757.251726

12、97.09973440.635812 宁 波 3349.52454.52385.50722018.86175891.209213 合 肥 2002.21426.23376.96521111.11881207.71914 福 州 1711.21700.12720.28271947.81022083.485615 厦 门 1297.11217.51128.0872800.27794210.000216 南 昌 1579.3974.72022.3297928.3438473.022617 济 南 1829.02339.51934.33892114.2868624.12318 青 岛 3150.731

13、58.53502.53822302.37034329.130219 州 2874.21974.03002.51987.1147959.735420 武 汉 3254.03309.54255.16213031.78851367.374821 长 沙 3151.72224.33510.24252201.6112449.360422 广 州 4577.07641.93412.252436970.2623 深 圳 5343.36155.72136.38823520.873624845.98224 海 口 177.9487.7395.0408387.1804236.175625 重 庆 5543.036

14、23.87579.44543487.8071753.071626 成 都 3143.83383.44944.01572861.28352274.379827 贵 阳 586.8733.71600.5898584.3292392.979628 拉 萨75.2137.2220.5031102.594878.445229 西 安 1697.21993.93352.121965.9774754.7430 兰 州 656.5663.5950.5758639.7231112.765831 西 宁 411.3332.0528.0052271.287348.937832 乌木齐 759.1908.9427.6

15、221695.0278541.790433 杭 州 3323.83458.53100.02182548.35993838.30834 南 宁 829.61076.31950.86281073.1541150.625235 昆 明 1161.21214.62275.52861271.7298721.322436 银 川 525.2414.4720.5627274.470572.64，数据结果及分析4.1 聚类分析4.1.1 聚类分析过程采用统计软件SPSS可以快速便的将样本分类，“K-均值聚类”将样本分为设定好的三类，分类结果如下：（1）K-均值聚类初始聚类中心初始聚类中心聚类123第二产业（亿

16、元）7928592875第三产业（亿元）111435219137固定资产投资总额（亿元）5064.17483.1220.0社会消费品零售总额（亿元）6814.80003395.0600102.5948货物进出口总额（亿元）26246.06203.078.0（2）样本聚类 聚类成员案例号城市聚类距离案例号城市聚类距离1 北 京13937.77217 济 南31347.1542 天 津24379.85018 青 岛21710.0433 家 庄31259.02619 州31969.2614 太 原31214.20 武 汉22771.8345 呼和浩特31414.69721 长 沙32607.5836

17、 阳33452.67422 广 州25518.2357 大 连21842.87323 深 圳14887.3768 长 春3837.81124 海 口32474.7509 哈 尔 滨31584.29125 重 庆24072.60110 上 海13214.67326 成 都21942.91011 南 京21681.20527 贵 阳31402.62012 宁 波23455.97928 拉 萨32918.19013 合 肥31536.88129 西 安31652.62514 福 州31682.56330 兰 州31767.33415 厦 门33577.16931 西 宁32433.50316 南 昌

18、3617.36732 乌木齐31898.368（3）最终聚类中心最终聚类中心聚类123第二产业（亿元）567538791467第三产业（亿元）988738401375固定资产投资总额（亿元）4350.84655.02009.0社会消费品零售总额（亿元）5745.33272995.78501263.8072货物进出口总额（亿元）24822.03984.0782.0最终聚类中心间的距离聚类123121946.79726337.272221946.7975669.843326337.2725669.843（4）聚类差分析差分析聚类误差FSig.均df均df第二产业（亿元）34313207.73521

19、231856.4792927.855.000第三产业（亿元）100446019.01321811059.4072955.463.000固定资产投资总额（亿元）24862358.67322205819.3762911.271.000社会消费品零售总额（亿元）30454986.0502887338.5312934.322.000货物进出口总额（亿元）753836973.38321848036.99229407.912.000F 检验应仅用于描述性目的，因为选中的聚类将被用来最大化不同聚类中的案例间的差别。观测到的显著性水平并未据此进行更正，因此无法将其解释为是对聚类均值相等这一假设的检验。4.1.

20、2 聚类结果分析从上述聚类分析过程可知，样本完全有效，32个个体被分成三大类：G1（发达城市）：北京，上海，。G2（中度发达城市）：天津，程度。G3（欠发达城市）：，呼和浩特，乌木齐。从城市分类结果可知，北上深作为国际化城市发展代表，其经济发展水平远超其他沿海城市及陆城市；沿海开放城市以及陆主要枢纽城市的发展水平高于其他城市；中部地区级西部城市发展水平受限于地理、资源和资本等因素，经济发展表现不强劲。从最后的差分析中可知，分类检验水平显著，分类结果值得借鉴。4.2 判别分析4.2.1 判别结果及分析一般来讲，利用判别分析首先要明确变量测量尺度及变量的类型和关系；因变量（dependent va

21、riable)： 分组变量定性数据（个体、产品/品牌、特征，定类变量）。自变量 (independent variable)：判别变量定量数据（属性的评价得分，数量型变量）。（1）判别图图 1 典则判别函数从图中很明显，看到三个组中心也就是经济发展水平，以及围绕着组中心的样本，说明直观上分组判别式完全可以接受的。（2）典型判别式函数摘要特征值函数特征值差的 %累积 %正则相关性137.790a98.098.0.9872.765a2.0100.0.658a. 分析中使用了前 2 个典型判别式函数。Wilks 的 Lambda函数检验Wilks 的 Lambda卡dfSig.1 到 2.01511

22、4.10610.0002.56715.3364.004标准化的典型判别式函数系数函数12第二产业（亿元）-.974.940第三产业（亿元）-1.198.773固定资产投资总额（亿元）.752.211社会消费品零售总额（亿元）1.190-.675货物进出口总额（亿元）1.409-.314结构矩阵函数12货物进出口总额（亿元）.863*.090第二产业（亿元）.184.920*固定资产投资总额（亿元）.076.857*社会消费品零售总额（亿元）.231.675*第三产业（亿元）.308.567*判别变量和标准化典型判别式函数之间的汇聚组间相关性 按函数相关性的绝对大小排序的变量。*. 每个变量和任

23、意判别式函数间最大的绝对相关性从表中我们看到，因为分组变量是三类，所以我们得到两个判别函数，其中第一判别函数解释了数据的98%，第二判别函数解释了2%；两个判别函数解释了100%；当然，两个判别函数直接具有显著的差异和判别力。（3）分类统计量组的先验概率类别先验用于分析的案例未加权的已加权的1.33333.0002.33399.0003.3332020.000合计1.0003232.000分类函数系数类别123第二产业（亿元）-.019-.001-.001第三产业（亿元）-.-.004-.003固定资产投资总额（亿元）.011.002.001社会消费品零售总额（亿元）.030.006.004货

24、物进出口总额（亿元）.003.001(常量)-231.519-12.269-2.727Fisher 的线性判别式函数Fisher线性判别函数，我们主要用来构建判别程，理论上说：如果我们知道某个城市在5经济指标的发展水平值，我们就可以估计出该城市应该是哪种类型的。4.2.2 判别检验判别变量是数量型测量尺度变量，分析样本个数至少比判别变量多两个，我们为了得到判别函数，经常需要把样本随机分成训练样本和检验样本等工作，如本文最后四个（序号33-36）个体就可作为检验样本，也成待判样本。由上表可知Fisher判别程：判别规则：， 则判别结果：1 杭 州21583.3912 南 宁3842.7743 昆

25、 明3401.5674 银 川32095.787直观上，作为沿海省会城市，虽然达不到北上深的经济发展地位，但其良好的地理位置以及投资开发环境，使得其经济发展水平非常迅速，归类到第二类经济发展城市是可以理解与接收的。其余三个城市虽然都是省会城市，但都属于西部城市，地理位置以及资源物产相对欠缺，得力于西部大开发政策影响，经济发展面距第三类城市相近。5，结论从本文所建立的模型对我国部分主要城市经济发展水平进行了聚类分析与判别分析，并运用模型判断最后四个城市，验证模型的有效性。从相关结果及分析可以得到一些直观的结论。北京作为我国首都，毋庸置疑具有天然的发展优势，其政治中心，交通中心，文化中心的地位吸引

26、了国外大量的投资建设，一批高科技产业带动的行业发展极大的推动了北京的经济发展；上海作为中国的经济金融中心，加之其边江浙地带发达的工业基础，都为上海的经济发展增加了强劲的力量；的发展是中国改革开放以来经济发展的典，开放的力量使得这个沿海城市一举成为中国发达城市的先驱。判别图里清晰的表明北上广的发展远远超过二三类型的城市发展水平。沿海主要城市以及陆枢纽城市的发展得益于丰富的资源以及便利交通带来的大量投资，或者传统的工业基础，这些因素都使得这一类的城市发展迅速，势头强劲。其余中西部城市的发展各有其优劣，但总体上西部城市受西部大开发政策影响，变现出新兴的发展势头。判别图分析可见二三类型经济发展水平相差不大。由此，所建立模型直观上符合我国部分主要城市经济发展水平类型，最后的四个城市判别再次说明了模型的有效性。参考文献1海燕,梦,卫国,伟. 应用数理统计M. 北京:北京航空航天大学数学系, 2011.2建同,昌言. 以Excel和SPSS为工具的管理统计M. 北京:清华大学出版社，2002.3 统计局.2012年中国统计年鉴M. 中国统计出版社，20124 中国城市经济学会. 中国城市经济J.北京：中国社会科学院