如何利用SPSS进行主成分分析.docx
- 文档编号:27292542
- 上传时间:2023-06-28
- 格式:DOCX
- 页数:26
- 大小:331.19KB
如何利用SPSS进行主成分分析.docx
《如何利用SPSS进行主成分分析.docx》由会员分享,可在线阅读,更多相关《如何利用SPSS进行主成分分析.docx(26页珍藏版)》请在冰豆网上搜索。
如何利用SPSS进行主成分分析
利用SPSS进行主成分分析
【例子】以全国31个省市的8项经济指标为例,进行主成分分析。
第一步:
录入或调入数据(图1)。
图1原始数据(未经标准化)
第二步:
打开“因子分析”对话框。
沿着主菜单的“Analyze→DataReduction→Factor
”的路径(图2)打开因子分析选项框(图3)。
图2打开因子分析对话框的路径
图3因子分析选项框
第三步:
选项设置。
首先,在源变量框中选中需要进行分析的变量,点击右边的箭头符号,将需要的变量调入变量(Variables)栏中(图3)。
在本例中,全部8个变量都要用上,故全部调入(图4)。
因无特殊需要,故不必理会“Value
”栏。
下面逐项设置。
图4将变量移到变量栏以后
设置Descriptives选项。
单击Descriptives按钮(图4),弹出Descriptives对话框(图5)。
图5描述选项框
在Statistics栏中选中Univariatedescriptives复选项,则输出结果中将会给出原始数据的抽样均值、方差和样本数目(这一栏结果可供检验参考);选中Initialsolution复选项,则会给出主成分载荷的公因子方差(这一栏数据分析时有用)。
在CorrelationMatrix栏中,选中Coefficients复选项,则会给出原始变量的相关系数矩阵(分析时可参考);选中Determinant复选项,则会给出相关系数矩阵的行列式,如果希望在Excel中对某些计算过程进行了解,可选此项,否则用途不大。
其它复选项一般不用,但在特殊情况下可以用到(本例不选)。
设置完成以后,单击Continue按钮完成设置(图5)。
设置Extraction选项。
打开Extraction对话框(图6)。
因子提取方法主要有7种,在Method栏中可以看到,系统默认的提取方法是主成分(),因此对此栏不作变动,就是认可了主成分分析方法。
在Analyze栏中,选中Correlationmatirx复选项,则因子分析基于数据的相关系数矩阵进行分析;如果选中Covariancematrix复选项,则因子分析基于数据的协方差矩阵进行分析。
对于主成分分析而言,由于数据标准化了,这两个结果没有分别,因此任选其一即可。
在Display栏中,选中Unrotatedfactorsolution(非旋转因子解)复选项,则在分析结果中给出未经旋转的因子提取结果。
对于主成分分析而言,这一项选择与否都一样;对于旋转因子分析,选择此项,可将旋转前后的结果同时给出,以便对比。
选中ScreePlot(“山麓”图),则在分析结果中给出特征根按大小分布的折线图(形如山麓截面,故得名),以便我们直观地判定因子的提取数量是否准确。
在Extract栏中,有两种方法可以决定提取主成分(因子)的数目。
一是根据特征根(Eigenvalues)的数值,系统默认的是
。
我们知道,在主成分分析中,主成分得分的方差就是对应的特征根数值。
如果默认
,则所有方差大于等于1的主成分将被保留,其余舍弃。
如果觉得最后选取的主成分数量不足,可以将
值降低,例如取
;如果认为最后的提取的主成分数量偏多,则可以提高
值,例如取
。
主成分数目是否合适,要在进行一轮分析以后才能肯定。
因此,特征根数值的设定,要在反复试验以后才能决定。
一般而言,在初次分析时,最好降低特征根的临界值(如取
),这样提取的主成分将会偏多,根据初次分析的结果,在第二轮分析过程中可以调整特征根的大小。
第二种方法是直接指定主成分的数目即因子数目,这要选中Numberoffactors复选项。
主成分的数目选多少合适?
开始我们并不十分清楚。
因此,首次不妨将数值设大一些,但不能超过变量数目。
本例有8个变量,因此,最大的主成分提取数目为8,不得超过此数。
在我们第一轮分析中,采用系统默认的方法提取主成分。
图6提取对话框
需要注意的是:
主成分计算是利用迭代(Iterations)方法,系统默认的迭代次数是25次。
但是,当数据量较大时,25次迭代是不够的,需要改为50次、100次乃至更多。
对于本例而言,变量较少,25次迭代足够,故无需改动。
设置完成以后,单击Continue按钮完成设置(图6)。
设置Scores设置。
选中Saveasvariables栏,则分析结果中给出标准化的主成分得分(在数据表的后面)。
至于方法复选项,对主成分分析而言,三种方法没有分别,采用系统默认的“回归”(Regression)法即可。
图7因子得分对话框
选中Displayfactorscorecoefficientmatrix,则在分析结果中给出因子得分系数矩阵及其相关矩阵。
设置完成以后,单击Continue按钮完成设置(图7)。
其它。
对于主成分分析而言,旋转项(Rotation)可以不必设置;对于数据没有缺失的情况下,Option项可以不必理会。
全部设置完成以后,点击OK确定,SPSS很快给出计算结果(图8)。
图8主成分分析的结果
第四步,结果解读。
在因子分析结果(Output)中,首先给出的DescriptiveStatistics,第一列Mean对应的变量的算术平均值,计算公式为
第二列Std.Deviation对应的是样本标准差,计算公式为
第三列AnalysisN对应是样本数目。
这一组数据在分析过程中可作参考。
接下来是CorrelationMatrix(相关系数矩阵),一般而言,相关系数高的变量,大多会进入同一个主成分,但不尽然,除了相关系数外,决定变量在主成分中分布地位的因素还有数据的结构。
相关系数矩阵对主成分分析具有参考价值,毕竟主成分分析是从计算相关系数矩阵的特征根开始的。
相关系数阵下面的Determinant=1.133E-0.4是相关矩阵的行列式值,根据关系式
可知,det(λI)=det(R),从而Determinant=1.133E-0.4=λ1*λ2*λ3*λ4*λ5*λ6*λ7*λ8。
这一点在后面将会得到验证。
在Communalities(公因子方差)中,给出了因子载荷阵的初始公因子方差(Initial)和提取公因子方差(Extraction),后面将会看到它们的含义。
在TotalVarianceExplained(全部解释方差)表的InitialEigenvalues(初始特征根)中,给出了按顺序排列的主成分得分的方差(Total),在数值上等于相关系数矩阵的各个特征根λ,因此可以直接根据特征根计算每一个主成分的方差百分比(%ofVariance)。
由于全部特征根的总和等于变量数目,即有m=∑λi=8,故第一个特征根的方差百分比为λ1/m=3.755/8=46.939,第二个特征根的百分比为λ2/m=2.197/8=27.459,……,其余依此类推。
然后可以算出方差累计值(Cumulative%)。
在ExtractionSumsofSquaredLoadings,给出了从左边栏目中提取的三个主成分及有关参数,提取的原则是满足λ>1,这一点我们在图6所示的对话框中进行了限定。
图8特征根数值衰减折线图(山麓图)
主成分的数目可以根据相关系数矩阵的特征根来判定,如前所说,相关系数矩阵的特征根刚好等于主成分的方差,而方差是变量数据蕴涵信息的重要判据之一。
根据λ值决定主成分数目的准则有三:
只取λ>1的特征根对应的主成分
从TotalVarianceExplained表中可见,第一、第二和第三个主成分对应的λ值都大于1,这意味着这三个主成分得分的方差都大于1。
本例正是根据这条准则提取主成分的。
累计百分比达到80%~85%以上的λ值对应的主成分
在TotalVarianceExplained表可以看出,前三个主成分对应的λ值累计百分比达到89.584%,这暗示只要选取三个主成分,信息量就够了。
根据特征根变化的突变点决定主成分的数量
从特征根分布的折线图(ScreePlot)上可以看到,第4个λ值是一个明显的折点,这暗示选取的主成分数目应有p≤4(图8)。
那么,究竟是3个还是4个呢?
根据前面两条准则,选3个大致合适(但小有问题)。
在ComponentMatrix(成分矩阵)中,给出了主成分载荷矩阵,每一列载荷值都显示了各个变量与有关主成分的相关系数。
以第一列为例,0.885实际上是国内生产总值(GDP)与第一个主成分的相关系数。
将标准化的GDP数据与第一主成分得分进行回归,决定系数R2=0.783(图9),容易算出R=0.885,这正是GDP在第一个主成分上的载荷。
下面将主成分载荷矩阵拷贝到Excel上面作进一步的处理:
计算公因子方差和方差贡献。
首先求行平方和,例如,第一行的平方和为
h12=0.88492+0.38362+0.12092=0.9449
这是公因子方差。
然后求列平方和,例如,第一列的平方和为
s12=0.88492+0.60672+…+0.82272=3.7551
这便是方差贡献(图10)。
在Excel中有一个计算平方和的命令sumsq,可以方便地算出一组数据的平方和。
显然,列平方和即方差贡献。
事实上,有如下关系成立:
相关系数矩阵的特征根=方差贡献=主成分得分的方差
至于行平方和,显然与前面公因子方差(Communalities)表中的Extraction列对应的数据一样。
如果我们将8个主成分全部提取,则主成分载荷的行平方和都等于1(图11),即有hi=1,sj=λj。
到此可以明白:
在Communalities中,Initial对应的是初始公因子方差,实际上是全部主成分的公因子方差;Extraction对应的是提取的主成分的公因子方差,我们提取了3个主成分,故计算公因子方差时只考虑3个主成分。
图9国内生产总值(GDP)的与第一主成分的相关关系(标准化数据)
图10主成分方差与方差贡献
图11全部主成分的公因子方差和方差贡献
提取主成分的原则上要求公因子方差的各个数值尽可能接近,亦即要求它们的方差极小,当公因子方差完全相等时,它们的方差为0,这就达到完美状态。
实际应用中,只要公因子方差数值彼此接近(不相差太远)就行了。
从上面给出的结果可以看出:
提取3个主成分的时候,居民消费的公因子方差偏小,这暗示提取3个主成分,居民消费方面的信息可能有较多的损失。
至于方差贡献,反映对应主成分的重要程度,这一点从方差的统计学意义可以得到理解。
在图11中,将最后一行的特征根全部乘到一起,得0.0001133,这正是相关系数矩阵的行列式数值(在Excel中,求一组数据的乘积之和的命令是product)。
最后说明ComponentScoreCoefficientMatrix(成分得分系数矩阵)和ComponentScoreCovarianceMatrix(成分得分协方差矩阵),前者是主成分得分系数,后者是主成分得分的协方差即相关系数。
从ComponentScoreCovarianceMatrix可以看出,标准化主成分得分之间的协方差即相关系数为0(j≠k)或1(j=k),这意味着主成分之间彼此正交即垂直。
初学者常将ComponentScoreCoefficientMatrix表中的数据当成主成分得分或因子得分,这是误会。
成分得分系数矩阵的数值是主成分载荷除以相应的特征根得到的结果。
在ComponentMatrix表中,将第一列数据分别除以λ1=3.755,第二列数值分别除以λ2=2.197,…,立即得到ComponentScoreCoefficient;反过来,如果将ComponentScoreCoefficientMatrix表中的各列数据分别乘以λ1=3.755,λ2=2.197,…,则可将其还原为主成分载荷即ComponentMatrix中的数据。
实际上,主成分得分在原始数据所在的SPSS当前数据栏中给出,不过给出的都是标准化的主成分得分(图12a);将各个主成分乘以相应的√λ即特征根的二次方根可以将其还原为未经标准化的主成分得分。
a.标准化的主成分得分b.非标准化的主成分得分
图12两种主成分得分
计算标准化主成分得分的协方差或相关系数,结果与ComponentScoreCovarianceMatrix表中的给出的结果一致(见图13)。
第一因子
第二因子
第三因子
第一因子
1
第二因子
0.00000
1
第三因子
0.00000
0.00000
1
图13主成分(得分)之间的相关系数矩阵
第五步,计算结果分析。
从ComponentMatrix即主成分载荷表中可以看出,国内生产总值、固定资产投资和工业产值在第一主成分上载荷较大,亦即与第一主成分的相关系数较高;职工工资和货物周转量在第二主成分上的载荷绝对值较大,即负相关程度较高;消费价格指数在第三主成分上的载荷较大,即相关程度较高。
因此可将主成分命名如下:
第一主成分:
投入-产出主成分;
第二主成分:
工资-物流主成分;
第三主成分:
消费价格主成分。
问题在于:
一方面,居民消费和商品零售价格指数的归类比较含混;另一方面,主成分的命名结构不清。
因此,有必要作进一步的因子分析。
至于因子旋转之类,留待“因子分析”部分说明;计算结果的系统分析不属于软件操作范围,预备课堂讲解。
【说明】本人计算机是双系统,现在常用的WinMe系统出了毛病,SPSS10.0在WinMe系统中;故这次改用本人Win2000系统中的SPSS11.0。
对于因子分析之类,SPSS11.0与SPSS10.0基本没有什么差别。
表2DescriptiveStatistics
Mean
Std.Deviation
AnalysisN
水资源总
60.1120
15.30739
10
水资源开
51.8570
18.97639
10
总人口(
584.7250
6.62343
10
城镇化率
39.5950
9.15229
10
农业用水
13.7520
1.90408
10
工业用水
13.7260
3.11251
10
生活用水
2.8020
.50758
10
人均用水
525.50
68.430
10
全市生产
1149.9290
698.85881
10
固定资产
5372749.50
5046080.033
10
农田灌溉
422.30
79.452
10
万元工业
324.90
123.457
10
生态需水
.009700
.0056046
10
用户废污
53225.30
4569.898
10
表2相关系数矩阵CorrelationMatrix(a,b)
Statistics
Variables
Variables
水资源总
水资源开
总人口(
城镇化率
农业用水
工业用水
生活用水
人均用水
全市生产
固定资产
农田灌溉
万元工业
生态需水
用户废污
Correlation
水资源总
1.000
-.826
-.498
-.715
-.399
-.315
-.179
-.567
-.594
-.617
-.033
.544
-.721
.281
水资源开
-.826
1.000
.768
.820
.245
.634
.406
.813
.857
.867
-.291
-.819
.897
-.181
总人口(
-.498
.768
1.000
.912
-.223
.854
.472
.743
.946
.930
-.690
-.986
.938
-.173
城镇化率
-.715
.820
.912
1.000
.018
.814
.451
.834
.935
.930
-.464
-.894
.955
-.209
农业用水
-.399
.245
-.223
.018
1.000
-.234
-.459
.325
-.021
-.048
.773
.155
-.090
-.412
工业用水
-.315
.634
.854
.814
-.234
1.000
.676
.825
.821
.810
-.667
-.788
.788
.159
生活用水
-.179
.406
.472
.451
-.459
.676
1.000
.435
.498
.554
-.576
-.426
.544
.646
人均用水
-.567
.813
.743
.834
.325
.825
.435
1.000
.858
.839
-.221
-.740
.775
-.135
全市生产
-.594
.857
.946
.935
-.021
.821
.498
.858
1.000
.994
-.499
-.960
.953
-.236
固定资产
-.617
.867
.930
.930
-.048
.810
.554
.839
.994
1.000
-.499
-.945
.962
-.185
农田灌溉
-.033
-.291
-.690
-.464
.773
-.667
-.576
-.221
-.499
-.499
1.000
.630
-.571
-.320
万元工业
.544
-.819
-.986
-.894
.155
-.788
-.426
-.740
-.960
-.945
.630
1.000
-.945
.254
生态需水
-.721
.897
.938
.955
-.090
.788
.544
.775
.953
.962
-.571
-.945
1.000
-.129
用户废污
.281
-.181
-.173
-.209
-.412
.159
.646
-.135
-.236
-.185
-.320
.254
-.129
1.000
TableCaption
aDeterminant=.000
bThismatrixisnotpositivedefinite.
表3公因子方差Communalities
Initial
Extraction
水资源总
1.000
.742
水资源开
1.000
.899
总人口(
1.000
.983
城镇化率
1.000
.932
农业用水
1.000
.935
工业用水
1.000
.850
生活用水
1.000
.927
人均用水
1.000
.857
全市生产
1.000
.971
固定资产
1.000
.961
农田灌溉
1.000
.933
万元工业
1.000
.978
生态需水
1.000
.968
用户废污
1.000
.957
ExtractionMethod:
PrincipalComponentAnalysis.
表4特征值﹑主成分贡献率和累计贡献率TotalVarianceExplained
Component
InitialEigenvalues
ExtractionSumsofSquaredLoadings
RotationSumsofSquaredLoadings
Total
%ofVariance
Cumulative%
Total
%ofVariance
Cumulative%
Total
%ofVariance
Cumulative%
1
8.988
64.201
64.201
8.988
64.201
64.201
8.767
62.625
62.625
2
2.769
19.779
83.980
2.769
19.779
83.980
2.342
16.729
79.353
3
1.136
8.116
92.096
1.136
8.116
92.096
1.784
12.742
92.096
4
.591
4.222
96.318
5
.227
1.620
97.938
6
.180
1.288
99.225
7
.065
.467
99.692
8
.031
.223
99.915
9
.012
.085
100.000
10
2.535E-16
1.811E-15
100.000
11
5.751E-17
4.108E-16
100.000
12
3.109E-17
2.221E-16
100.000
13
-1.376E-16
-9.829E-16
100.000
14
-2.310E-16
-1.650E-15
100.000
ExtractionMethod:
PrincipalComponentAnalysis.
表6载荷矩阵ComponentMatrix(a)
Component
Z1
Z2
Z3
水资源总
-.644
-.520
-.237
水资源开
.883
.295
.180
总人口(
.958
-.090
-.239
城镇化率
.957
.126
-.023
农业用水
-.077
.895
.359
工业用水
.871
-.285
.097
生活用水
.574
-.608
.477
人均用水
.852
.248
.263
全市生产
.978
.082
-.081
固定资产
.978
.053
-.030
农田灌溉
-.576
.717
.295
万元工业
-.955
.000
.256
生态需水
.983
.026
-.024
用户废污
-.110
-.742
.628
ExtractionMethod:
PrincipalComponentAnalysis.
a3componentsextracted.
成分得分协方差矩阵RotatedComponent
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 如何 利用 SPSS 进行 成分 分析