首先介绍一下主成分分析法:1 I4 M& l6 E( X1 O$ B
主成分分析是一种通过降维技术把多个指标约化为少数几个综合指标的综合统计分析方法,而这些综合指标能够反映原始指标的绝大部分信息,它们通常表现为原始几个指标的线性组合。( _+ R6 V& y4 c& Q% D! L3 ^
简单说就是 “你要将后宫佳丽三千来分个漂亮等级来赏赐,由于都很好看所以不知道怎么分,这时你练成了主成分分析大法,一下选出了20个不同等级的佳丽,其他2980个佳丽和这20个比就可以分开了。从此你后宫无忧”' ~+ s4 O/ _! q9 Z
[size=42.6667px] 它的基本思想是:具有一定相关性的众多指标重新组合成新的无相互关系的综合指标来代替。通常数学上的处理就是将这P个指标进行线性组合作为新的综合指标。问题是:这样的线性组合会很多,如何选择?
# ^, {. I' b7 T/ ?
* s0 ^9 I$ D( R1 @; Y# U$ i* @' H 如果选取的第一个线性组合即第一个综合指标为 F1, 希望它尽可能多的反映原来指标的信息,即方差(var(F1))越大,所包含的原来指标信息就越多,F1的方差应该是最大的,把它称作“第一主成分”,方差其次就是第二主成分...
/ Y& ^9 z; S' l6 k4 J, S$ l
* K9 q. d$ Y& f/ c S 但是我们不可能一直选到第100主成分,那样就是浪费劳动力,社会主义精神不允许我们这么做!" t1 Q/ G$ \! n3 A: X
因此,我们选到第几主成分就可以不选了?或者说,我们选几个主成分就可以代表原来那3000佳丽,难道一定得是20个?; K5 a R4 g' X1 w6 e
85%准则:% e) N% I9 U1 b: x1 B6 s
先准备一个概念:累计贡献率7 s* R; `( f3 d* b9 o5 ?* c
累积贡献率:前k个主成分共有多大的综合能力,用这k个主成分的方差和在全部方差中所占比重来描述,称为累积贡献率
- J. R; R+ T' q- y/ t4 U& B0 ?) f& _0 S7 e/ G. m% w
我们进行主成分分析的目的之一是希望用尽可能少的主成分F1,F2,…,Fk(k≤p)代替原来的P个指标。到底应该选择多少个主成分,在实际工作中,主成分个数的多少取决于能够反映原来变量85%以上的信息量为依据,即当累积贡献率≥85%时的主成分的个数就足够了。最常见的情况是主成分为2到3个。
9 x8 X- ^$ @$ G/ M7 _) R/ J9 w$ b; } `/ F9 Q1 A q7 T/ p5 @
但是2 Q! `4 P1 b0 X
你实际情况下面对的一群变量不会是一群乖孩子等着你去分类,它们可能会很狂躁,比如变量的单位不一致,变化方向不一致,这样就导致了各个变量分散程度加大,总体的方差就会受到方差较大的变量的控制。/ C2 d' D, o: \. B! D- }
因此嘞,我们要做一件事:标准化处理变量,学秦始皇一统天下
' W5 ]- | ^( Q! y! c 嘞原来变量我们记为X我们要做的就是将X标准化为X*:
$ _3 p3 k, D) k% O" B8 C8 i# ~
! z! e9 t5 S2 P* \% |; A! g( [9 @7 ~$ ~: u/ @
$ f4 ^3 ?$ q5 p/ O) K
0 d1 @5 b0 U+ R) v
c" Z( ]8 Z, R$ \( e& p
8 S, A/ Y. o) F- b( F% i% e1 X1 m4 S
/ l* W( K& n- D- Y7 r之后就可以进行我们的主成分分析操作了:' p. b L1 f0 W4 W. Y# S( |
STEP01:计算X*的相关系数矩阵R; STEP03:求相关系数矩阵R的特征根 λ1≥λ2 ≥‥‥≥ λp ≥0及相应的单位正交特征向量U1,U2,‥‥, Up; STEP04:计算方差累积贡献率,确定主成分的个数q; STEP05:写出主成分F=X*U,解释其实际经济意义并指导实践。 ( u# Y6 x4 U. s: L
在下一期我将会用实际例题来加深大家对主成分分析法的掌握,这期主要是让大家对这个有一个原理上的认识。
% a' `1 ~4 r) P1 x- B4 C5 g/ c' H" F
0 O: D- r# L' W) L3 l' x
" P# Y3 @0 P r下载附件了解主成分分析的详细内容!!!!快下载,别磨蹭
) ~) t- ^3 H9 T, v& @6 m
# z1 x, N' `6 G# V% z |