主成分分析法实例
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1、主成分法:
用主成分法寻找公共因子的方法如下:
假定从相关阵出发求解主成分,设有p 个变量,则可找出p 个主成分。将所得的p 个主成分按由大到小的顺序排列,记为1Y ,2Y ,…,P Y , 则主成分与原始变量之间存在如下关系:
11111221221122221122....................p p p p p
p p pp p Y X X X Y X X X Y X X X
γγγγγγγγγ=+++⎧⎪
=+++⎪⎨⎪
⎪=+++⎩ 式中,ij γ为随机向量X 的相关矩阵的特征值所对应的特征向量的分量,因为特征向量之间彼此正交,从X 到Y 得转换关系是可逆的,很容易得出由Y 到
X 得转换关系为:
11112121212122221122....................p p p p p
p p pp p X Y Y Y X Y Y Y X Y Y Y
γγγγγγγγγ=+++⎧⎪
=+++⎪⎨⎪
⎪=+++⎩ 对上面每一等式只保留钱m 个主成分而把后面的部分用i ε代替,则上式变为:
111121211
2121222221122.................
...m m m m p p p mp m p X Y Y Y X Y Y Y X Y Y Y γγγεγγγεγγγε=++++⎧⎪=++++⎪⎨⎪⎪=++++⎩
上式在形式上已经与因子模型相一致,且i Y (i=1,2,…,m )之间相互独立,且i Y 与i ε之间相互独立,为了把i Y 转化成合适的公因子,现在要做的工作只是把主成分i Y 变为方差为1的变量。为完成此变换,必须将i Y 除以其标准差,由主成分分析的知识知其标准差即为特征根的平方根
i λ/i i i F Y λ=,
1122m m λγλγλγ,则式子变为:
1111122112211222221122....................m m m m p p p pm m p X a F a F a F X a F a F a F X a F a F a F εεε=++++⎧⎪=++++⎪⎨⎪
⎪=++++⎩
这与因子模型完全一致,这样,就得到了载荷A 矩阵和 初始公因子(未旋转)。
一般设A ∧
为样本相关矩阵R 的特征根,12,,...,p γγγ为对应的标准正交化特征向量。设m
A ∧
=(1122,,...,m m λγλγλγ)
共同度的估计为:
2
2
2
2
12...i i i im h a a a ∧∧∧∧=+++
下面用主成分法分析以下数据:
步骤:
第一步,把Excel 中的数据导入到SPSS 中:File →Open →Data ; 第二步,数据标准化:Analyze →Descriptive Statistics →Descriptives 如图:
第三步,检验数据:
如图:
得到结果如下:
KMO 和Bartlett 的检验
取样足够度的Kaiser-Meyer-Olkin 度量。.754 Bartlett 的球形度检验近似卡方
df 36
Sig. .000
Sig小于,所以该数据可用;
第四步,用主成分法分析数据:Analyze→Dimension Reduction→Factor 如图:
得到结果如下图:
相关矩阵
Zscore: 100固定资产原值实现值(%)
Zscore:
100元固
定资产原
值实现利
税(%)
Zscore:
100元资
金实现利
税(%)
Zscore:
100元工
业总产值
实现利税
(%)
Zscore:
100元销
售收入实
现利税
(%)
Zscore(
每吨标准
煤实现工
业产值
(元))
Zscore(
每千瓦时
电力实现
工业产值
(元))
Zscore:
全员劳动
生产率
(元/人.
年)
Zscore:
100元流
动资金实
现产值
(元)
相关Zscore: 100固
定资产原值实现
值(%)
.869 .770 .211 .920 .899 .795 .896
Zscore: 100元
固定资产原值实
现利税(%)
.869 .978 .387 .472 .886 .804 .814 .849
Zscore: 100元
资金实现利税
(%)
.770 .978 .523 .531 .797 .736 .740 .811
Zscore: 100元
工业总产值实现
利税(%)
.387 .523 .323 .115 .125 .051
Zscore: 100元
销售收入实现利
税(%)
.211 .472 .531 .323 .175 .260 .371 .317
Zscore(每吨标
准煤实现工业产
值(元))
.920 .886 .797 .115 .175 .877 .815 .768
Zscore(每千瓦
时电力实现工业
产值(元))
.899 .804 .736 .260 .877 .757 .818
Zscore: 全员
劳动生产率(元/
人.年)
.795 .814 .740 .125 .371 .815 .757 .715
Zscore: 100元
流动资金实现产
值(元)
.896 .849 .811 .051 .317 .768 .818 .715