利用Matlab和SPSS软件进行主成分分析.docx

1、利用Matlab和SPSS软件进行主成分分析9. 利用Matlab和SPSS实现主成分分析1.直接调用Matlab软件实现在软件Matlab中实现主成分分析可以采取两种方式实现：一是通过编程来实现；二是直接调用Matlab中自带程序实现。通过直接调用Matlab中的程序可以实现主成分分析：式中：X为输入数据矩阵 （一般要求nm）输出变量：pc 主分量fi的系数，也叫因子系数；注意：pcTpc=单位阵score是主分量下的得分值；得分矩阵与数据矩阵X的阶数是一致的；variance是score对应列的方差向量，即A的特征值；容易计算方差所占的百分比 percent-v = 100*varianc

2、e/sum(variance);t2表示检验的t2-统计量（方差分析要用）计算过程中应用到计算模型: (要求pm)例：表1为某地区农业生态经济系统各区域单元相关指标数据，运用主成分分析方法可以用更少的指标信息较为精确地描述该地区农业生态经济的发展状况。表1 某农业生态经济系统各区域单元的有关数据样本序号x1：人口密度(人/km2)x 2：人均耕地面积(ha)x 3：森林覆盖率(%)x 4：农民人均纯收入(元/人)x 5：人均粮食产量 (kg/人)x 6：经济作物占农作物播面比例()x 7：耕地占土地面积比率()x 8：果园与林地面积之比()x 9：灌溉田占耕地面积之比()1363.9120.3

3、5216.101192.11295.3426.72418.4922.23126.2622141.5031.68424.3011 752.35452.2632.31414.4641.45527.0663100.6951.06765.6011 181.54270.1218.2660.1627.47412.4894143.7391.33633.2051 436.12354.2617.48611.8051.89217.5345131.4121.62316.6071 405.09586.5940.68314.4010.30322.932668.3372.03276.2041 540.29216.398.

4、1284.0650.0114.861795.4160.80171.106926.35291.528.1354.0630.0124.862862.9011.65273.3071 501.24225.2518.3522.6450.0343.201986.6240.84168.904897.36196.3716.8615.1760.0556.1671091.3940.81266.502911.24226.5118.2795.6430.0764.4771176.9120.85850.302103.52217.0919.7934.8810.0016.1651251.2741.04164.609968.3

5、3181.384.0054.0660.0155.4021368.8310.83662.804957.14194.049.1104.4840.0025.7901477.3010.62360.102824.37188.0919.4095.7215.0558.4131576.9481.02268.0011 255.42211.5511.1023.1330.0103.4251699.2650.65460.7021 251.03220.914.3834.6150.0115.59317118.5050.66163.3041 246.47242.1610.7066.0530.1548.70118141.47

6、30.73754.206814.21193.4611.4196.4420.01212.94519137.7610.59855.9011 124.05228.449.5217.8810.06912.65420117.6121.24554.503805.67175.2318.1065.7890.0488.46121122.7810.73149.1021 313.11236.2926.7247.1620.09210.078对于上述例子，Matlab进行主成分分析，可以得到如下结果。1 以及每一个主成分的贡献率和累计贡献率，如表2和图1。表2. 特征根及主成分贡献率 主成分 特征值 贡献率 累积贡献率

7、 14.66151.79151.791 22.08923.21675.007 31.04311.58986.596 40.5075.63892.234 50.3153.50295.736 60.1932.14097.876 70.1141.27199.147 84.533E-020.50499.650 93.147E-020.350100.000 图1 特征根2 前3几个主成分的载荷系数如表3所示。表3 前三个主成分在原变量上的载荷 前三个主成分 变量123 X10.158-0.255-0.059 X20.0260.424-0.027 X3-0.2070.0460.091 X40.0090.4

8、150.036 X50.1740.212-0.011 X60.1760.0860.120 X70.200-0.064-0.241 X80.042-0.0480.930 X90.207-0.0120.088 2. 直接调用SPSS软件实现多元分析处理的是多指标的问题。由于指标太多，使得分析的复杂性增加。观察指标的增加本来是为了使研究过程趋于完整，但反过来说，为使研究结果清晰明了而一味增加观察指标又让人陷入混乱不清。由于在实际工作中，指标间经常具备一定的相关性，故人们希望用较少的指标代替原来较多的指标，但依然能反映原有的全部信息，于是就产生了主成分分析、对应分析、典型相关分析和因子分析等方法。调用

9、Data Reduction菜单的Factor过程命令项，可对多指标或多因素资料进行因子分析。因子分析的基本目的就是用少数几个因子去描述许多指标或因素之间的联系，即将相关比较密切的几个变量归在同一类中，每一类变量就成为一个因子（之所以称其为因子，是因为它是不可观测的，即不是具体的变量，这与上一章的聚类分析不同），以较少的几个因子反映原资料的大部分信息。基本操作以以上例子来说明在SPSS中进行因子分析的整个过程。将以上数据导入到数据窗口中，先定义各变量为numberic型。激活Analysis菜单选Data Reduction的Factor.命令项，弹出Factor Analysis对话框（图1

10、）。在对话框左侧的变量列表中选变量X1至X9，点击钮使之进入Variables框。图 1 因子分析对话框点击Descriptives.钮，弹出Factor Analysis:Descriptives对话框（图2），在Statistics中选Univariate descriptives项要求输出各变量的均数与标准差，在Correlation Matrix栏内选Coefficients项要求计算相关系数矩阵，并选KMO and Bartletts test of sphericity项，要求对相关系数矩阵进行统计学检验。点击Continue钮返回Factor Analysis对话框。图2 描述性

11、指标选择对话框点击Extraction.钮，弹出Factor Analysis:Extraction对话框（图3），系统提供如下因子提取方法：图3 因子提取方法选择对话框Principal components：主成分分析法； Unweighted least squares：未加权最小平方法； Generalized least squares：综合最小平方法； Maximum likelihood：极大似然估计法； Principal axis factoring：主轴因子法； Alpha factoring：因子法； Image factoring：多元回归法。本例选用Principal

12、 components方法，之后点击Continue钮返回Factor Analysis对话框。点击Rotation.钮，弹出Factor Analysis:Rotation对话框，系统有5种因子旋转方法可选：None：不作因子旋转； Varimax：正交旋转； Equamax：全体旋转，对变量和因子均作旋转； Quartimax：四分旋转，对变量作旋转； Direct Oblimin：斜交旋转。 旋转的目的是为了获得简单结构，以帮助我们解释因子。本例选正交旋转法，之后点击Continue钮返回Factor Analysis对话框。点击Scores.钮，弹出弹出Factor Analysis:

13、Scores对话框，系统提供3种估计因子得分系数的方法，本例选Regression（回归因子得分），之后点击Continue钮返回Factor Analysis对话框，再点击OK钮即完成分析。结果解释在输出结果窗口中将看到如下统计数据：系统首先输出各变量的均数（Mean）与标准差（Std Dev），并显示共有21例观察单位进入分析；接着输出相关系数矩阵（Correlation Matrix），经Bartlett检验表明：Bartlett值 = 159.767，P0.0001，即相关矩阵不是一个单位矩阵，故考虑进行因子分析。Correlation Matrix人口密度人均耕地面积森林覆盖率农民人均纯收入人均粮食产量经济作物占农作物播面比例耕地占土地面积比率果园与林地面积之比灌溉田占耕地面积之比人口密度1.000-.327-.714-.336.309.408.790.156.744人均耕地面积-.3271.000-.035.644.420.255.009-.078.094森林覆盖率-.714-.0351.000.070-.740-.755-.930-.109-.