基于全国主要城市平均气温的统计分析.docx

1、基于全国主要城市平均气温的统计分析基于全国主要城市平均气温的统计分析摘要 本文主要采用2012年全国主要城市平均气温的统计数据，首先通过相关分析对全国主要城市的平均气温进行简单的相关分析，然后在控制平均气温的基础上对个月平均气温进行偏相关分析；其次运用聚类分析对各相似省份进行聚类分析；最后运用因子分析对各城市进行因子分析，计算出每个城市的综合得分，通过综合得分计算出每个城市平均气温的排名。关键词：相关分析；聚类分析；因子分析1、 研究背景及目的全球气温变暖为世人所瞩目，近几十年来的全球气候变暖是一个国内外都非常关注的重大问题。自从全球气温变暖的议题出现以来，关于气温的不正常变化仁者见仁，智者见

2、智。气温的变化对农作物，人们的生活及经济的运行都有重要的影响。全世界的气温研究工作者都研究出了许多关于气温的结论，但都有一个共性即全球气温变暖。然而我国跨北纬4度到北纬53度，从东经73度到东经135度，具有比较丰富的气候和比较多变的气温，因此对我国气温的研究具有很重要的意义。在此大背景下，对我国气温的研究不仅能说明我国气温的分布和变化，同时对世界气温的研究同样具有重要作用。本文的研究目的如下：通过对我国主要城镇平均气温的研究和分析，一方面可以得出我国各市平均气温与全国平均气温差异；另一方面也可以给我们以后的就业进行指导，在适宜的季节去适宜的城市工作。2、研究方法气温的变化不仅能促进经济的发展

3、，同时也能阻碍经济的发展。据此我们通过全国各主要城市每个月份的平均气温对我国气温的分布情况进行分析。本研究采用的数据是中国2012年主要城市平均气温的数据指标，数据摘自中国统计年鉴20137-7.采用的分析方法主要有相关分析，聚类分析，因子分析等。基本思路是：首先利用线图分析各个月份平均气温的分布，然后利用条形图分析全国主要城市平均气温的对比，得出平均气温最高和最低的城市；然后利用因子分析对各个月份的平均气温进行相关分析；其次利用聚类分析对各个城市在不同气温下的类别；再次利用因子分析对构成平均气温的各各月平均气温提取公因子；最后使用一些简单的SPSS数据处理技巧依据提出的公因子对各城市进行分类

4、及排序。3、实证分析由于从国家统计局网站下载的数据为EXCEL格式，可以将数据导入成SPSS数据，我们共设置了14个变量，分别是“城市”、“1月”、“2月”、“3月”、“4月”、“5月”、“6月”、“7月”、“8月”、“9月”、“10月”、“11月”、“12月”和“年平均气温”。 样本是中国2012年全国重要城市平均气温的主要数据。数据参见附表1。3.1、画基本图形图1通过图1可以看出全国平均的月份气温最低在0度左右，最高在25度左右，说明全国平均月份气温比较温和。图2条形图反映的是各个城市的年平均气温的情况，最高的为海口，最低的为哈尔滨，比较符合现实的气温分布。3.2、相关分析对于相关分析主

5、要有以下几部分；第一、 对每个季度的平均气温进行简单相关分析；第二、 在控制年平均气温的情况下，再分别对各季度的平均气温进行偏相关分析。3.2.1各季度平均气温的相关分析图3相关性1月2月3月1月Pearson 相关性1.989*.981*显著性（双侧）.000.000N3131312月Pearson 相关性.989*1.984*显著性（双侧）.000.000N3131313月Pearson 相关性.981*.984*1显著性（双侧）.000.000N313131*. 在 .01 水平（双侧）上显著相关。图3显示的是1月、2月和3月的平均气温的相关性，可以看出这三个月的相关性比较明显。通过这三

6、个月的相关分析我们容易得到其他三个季度的气温也具有比较强的相关性。3.2.2各个季度平均气温的偏相关分析图4相关性控制变量1月2月3月年平均1月相关性1.000.913.724显著性（双侧）.000.000df028282月相关性.9131.000.859显著性（双侧）.000.000df280283月相关性.724.8591.000显著性（双侧）.000.000.df28280在控制年平均气温的情况下第一季度个月平均气温的相关性有所下降，但是没有改变它们的相关性。说明各个季度各个月份的平均气温都具有非常强的相关性。图5图5表示的是各个月份平均气温的相关图，可以看出每个月份之间的相关性都非常好

7、。3.3聚类分析图6图7图8聚类成员案例号城 市聚类距离1北京 34.6862天津 34.2153石家庄 36.4704太原 35.9965呼和浩特 14.7726沈阳 15.2917长春 17.2958哈尔滨 110.3969上海 23.93310南京 24.35411杭州 23.89612合肥 24.86913福州 48.10614南昌 25.97715济南 37.86016郑州 28.48717武汉 24.32818长沙(望城) 24.95319广州 41.88520南宁 43.34621海口 410.97322重庆(沙坪坝)27.87023成都(温江) 25.00524贵阳 211.2

8、2325昆明 216.41926拉萨 317.50427西安(泾河) 37.07928兰州(皋兰) 18.33829西宁 113.20830银川 39.75631乌鲁木齐 17.466通过聚类分析，由图6到图8我们可以清晰的看出，在图6和图8中对各个城市的划分比较统一，而在图7中我们显然发现对于气温我们可以分成4类，分别是5、6、7、8、9月为第一类，4、10月为第二类，3、11月第三类，1、2、12月为第四类，很显然是根据气温的从高到低进行分类的，第一类温度最高，第二类次之，第三类温度稍低，第四类温度最低。3.4因子分析在上面的研究中，我们分别对温度和城市进行了分类，并且对气温进行了相关分析

9、，得出了气温的分布情况，并且研究了气温对城市的影响，但是由于各种变量的多重共线性问题，不太利于具体问题意义解释，因此我们可以利用因子分析对各个月份的平均气温变量提取公因子，对数据进行降维处理。首先利用主菜单分析中的降维工具条进行因子分析。在对话框中，依次将“1月”、“2月”、“3月”、“4月”、“5月”、“6月”、“7月”、“8月”、“9月”、“10月”、“11月”和“12月” 添加到“变量”列表中。接着在相应的按钮中选中输出相关系数矩阵和KMO和Bartlett的球形度检验的复选框，并且设置特征根大于1为提取公因子的要求，利用最大方差法进行因子旋转，最后保存因子得分变量，并显示因子得分系数矩

10、阵，输出结果参见图9至图15（1）KMO 和 Bartlett 的检验，如图9所示，KMO的取值为0.811，表明变量间有较强的相关性，数据很适合做因子分析。Bartlett检验的Sig.值为0.000，说明数据来自正态总体，适合进一步分析。图9KMO 和 Bartlett 的检验取样足够度的 Kaiser-Meyer-Olkin 度量。.811Bartlett 的球形度检验近似卡方1009.087df66Sig.000（2）变量共同度，指的是按照所选标准提取相应数量主成分后，各变量中信息分别被提取的比例。如图10所示，所有变量共同度都在90%以上，所以提取这几个公因子对各变量的解释力非常好。

11、图10公因子方差初始提取1月1.000.9842月1.000.9873月1.000.9904月1.000.9525月1.000.8906月1.000.9257月1.000.9628月1.000.9249月1.000.95210月1.000.98411月1.000.98712月1.000.986提取方法：主成份分析。（3）解释的总方差，由图11中可以看出，“初始特征值”一栏显示只有前两个特征值大于1，所以只选取了前两个公因子；“提取平方和载入”一栏显示第一公因子的方差贡献率是82.367%，前两个公因子的方差总和占所有主成分方差的96.014%，可见选取前两个因子已足够替代原来的变量。图11解释

12、的总方差成份初始特征值提取平方和载入旋转平方和载入合计方差的 %累积 %合计方差的 %累积 %合计方差的 %累积 %19.88482.36782.3679.88482.36782.3676.45853.81353.81321.63813.64696.0141.63813.64696.0145.06442.20196.0143.2602.16698.1804.074.61898.7985.065.54599.3436.032.26699.6107.025.20899.8188.008.07099.8889.007.05699.94410.005.04099.98411.001.01099.994

13、12.001.006100.000提取方法：主成份分析。（4）碎石图，是按照特征根大小排列的主成分散点图。如图12所示，我们可以看到有；两个成分的特征值超过1。图12（5）成分矩阵，如表13可见，所有的因子都聚集在第一主成分中，因此有必要进行因子旋转。我们利用最大方差法进行因子旋转得到的旋转成分矩阵如图13所示。图13成份矩阵a成份1210月.992-.0324月.975-.00211月.958-.2629月.951.2163月.931-.3501月.927-.3525月.916.22712月.906-.4072月.895-.4328月.838.4716月.821.5007月.750.632提取方法 :主成份。a. 已提取了 2 个成份。图14旋转成份矩阵a成份122月.962.24612月.955.2733月.937.3331月.936.32911月.902.41710月.779.6154月.747.6277月.166.9666月.305.9128月.337.9009月.588.7785月.554.763提取方法 :主成份。 旋转法 :具有 Kaiser 标准化的正交旋转法。a. 旋转在 3 次迭代后收敛。这样每个因子都很明确了，第一主成分包括1、2、3、4、10、11、12月，可以命名为较冷气温；第二主成分包括5、6、7、8、9月，可以命名为叫人气温。