数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
查看: 6618|回复: 2

请教陆老师,作PCA时,主成分系数矩阵列向量的符号是否还有特殊的确定方法?

[复制链接]
发表于 2017-4-6 02:21 | 显示全部楼层 |阅读模式
本帖最后由 Ysu2008 于 2017-4-5 18:36 编辑

多元统计书上的一个例子,对20个样品,4个指标的数据作主成分分析。

4个指标相关系数矩阵 R(与书上的输出相同)=
    1.0000    0.6950    0.2195    0.0249
    0.6950    1.0000   -0.1480    0.1351
    0.2195   -0.1480    1.0000    0.0713
    0.0249    0.1351    0.0713    1.0000

R 的特征值(与书上相同) = [1.7183    1.0935    0.9813    0.2069]

然后根据 R 和特征值构造4个基本方程组,依次求解并作斯密特正交,得到主成分系数矩阵:
    0.7000    0.0950   -0.2400    0.6659
    0.6898   -0.2836    0.0585   -0.6636
    0.0879    0.9042   -0.2703   -0.3189
    0.1628    0.3050    0.9305    0.1208
但书上的输出是:
    0.7000    0.0950   -0.2400   -0.6659
    0.6898   -0.2836    0.0585    0.6636
    0.0879    0.9042   -0.2703    0.3189
    0.1628    0.3050    0.9305   -0.1208
——第4个列向量符号相反。

该书是用 SAS 作的计算,给出了计算公式以及每一步的计算结果。我按照该书的方法用 MATLAB 编的程序,在解4个齐次方程组时,只计算了基础解系,没作其他处理。我还将这个例子用 SPSS 复算了一次,SPSS的输出与书上的一致,难道这些商业软件有未公开的秘密算法?
发表于 2017-4-6 06:05 | 显示全部楼层
计算出来的特征向量,只要互相正交,本身取正号或负号,没有关系,都是正确的。

点评

哦,多谢陆老师。  发表于 2017-4-6 19:32
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2026-5-16 18:54 , Processed in 0.309533 second(s), 16 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表