因子分析方法

小丫头妮妮呵呵

第十章 因子分析法  

☆　引 言

● 因子分析是研究相关矩阵的内部依赖关系，它将多个变量综合为少数几个“因子”，以再现原始变量与“因子”之间的相关关系。在统计学领域内，它是属于多元分析的范畴。本章包括：
$ T' w2 \) W$ W( |9 o7 A10.1 通俗概述  10.2 方法原理
4 C/ d4 i5 w& V3 s10.3 因子载荷的统计意义  10.4 方差最大正交旋转
10.5 Promax斜旋转 10.6 因子得分  10.7 应用实例
[img=584,5]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
10.1 通俗概述
% }5 Q) t- ?' q( S! _0 `. L6 c+ b0 {    因子分析主要是由心理学家发展起来的， 1904年Chales Spearman提出这种方法用来解决智力测验得分的统计分析。目前因子分析在心理学、社会学、经济学、人口学、地质学、生理学，甚至在化学和物理学中都取得了功的应用，它的主要应用有二个方面，一是寻求基本结构、简化观测系统。给定一组变量或观测数据，我们要问，变量的维数是否一定需要这么多，是否存在一个子集，特别是一个加权子集，来解释整个问题。通常采用因子分析的方法将为数众多的变量减少为几个新因子，以再现它们之间的内在联系，二是用于分类，将变量或者样本进行分类，根据因子得分值，在因子轴所构成的空间中进行分类处理。
+ ?* [( b' g4 \1 B9 j    因子分析与主成分分析有很大的区别，主成分分析不能作为一个模型来描述，它只能作为通常的变量变换，而因子分析需要构造因子模型；主成分分析的主分量数M和变量数P相等，它是将一组具有相关性的变量变换为一组独立的变量，而因子分析的目的是要使M比P小，而且要尽可能地选取小的M，以便尽可能地构造一个结构简单的模型，主成分分析是将主分量表示为原观测变量的线住组合，而因子分析是将原观测变量表示为新因子的线性组合，即为新因子的综合指标。所以因子分析是主成分分析的发展。
[img=584,5]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]

10.2 方法原理
+ U- \# J2 g& ?7 v6 x, |     因子分析的结果不仅是要给出因子模型，而且还要得出变量和因子间的相关系数，这些相关系数构成“因子结构”。
    一个完全的因子解包括因子模型和因子结构两个方面，因子结构反映变量与因子间的相关关系，而因子模型则是以回归方程的形式将变量表为因子的线性组合。
    因子分析的基本问题是用变量之间的相关系数未决定因子载荷。
    因子模型的求解过程如下：
9 M' p3 |2 t8 N' t    设原始数据矩阵为：
[img=280,124]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
7 ?& G* B2 Y8 a% np表示样本数，n表示变量数。
; Z# N' b: e+ Q+ o  {) d7 U& F; O将原始数据进行标准化变换。
    [img=288,83]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
! l. y! `1 E! E" @其中：
    [img=400,128]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
经标准化变换后xij'的均值为0，方差为1，这样相关矩阵R和协方差矩阵S完全一样。这时相关矩阵：
    R=X·X'        (10.2.2)
% D: [. v' [% K/ @0 {1 Z（为方便计，假定经标准化处理后的矩阵仍记为X）。
; k0 U( ^6 A, ?8 U: O5 `$ t4 n求解R矩阵的特征方程组|R-λI|=0，记特征值为λ1>λ2>...>λp≥0，特征向量矩阵为
    [img=264,130]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
6 ^9 V+ C3 Y$ y1 z4 w9 K: K这样有关系：
    [img=256,128]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
( a# }, z" a, C) vU为正交矩阵，并且满足：
: M8 |! Q4 \$ f3 V1 _6 q8 _( k+ NU'U=UU'=I
即有：
- I( Z7 Q: N* w4 O# o$ J1 ^    [img=288,128]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]   （10.2.3)
将（10.2.3)式两边左乘以U’，右乘以U得：
' ^. q! z& j* Z# {   [img=400,256]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]     [img=80,77]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
1 {( Q/ L! H# M7 C( i令F=U'X
于是(10.2.4)式变为
[img=240,127]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
! Z+ d) P+ r9 YF为主因子阵，并且Fα＝U'X。(α=1，2,...n)，即每一个Fα为第α个样品主因子观测值。
    在因子分析中，通常只选其中m个（m＜p)主因子。就是根据变量的相关选出第一主因子F1，使其在各变量的公共因子方差中所占的方差贡献为最大，然后消去这个因子的影响，而从剩余的相关中，选出与F1不相关的因子F2，使其在各个变量的剩余因子方差贡献中为最大，……这样直到各个变量公共因子方差被分解完毕为止。
  ?: m$ a# t: q, @! W2 \. f    例如我们按所选取的主因子的信息量的和占总体信息量的85％，即有:
* J- _1 E7 v) P这样的M，使得：
! y) J2 ], B& K  o% t4 N( [   [img=244,80]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
这们个主因子将U矩阵剖分为：
: J' q1 u% i; X   [img=284,68]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
) R" V) V6 B( D. u0 ]2 c5 I由F=U'X将此式两端左乘U得之：
    [img=312,94]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
其中U(1)F(1)为m个主因子所能解释的部分，而U(2)F(2)为其残余部分。
. J" r# r+ v: N- O5 D  e# }记残余部分为：
  [img=216,122]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
+ w" b) b1 d6 [; Q则有:
   X=U(1)F(1)+ε             （10.2.5）
（10·2·5）式称为因子模型，U(1)称为因子载荷矩阵，F(1)称为主要因子，ε为特殊因子。
    因子模型的表达式即为（略去特殊因子）：
" Y# T! }5 r4 `& W  O5 a' Q5 w      [img=328,106]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
特征向量ui通常用单位向量表示，须进行规格化处理，即
    [img=132,24]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
( N* `) y8 J  {% o所以，因子载荷矩阵为。
$ `9 j, Y8 A7 W+ e: t    [img=480,129]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
# l: u1 i- C. c1 Y% a2 p9 j因此，R型的因子模型即为。
    [img=360,108]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
: l; @" k+ W& b$ A0 U/ v在因子模型中：
（1）F1、F2、...Fm叫做公共因子，它是在各个变量中共同出现的因子，它为高维空间中所张起的互相垂直的M个坐标轴。
- u+ X/ b: X8 Q4 k& H; Z（2）aij叫因子载荷，它是第i个变量在第j个主因子上的负荷，或者叫第i个变量在第j个主因子上的权，它反映了第i个变量在第j个主因子上的相对重要性。如果把xi看成m维因子空间上的一个向量，则aij表示xi在坐标轴Fj上的投影。
（3）εi为特殊因子，它相互独立地遵从正态分布N（0,σ2)，ai为特殊固子的载荷。
    在主成分分析中，公共因子数等于变量数，没有特殊因子，而在因于分析中，公共因子数少于变量数，有各变量的特殊因子，并且因子模型是将原始变量表为新因子的线性组合。
例10-1研究影响火柴销售的变量之间的关系。火柴是人民生活的必需品，商品虽小，但与日常生活关系十分密切，该商品易燃、易潮，较难储存，经营面广、量大，价格低，易受其他日用市场冲击的影响，供应脱节并影响其他旧用工业品销售，因此为了研究火柴需求量的变化趋势，需研究影响火柴销售量的变量之间的关系。
) v+ h6 Z' o# @; Y4 [; ^& F    影响火柴销售量的主要指标有：煤气、液化气户数、卷烟销售量、蚁香销售量、打火石钥售量。调查从1963年至1982年共20个年头的数据，即选择的变量数p=4、样本数N＝20，进行因子分析。
6 @9 q; q' e1 M0 ?9 Y原始数据见表10-1。
1 e0 r& r$ j" S   表10-1  影响火柴销售量的原始数据

编号	煤气液化气户数	卷烟销售量	蚊香销售量	火石销售量
1	11.11	14.34	5.10	20.41 [/td]
2	13.63	17.15	4.90	26.94
3	16.99	17.45	5.06	25.75
4	21.12	18.05	5.32	25.98
5	24.86	20.42	7.48	16.18
6	25.68	23.60	10.10	4.18
7	25.77	23.42	13.31	2.43
8	25.88	22.09	9.49	6.50
9	27.43	21.43	11.09	25.78
10	29.95	24.96	14.48	28.16
11	33.53	28.37	16.97	24.26
12	37.31	42.57	20.16	30.18
13	41.16	45.16	26.39	17.08
14	45.73	52.46	27.04	7.39
15	50.59	45.30	23.08	3.88
16	58.82	46.80	24.46	10.53
17	65.28	51.11	33.82	20.09
18	71.25	53.29	33.57	21.22
19	73.37	55.36	39.59	12.63
20	76.68	45.00	48.49	11.17

原始数据经标准化处理后，计算其相关系数矩阵，如表10.2
) t* Y5 ^+ T! p# e1 n. z% B, k  S表10.2 相关系数矩阵

	1	2	3	4
1	1	0.939	0.964	-0.2661
2	0.939	1	0.9407	-0.2754
3	0.964	0.9407	1	-0.2568
4	-0.2661	-0.2754	-0.2568	1

由表10.2,变量x1、x2、x3三个变量彼此间的相关性很强，而变量x4与前三个变量
表现为负相关。
    求解相关矩阵的特征方程，得到特征向量矩阵为：
8 M0 I8 [- g: T8 y5 o( V [img=512,124]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
4 y3 ]1 X0 k: u: G: ~( Y相应的特征值、占总体的百分比、累计百分比如表10.3所示。
1 Y" O* j8 n+ X) L2 H0 p表10.3特征值及百分比表

序号	特征值	占总体的百分比	累计百分比
1	3.0020	75.05	75.05
2	0.8941	22.35	97.40
3	0.0679	1.69	99.10
4	0.0358	0.89	100

根据表10.3选择两个主因子（第一主因子和第二主因子），它们反映的信息量已占
总体信息量的97.4％，即几乎完全反应了四个变量的总信息量。由此选出的两个主因子，得到因子在和矩阵为·
8 W9 s9 s1 |  W[img=224,120]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]       [img=76,56]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
因子模型（对特殊因子忽略不计)为:
; V" m; q$ k5 D: }+ w[img=296,124]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
- H) P2 z& o4 D4 {" D, S从因子载荷矩阵看出，第一主因子主要由前三个变量所决定,而前三个变量彼此间相关性很强，因此在作火柴销售预测时，只需选取三个变量中的任一个变量就可以了。第二个主因子主要由第四个变量所决定，该变量在主因子上的载荷为0.9213。
附SPSS操作说明：
    在Statistics菜单下选择Data Deduction，然后选择Factor过程。Factor过程的选择项共分为五类，分别使用主对话框中的五个子对话框按钮展开相应的子对话框进行选择。
  ☆Descriptives展开相应的子对话框可以选择单变量的描述统计量和初始分析结果。
  ☆Extraction展开相应的子对话框可以选择不同的提取公因子的方法和控制提取结果的判据。
  ☆Rotation展开相应的子对话框可以选择因子旋转方法。
- z6 r+ v8 |5 u! t7 A4 }3 ^  ☆Scores展开相应的子对话框可以要求计算因子得分，选择显示或作为新变量保存。
7 e* p( I1 p; b9 [& i" c8 H  ☆Option展开相应的子对话框可以进一步选择各种输出项。
; y- ^4 K; n) d( k6 v9 c/ ^$ u: H    下面就这几个子对话框分别说明其中的选择项的含义或功能。
8 w5 J6 ?! |7 b7 M8 a$ s1.Descriptives 子对话框中的选择项
! r9 i$ w8 V) S8 L6 e3 B( i; e; M    描述统计量分的选择项为两组:
7 J$ ^5 ?" H( m- p4 V( M3 C6 D. E（1）Statistics 统计量组，共有两项供选择：
. a2 _, S6 l4 v0 [% `) R   ①Univariate descriptives 单变量描述统计量，选择此项可以输出参与分析的各原始变量的均值、标准差等。
8 S* z2 Q) S3 q% D+ H: Q. `9 z   ②Initial solution 初始分析结果，选择此项可以给出原始变量的公因子方差、与变量数目相等的因子、各因子的特征值、各因子特征值占总方差的百分比以及累积百分比。
4 K% l: N, d! B( K$ g5 B. l, x（2） Corelation Matrix 相关矩阵组
   ①Coefficients 相关系数选择此项给出原始变量间的相关系数矩阵。这是分析的基础。
   ②Significance levels 显著性水平选择此项给出每个相关系数相对于相关系数为0的假设检验的概率水平。
   ③Determinant 相关系数矩阵的行列式。
* e) R9 S/ @- L# |/ f   ④Inverse 相关系数矩阵的逆矩阵。
  V3 k8 K5 r3 Y9 }: q: i   ⑤Reproduced 再生相关阵，选择此项给出因子分析后的相关阵，还给出残差，即原始相关与再生相关之间的差值。
   ⑥Anti-image 反映象相关阵。包括偏相关系数的负数；反映象协方差阵，包括偏协方差的负数；在一个好的因子模型中除对角线上的系数较大外，远离对角线的元素应该比较小。
" ~6 B" f) I2 k1 ?; Q, R* ^   ⑦KMO and Bartlett's test of sphericity KMO和球形Bartlett检验。选择此项给出对采样充足度的Kaisex-Meyer-Olkin测度。检验变量间的偏相关是否很小。Barllet球形检验，检验的是相关阵是否是单位阵。它表明因子模型是否是不合适宜的。
    Continue按钮确认选择，返回主对话框；Cancel按钮使选择作废，返回主对话框。
2.Extraction 因子提取子对话框
    有关因子提取的选择项共分四部分。
（1）因子提取方法选择项
" g% a+ @3 j1 ^, A    子对话框第一项Methed：是一组指定提取方法的选择项。单击矩形框右面的箭头可以展开提取方法选择项表，提供七种提取方法可以选择：
   ①Principal components 主成份法。该方法假设变量是因子的纯线性组合。
3 N+ l2 [& V/ D9 F6 E4 S$ u, y   ②Unweighted least square 不加权最小平方法。该方法使观测的和再生的相关阵之差的平方最小。
8 ?9 e' u  G! U% b! x: I3 a   ③Generalized least squares用变量的单值加权，使观测的和再生的相关阵之差的平方最小。
/ t" }  h  Q  q3 q- @  Q0 u   ④Maximum likelihoud 最大似然法。此方法不要求多元正态分布。给出参数估计，如果样本来自多元正态总体它们与原始变量的相关阵极为相似。
! F7 H/ U" ^- m   ⑤Principal Axis factoring 使用多元相关的平方作为对公因子方差的初始估计。
   ⑥Alpha α因子提取法
! }7 |) o) c! i  }! M! ~   ⑦Image 映象因子提取法。是根据变量映象的概念提取公因子的方法。把一个变量看
作其它各变量的多元回归。
（2） Extract 控制提取进程和提取结果的选择项
    理论上因子数目与原始变量数目相等，但因子分析的目的是用少量因子代替多个原始变量，选择提取多少个因子由本组选择项决定。
+ o1 t  `- R5 T   ①Eigenvalues over 该选择项指定提取的因子的特征值。在此项后面的矩形框中给出系统默认值为1即要求提取那些特征值大于1的因子。指定特征值决定提取因子数目的方法是系统默认的方法。
( E. m' j5 z/ j9 }2 o% e   ②Number of factors 该选择项指定提取公因子的数目。用鼠标单击选择此项后，将指定的数目键人到该选择项后面的矩形框中。
! A6 D5 e3 g$ T( e' E（3）Display 指定与因子提取有关的输出项
& W; h1 w0 L# k3 l4 R3 `& e/ c% _: a    ①Unrotated factor solution 要求显示未经旋转的因子提取结果，此项为系统默认的输出项。
' |. s0 Z& {4 _1 z& n+ ]    ②Screen plot 要求显示按特征值大小排列的因子序号与特征值为两个坐标轴的碎石图。
（4）Maximum iterations for Convergence
( x; R8 Q; J" Y6 K    因子分析收敛的最大迭代次数，系统默认的最大迭代次数为25。
    Continue 按钮确认所有选择，返回主对话框；Cancel按钮作废本次所有选择，返回主对话框。
& d6 w  L: v- O- h  V1 H% G3. Rotation 旋转方法选择子对话框
" H+ h% W) q, W1 k1 s6 Q    Rotation旋转方法选择子对话框选择项分三组。
（1） Methed 旋转方法选择项
( N7 {5 ~1 p# u. [! a   ①None不进行旋转。此为系统默认的选择项。
   ②Varimax 方差最大旋转
   ③Equamax 平均正交旋转
$ ^" t  u( N" |0 V2 k9 F   ④Quartmax 四次方最大正交旋转
! j# ?2 W: e7 c1 n8 v  f. ?   ⑤Direct  Oblimin 斜交旋转，指定此项可以在下面的矩形框中键入Delta值，该值应该在0～1之间，是因子映象自相关的范围。0值产生最高相关因子，大负数产生旋转的结果与正交接近。
0 z, }9 L' h# r2 ^（2）Display 有关输出显示的选择项
   ①Rotated solution 旋转结果。指定此项将对正交旋转显示旋转后的因子矩阵模式、因子转换矩阵；对斜交旋转显示旋转后的因子矩阵模式、因子结构矩阵和因子间的相关阵。
* d* i+ s0 j, v: `$ x. d   ②Loding plot(s) 因子载荷散点图。指定此项将给出以两两因子为坐标轴的各变量的载荷散点图。如果有两个因子，给出各原始变量在Factorl-Factor2坐标系中的散点图，如果多于两个因子，则给出三维因子载荷散点图。
  ~# M. D5 V. W) Z（3） Maximum iterations for Convergence：旋转收敛的最大迭代次数。系统默认值为2。可以在此项后面的矩形框中键入指定值。
4.Factor scores 有关因子得分的选择项
    有关因子得分的选择项在Factor scores子对话框中。
    选择项共分三组：
, N: Y+ j/ d3 u7 a* J; o$ Q（1） Save as variables 选择此项，将因子得分作为新变量保存在数据文件中。程序运行结束后，在数据窗中显示出新变量。对每一次分析产生一组新变量，用最后一个数字字符表示分析的顺序号。每次分析中产生多少个因子，就生成多少个新变量，因子序号占倒数第三个字符的位置。倒数第二个字符为“_”。在输出窗中给出对因子得分的命名解释。
（2） Methed 指定计算因子得分的方法。可供选择的方法有：
   ①Regression 回归法。其因子得分的均值为0。方差等于估计因子得分与实际因子得分之间的多元相关的平方。
  H/ v7 L, g- X' R0 c; A: I   ②Bartlett巴特利特法。因子得分均值为0。
3 z+ _9 c3 }3 w: l2 @& D   ③Andeson-Rubin 安德森一鲁宾法。是为了保证因子的正交性而对巴特利特因子得分的调整。其因子得分的均值为0，标准差为1。
" ^6 P. g8 G+ Y+ I4 U3 r（3）Display factor score coefficient matrix 选择此项将在输出窗中显示因子得分系数矩阵。是标准化的得分系数，原始变量值进行标准化后，可以根据该矩阵给出的系数计算各观测量的因子得分。
0 ]& Q  y) N2 I4 _: X$ W5.Options 有关输出的选择项
    有关输出方式以及对带有缺失值的观测量处理方法的选择项均在Options子对话框中。
3 ]0 p( Y( m- g5 J! |" g4 |（1）Missing Values 缺失值处理方法选择项共有三个：
   ①Exclude cases listwise 选择此项，在分析过程中对那些变量中有缺失值的观测量一律剔除。所有带有缺失值的观测量都不参与分析。
   ②Exclude cases pairwise 选择此项，成对剔除带有缺失值的观测量。换句话说，在计算两个变量的相关系数时，只把这两个变量中带肴缺失值的观测量剔除。如果一个观测量在正进行相关系数的变量中没有缺失值，其它变量中带有缺失值，那么不影响当前相关系数的计算。
   ③Replace wlth mean 选择此项，用该变量的均值代替工作变量的所有缺失值。
（2）Coefficient DiDisplay Format这一组选择项决定相关系数的显示格式：
  E7 p2 z1 Q' W& z8 h' J9 S, K" O   ①Sorted by size 选择此项，系数按其数值的大小排列。
/ y: w; n8 H6 U# K' ~   ②Suppress absolute values less than：选择此项，不显示那些绝对值小于指定值的相关系数,选择此项需要在该项右面的矩形框中键人临界值。系统默认的临界值为0.10。
[img=584,5]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
10.3 因子载荷的统计意义
    在因子模型中，假定原始数据已标准化，即原始变量的均值为0，方差为1。并且假定各公共因子和特殊因子都已标准化，即平均值为0，方差为1，在此基础上，我们可进一步讨论与因子载荷有关的一些量的统计意义。
1.因子载荷的统计意义
. b; M; M  }' _- A) H' f因子载荷矩阵：
   [img=272,130]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
因子载荷aij表示第i个变量和第j个公共因子的相关系数。
即 rxiFj=aij
aij越大表示公因子Fj与变量的关系越密切，如(10.2.7)式中所表示的载荷矩阵，第一主因子F1与变量x1、x2、x3都较留切，而F1与x4的关系不够密切。
- n$ s: `5 Q; }' i2．变量共同度的统计定义
- {$ b8 Y" y& ^9 j5 [, M/ j    因子载荷矩阵中各行元素的平方和。
[img=256,126]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
    称为变量x1、x2、...xp的共同度。
    计算每一个变量xi的方差：
  Dxi=ai12DF1+ai22DF2+...+aim2DFm+ai2Dεi
9 N5 R. y& L  {  N& o" @由于假设原始变量xi和主因子、特殊因子也都进行标准化处理，所以有：
( I3 c/ ^0 J) Z" e- B# v    1=hi2+ai2
即变量xi的方差由两部分组成，第一部分为共同度hi2，它是全部公因子对变量xi的总方差的贡献，如公因子的方差接近于1，则说明该变量的几乎全部原始信息都被所选取的因子说明了。如载荷矩阵中h1=0.97，即变量x1的97％的信息量被第一主因子和第二主因子说明。
# W4 p+ h3 O% r$ T6 \. r    第二部分是特殊因子的方差，仅与变量xi本身的变化有关，它是变量xi的方差为1的补充值，当不考虑特殊因子时，公因子方差hi＜1(i=1,2,...p)。
  `1 C8 R& l) N6 ]6 ]    公因子方差的意义在于说明当由原始变量空间转为因子空间后，保留原来各变量的信息有多少，当hi越接近于1时，说明空间转化性质越好。
1 K7 G: Q& u' H6 b  I0 X3．公因子Fj的方差贡献统计意义
2 m# B! O1 @9 ~3 N; d/ o0 C    各列元素的平方和Sj为公因子Fj的方差贡献。
! U0 E3 u% `. o4 ~& Q& ?   [img=276,122]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
$ q: P4 @' c: ]3 f2 C, n# s/ @Sj是公因子Fj对诸原始变量所提供方差贡献的总和。它是衡量公因子相对重要性的指标，它等于公因子Fj所对应的特征值，即
   [img=108,52]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
6 ]* L4 k+ Y( j[img=584,5]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
8 E. `% U! \( i10.4 方差最大正交旋转
     主因子分析的目的不仅是找出主因子，更重要的是知道每个主因子的意义。但是用上述方法所求出的主因子解，初始因子载荷矩阵并不满足“简单结构准则”，各因子的典型代表变量不很突出，因而容易使因子的意义含糊不清，不便于对因子进行解释。为此须对因子载荷矩阵施行旋转，使得因子载荷的平方按列向0和1两极转化，达到其结构简化的目的。
* l9 E! l4 f$ S6 b, J    采用方差最大旋转（Varimax rotation)，使因子载荷矩阵中，各因子载荷值的总方差达到最大作为因子载荷矩阵简化的准则。这里总方差最大，而不是某一因子方差极大。即是说，如果第i个变量在第j个公因子上的载荷经过“方差极大”旋转后其值增大或减小，意味着这变量在另一些因子上的载荷要缩小或增大。所以，“方差极大”旋转是使载荷按列向0、 1两极分化，同时也包含着按行向两极分化。
设因子载荷矩阵为：
& J; a% h3 G6 X/ g; S6 I) r/ m  [img=276,123]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
bij为经过方差极大正交旋转后所得到的因子载荷矩阵B的元素。
某一因子Fj的简化可由因子载荷平方的方差来表示。
! l! X4 I* b3 H( i8 d. e2 G+ S; F  [img=364,130]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
' \: M9 [- }! S% t. ~9 C- D如果Vj为极大值，则此Fj因于具有简化性，它的载荷值或是趋于1或是趋于0。
    例10-1中的载荷矩阵进行方差极大旋转，经过二次正交旋转即收敛。
    第一次旋转的总方差值:V1=0.260395
    第二次旋转的总方差值:V2=0.347423
! c) K. E- w/ \    旋转后的正交载荷矩阵为：
    [img=248,129]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
    旋转结果明显将因子载荷向1和0两极方向分化，这样便于给予明确的经济意义上的解释。
) K! n/ R% F& {5 o: M3 U: x% i[img=584,5]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
' e5 L; B$ C: {# d+ j% i+ G10.5 Promax斜旋转
      在方差极大旋转过程中，因子轴互相正交，始终保持初始解中因子间不相关的特点，然而在社会学、经济学、心理学等科学领域内，斜交因子是普遍规律，即相互影响的各种因素不大可能是彼此无关的，各种事物变化的各种内在因素之间始终存在着错综复杂的联系，因此需要介绍斜交因子解，即将变量用相关因子进行线性描述，使得到的新因子模型最大程度地模拟自然模型。
1．斜交因子模型和斜交因子解
    如果p个相关变量x1、x2、...xp，用M个相关斜交公因子T1、T2、...Tm表示。
    斜交因子模型为：
0 y9 ^$ J; Z; Y* _0 x   [img=324,118]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
    T1为斜交公因子，bij为斜交因子载荷，它表示变量xi的向量pi在斜因子轴Tj上的坐标，在斜交因子轴情况下，坐标和投影，即模型和结构之间有差别．如图10-1所示，T2是斜因子轴，pi为因子空间中表示变量xi的向量(pi的长度|Opi|=1)，pi在OT1T2上的坐标为OQ、OR，坐标值可正、可负，且绝对值可以大于1。而pi在OT1T2上的投影为OM、ON，坐标值也可正、可负，但绝对值不能大于1。而在正交因子空间中，坐标值和投影值一致，即因子模型和因子结构是一回事，无需区分这两者。所以一个完全斜交因子解必须包括因子模型和因子结构两大部分。
! r6 y' l8 [8 ^4 G    [img=264,226]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
  图10-1
例1中两个正交因子（M＝2）经过Promax斜旋转，得到斜交因子相关阵为：
0 ?1 J8 ]( O3 v  X& s9 B. H[img=240,55]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
斜交因子结构如表10.4所示：
. v" _5 l. r* z: {' s, m# S6 |   表10.4 斜交因子结构
2 ^- b0 m; a& m# O* E. Y

变量	第一斜交因子	第二斜交因子
煤气液化气户数	0.9999	-0.2610
卷烟销售量	0.9999	-0.2740
蚊香销售量	0.9999	-0.2511
火石销售量	-0.2709	0.9999

    由表10.4知，斜旋转的结果与正文旋转结果基本一致，第一主因子仍主要由变量x1、x2、x3表征，第二主因子仍主要由变量x4所表征。
- i* g* J, O: y* @5 W. f5 {. p, |! A斜交因子模型如表10.5所示：
/ V& x# x( G% }0 k) C# a& Y7 L    表10.5 斜交因子模型
4 `, I. r$ s3 n" g* i$ u

变量	第一斜交因子	第二斜交因子
煤气液化气户数	1.00002	0.0011
卷烟销售量	0.9965	-0.0128
蚊香销售量	1.0030	0.0117
火石销售量	-0.0094	0.9974

: b: V. Q) }/ P5 x    因子模型矩阵是把变量表为斜因子的线性组合，所以例1的斜交因子模型为。
. e+ s/ b. V5 L6 Q9 V  ^+ f* ]5 qx1＝1.0002T1＋0.0011T2
x2＝0.0095T1-0.0128T2
x3＝1.0030T1＋0.0117T2
" r' s8 F  K+ ?7 C3 Wx4＝-0.0094T1＋0.9974T2
2 K6 w" f1 W% y0 @    由此说明：煤气、液化气户数主要由斜因子1提供，而打火石销售量主要由斜因子2提供。
2 Z  X2 Y' r* ~- o3 ]. J/ ?10.6 因子得分
( @' i& Q% D! T3 }      因子分析是将变量表为公共因子的线性组合。由于公因子能充分反映原始变量的相关关系，用公因子代表原始变量时，更有利于描述研究对象的特征，因而、反过来将公共因子表示为变量的线性组合（如同主成分分析一样），即用：
8 d% S, k+ }. d% r6 R+ l/ H0 xFj=βj1x1+βj2x2+...βjpxp
2 z  F1 Z% D8 Q1 }! E, s) z* B    (j=1,2,...,m)
来计量各个样品的公因子得分。上式为因子得分函数。
    由于上式中方程的个数少于变量个数，因此，只能在最小二乘意义下对因子得分进行估计。Thomson（1939年）假设因子可以对p个变量作回归，即建立回归方程：
由于变量和主因子均已标准化，所以有：
) A8 ?5 c6 l: E  Z4 B$ Z: |- S# x   [img=332,68]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
1 V7 W1 O! U" l: u9 B. |由最小二乘估计有：
4 Y7 |- S. _' b1 R) m; z1 [   [img=148,34]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
这里
   [img=280,117]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
为原变量的相关系数矩阵。
   A'=(aij)mxn
& [6 |4 M6 s) {4 v# x% h! L) s/ {当因子正交时，A’为旋转后的因子载荷矩阵A的转置，当因子斜交时,A'为因子结构矩阵的转置。
: X6 s; J4 E3 Z( F; ]   [img=212,63]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
例1中数据的正交因子得分表如表10.6所示。
, _( t; K) q+ M3 {! g% R# ]( |5 H表10.6 正交因子得分表

样本编号	[img=24,33]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]	[img=24,33]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
1	-1.3101	0.2160
2	-1.0817	0.9823
3	-1.0342	0.8543
4	-0.9392	0.8939
5	-0.9200	-0.2179
6	-0.9501	-1.5888
7	-0.9044	-1.7735
8	-0.9673	-1.3245
9	-0.5987	0.9241
10	-0.3305	1.2305
11	-0.1762	0.8058
12	0.4812	1.5381
13	0.5532	0.0672
14	0.6959	-1.0321
15	0.4078	-1.4581
16	0.7313	-0.6557
17	1.3537	0.5286
18	1.5278	0.6765
19	1.6269	-0.2774
20	1.8347	-0.3891

    根据因子得分表，以F1和F2为因子轴做因子得分图、如图10-2所示。它将20个抽样样本分为三类、一类以F1得分是负值、F2得分是正值为特征；另一类是以F1得分是正值为主要特征；再一类是以F1和F2得分都是负值为特征。由因子得分图可提供一个样本分类的参考图形。
[img=536,520]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]

10.7 应用实例
! b$ O) E$ w5 N0 L 例10-2 影响小春粮食产量的指标有九个：小春粮食播种面积、小麦播种面积、小麦良种推广比例、化肥用量、肥猪出栏、农业人口、耕牛数、小麦抽穗扬花期气温、小麦抽穗扬花期降水量。抽取从1952年～1957年和从1960年～1981年共28个年头的数据，采用因子分析研究变量之间的相关关系。原始数据如表10.7所示。
- u/ T! p$ e. B) H$ [    表10.7 原始数据表
! K( X" T1 m" f& u1 o% N( \

编号	小春播种面积	小麦播种面积	小麦良种比例	化肥用量	肥猪出栏	农业人口	耕牛数	抽扬期气温	抽扬期水量
1	4088	1569	26.9	.37	648	5694.0	459.0	17.0	46.9 [/td]
2	4216	1685	28.0	1.03	1146	5839.0	465.1	15.1	69.7
3	4233	1702	28.1	2.30	1308	5894.3	490.5	15.2	73.1
4	4284	1863	30.4	2.47	1219	6008.9	492.7	14.7	58.3
5	4639	2034	32.9	5.46	1053	6096.1	509.7	16.4	63.6
6	4786	2073	32.5	6.86	1080	6188.5	515.6	16.1	52.0
7	4995	2708	43.4	17.16	324	5552.3	503.5	15.8	40.5
8	4919	2608	42.4	18.28	214	5487.4	473.2	15.5	64.1
9	4672	2351	40.3	24.55	383	5672.4	473.8	15.2	59.7
10	4328	2097	38.8	35.56	754	5849.7	502.0	16.1	69.4
11	4295	2067	38.8	44.65	1256	6042.9	523.9	16.8	61.9
12	4170	1920	36.8	67.97	1524	6227.7	552.4	15.8	55.7
13	3909	1806	37.0	119.54	1698	6413.0	574.8	16.3	28.3
14	4063	2003	39.4	105.25	1728	6627.2	587.4	15.0	62.0
15	3887	1900	39.1	58.72	1584	6875.8	600.1	14.2	85.6
16	3872	1890	39.1	98.53	1629	7115.7	605.2	16.4	74.6
17	3937	1964	39.9	123.30	1610	7356.9	617.5	13.6	49.4
18	4039	2094	41.5	168.40	1539	7567.7	641.7	16.3	33.6
19	4354	2496	45.9	168.15	1898	7770.4	641.0	15.4	71.6
20	4465	2604	46.7	208.38	1974	8002.2	633.0	17.1	83.0
21	4436	2605	49.9	205.81	1895	8205.5	625.1	15.8	62.7
22	4367	2580	50.2	231.38	1807	8402.5	615.2	16.1	38.5
23	4590	2784	51.6	243.96	1754	8504.4	601.9	13.8	42.1
24	4663	2808	54.2	365.28	1831	8566.0	582.3	15.4	81.0
25	5041	3060	54.6	479.23	2155	8575.2	594.8	15.9	38.4
26	5257	3428	58.7	542.86	2736	8594.9	626.0	15.4	51.5
27	4819	3261	60.9	570.00	3100	8615.0	627.0	15.2	39.2
28	4486	3055	80.0	580.00	3200	8670.0	627.0	13.5	53.9

原始数据经标准化处理后，得到变量的相关系数矩阵如表10.8所示。
- k" [; ?$ N; t/ j4 W$ l    表10.8 相关系数矩阵
9 M' L/ k# F0 p

	小春播种面积	小麦播种面积	小麦良种比例	化肥用量	肥猪出栏	农业人口	耕牛数	抽扬期气温	抽扬期水量
小春播种面积	1.000	.777	.450	.441	.081	.228	-.069	.021	-.178
小麦播种面积	.777	1.000	.877	.838	.561	.716	.496	-.193	-.212
小麦良种比例	.450	.877	1.000	.905	.741	.804	.654	-.352	-.186
化肥用量	.441	.838	.905	1.000	.873	.876	.680	-.263	-.261
肥猪出栏	.081	.561	.741	.873	1.000	.857	.802	-.299	-.129
农业人口	.228	.716	.804	.876	.857	1.000	.869	-.249	-.148
耕牛数	-.069	.496	.654	.680	.802	.869	1.000	-.189	-.119
抽扬期气温	.021	-.193	-.352	-.263	-.299	-.249	-.189	1.000	-.022
抽扬期水量	-.178	-.212	-.186	-.261	-.129	-.148	-.119	-.022	1.000

    计算相关系数矩阵的特征值和特征向量。由表10.9，按累计比例数0.94选取主因子，可选取4个主因子，即4个主因子所包含的信息量占总体信息量的94.26％。
  表10.9 特征值及占总体的比例数表
[img=608,234]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
9 {5 w# ?3 t' t$ H% j% M相应地得到因子载荷矩阵如表5.10所示。
  a/ M6 [: [) s$ E+ K" U* q7 k    表10.10 因子载荷矩阵
7 v( S4 N& m; o6 M- ][img=608,234]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
+ ~- P  Y+ D- |7 U. G8 p0 U$ ^' J2 E  因子模型为：
[img=464,88]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
/ Q4 z5 }1 B; K0 g0 ^) r* |) s[img=492,182]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
/ R) b! M- d! o8 ?1 v' t5 I# p& Z    为了进一步简化结构，进行方差极大旋转。经过4次旋转，得到因子载荷矩阵的方差最大值为：
2 ^9 C# p& M$ P6 I+ G    V=0.457350387
% j) H. G  ^$ V  S" B  ~    得到正交旋转矩阵如表10.11 所示。
% S  b6 Y: c5 L) Z    表10.11 方差极大正交旋转结果表
[img=776,296]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
1 [' X* @4 J' R
/ O- p6 E) ?1 ^  从表10.11看出，第一主因子主要由变量4、变量5、变量6、变量7所决定，这四
个变量所反映的信息量几乎占总体信息量的57％；第二主因子主要由变量1所决定；第三主因子主要由变量9所决定、第四主因子主要由变量8所决定。这样，将9个变量简化为四个彼此独立的因子，即将9个变量分为四大类，第一类的代表变量是耕牛数，该变量有最大的载荷值（0.9955），第二类的代表变量是小春播种面积，该变量的载荷值为0.9957；第三类的代表变量是小麦抽穗扬花期降水量，其载荷值为0.9904；第四类的代表变量是小麦抽穗扬花期气温，其载荷值为0.9860。
    为了对样本分类和进一步解释，对前四个主因子计算其因子得分，得到因子得分表，如表10.12所示。
! B* ~$ T; l) J. b- Y    根据因子得分表，取F1和F2两个因子轴作图，得到因子得分图(图10-3)。
    由因子得分图，将样本分为三个类，第一类主要是以F2为正值，它包含了13个样本；第二类以F1为正值，它包含了9个样本；第三类是以F1和F2都为负值，它包含了5个样本。
* [# R  G  W! _$ e' S2 _   表10.12 因子得分表
[img=808,760]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
[img=476,426]mhtml:file://\\192.168.41.240\2010数学建模资料汇总\光盘4\729－补充材料多元\多元讲义\5-因子分析\第十章 因子分析.mht![/img]
    图10-3 因子载荷平面图
! X" }6 V8 p' K7 @$ e* [" Z1 c
这是在网上找的。呵呵，和大家分享啦！呵呵
& X7 l" y: U* T

web_star

留个脚印```````

9112001go

我回不抢呢 考虑再三 还是不抢了吧 ^_^

qysjg

顶顶更健康，越顶吃的越香。

lihomlh

很牛。。。。。。。。。。。。。

mengxue11180611

顶顶更健康，越顶吃的越香。

xy178836

楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  

hexingfeng

不错。。。。。

networm

楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   

andrew01

强烈支持。楼主万岁