1、1因子分析因子分析2 1 1 引言引言 因子分析(factor analysis)是一种数据简化的技术。它通过研究众多变量之间的内部依赖关系，探求观测数据中的基本结构，并用少数几个假想变量来表示其基本的数据结构。这几个假想变量能够反映原来众多变量的主要信息。原始的变量是可观测的显在变量，而假想变量是不可观测的潜在变量，称为因子。 例如，在企业形象或品牌形象的研究中，消费者可以通过一个有24个指标构成的评价体系，评价百货商场的24个方面的优劣。3 但消费者主要关心的是三个方面，即商店的环境、商店的服务和商品的价格。因子分析方法可以通过24个变量，找出反映商店环境、商店服务水平和商品价格的三个潜在
2、的因子，对商店进行综合评价。而这三个公共因子可以表示为： iiiiiiFFFx33221124, 1i 称 是不可观测的潜在因子。24个变量共享这三个因子，但是每个变量又有自己的个性，不被包含的部分 ，称为特殊因子。321FFF、i4注：注： 因子分析与回归分析不同，因子分析中的因因子分析与回归分析不同，因子分析中的因子是一个比较抽象的概念，而回归因子有非常明子是一个比较抽象的概念，而回归因子有非常明确的实际意义；确的实际意义； 主成分分析分析与因子分析也有不同，主成主成分分析分析与因子分析也有不同，主成分分析仅仅是变量变换，而因子分析需要构造因分分析仅仅是变量变换，而因子分析需要构造因子模型
3、。子模型。 主成分分析主成分分析: :原始变量的线性组合表示新的原始变量的线性组合表示新的综合变量，即主成分；综合变量，即主成分； 因子分析：潜在的假想变量和随机影响变因子分析：潜在的假想变量和随机影响变量的线性组合表示原始变量。量的线性组合表示原始变量。5 2 因子分析模型因子分析模型 一、数学模型一、数学模型 设 个变量，如果表示为iX), 2 , 1(pip11iiiimmiXa Fa F)(pm 11111211122212222212mmpppppmpmXFXFXF或XAF或6 称为 公共因子，是不可观测的变量，他们的系数称为因子载荷。 是特殊因子，是不能被前m个公共因子包含的部分。
4、并且满足：mFFF,21iIFD111)(cov( , )0,F,F即不相关；mFFF,21即 互不相关，方差为1。722221)(pD即互不相关，方差不一定相等， 。), 0(2iiN8用矩阵的表达方式X- = AF+( )EF0( )E0( )VarFI22212( )(,)pVardiag1 11212 122212()()()()()()cov()()()()()ppppppE FE FE FE FE FE FEE FE FE FF,F09二、因子分析模型的性质 1、原始变量X的协方差矩阵的分解X- = AF+()( )( )VarVarVarX- = AF A+x = AA +DA是
5、因子模型的系数22212( )(,)pVardiagD D的主对角线上的元素值越小，则公共因子共享的成分越多。10 2、模型不受计量单位的影响 将原始变量X做变换X*=CX,这里Cdiag(c1,c2,cn),ci0。)C(X-) = C(AF+CXC+CAF+C*XC+CAF+C*X + A F +*FF11*()EF0*( )E0*()VarFI*2222221122( )(,)ppVardiag ccc* *cov()()EF ,F 012 3、因子载荷不是惟一的 设T为一个pp的正交矩阵，令A*=AT，F*=TF，则模型可以表示为*X+ A F +()ET F0( )E0*()()(
6、)VarVarVarFT FTF TI22212( )(,)pVardiag*cov()()EF ,F 0且满足条件因子模型的条件13 三、三、 因子载荷矩阵中的几个统计特征因子载荷矩阵中的几个统计特征 1 1、因子载荷、因子载荷a aijij的统计意义的统计意义 因子载荷 是第i个变量与第j个公共因子的相关系数 ija模型为 imimiiFaFaX11 在上式的左右两边乘以 jF,再求数学期望 )()()()()(11jijmimjjijjijiFEFFEaFFEFFEaFXE 根据公共因子的模型性质，有ijFxji （载荷矩阵中第i行，第j列的元素）反映了第i个变量与第j个公共因子的相关重
7、要性。绝对值越大，相关的密切程度越高。14 2 2、变量共同度的统计意义、变量共同度的统计意义定义：定义：变量 的共同度是因子载荷矩阵的第i行的元素的平方和。记为iX统计意义统计意义：imimiiFaFaX11两边求方差 )()()()(2112imimiiVarFVaraFVaraXVarmjiija1221 所有的公共因子和特殊因子对变量 的贡献为1。如果 非常靠近1， 非常小，则因子分析的效果好，从原变量空间到公共因子空间的转化性质好。iXmjija122imjija12。mjijiah12215 3 3、公共因子、公共因子 方差贡献的统计意义方差贡献的统计意义jF因子载荷矩阵中各列元素
8、的平方和 称为某一公共因子 对诸变量所提供的方差贡献和。衡量的相对重要性。piijjaS12jFjF16 3 3 因子载荷矩阵的估计方法因子载荷矩阵的估计方法 设随机向量 的均值为 ，协方差为 , 为的特征根， 为对应的标准化特征向量，则pxxx,21x021pp21u,u,u12p = UUAA +D主成分分析法主成分分析法17 上式给出的 表达式是精确的，然而，它实际上是毫无价值的，因为我们的目的是寻求用少数几个公共因子解释，故略去后面的p-m项的贡献，有21111mmmmmmp1122ppu uu uu uuuu up2uuuuuuppp21122111100p212ppuuuuuu18
9、 12 mmm1122AA +Du uu uu uD1121122 mmp mpmm p2uuuuuDAADu 上式有一个假定，模型中的特殊因子是不重要的，因而从 的分解中忽略了特殊因子的方差。22212(,)pdiagD其中221miiiijjsa19 例例 假定某地固定资产投资率 ，通货膨胀率 ，失业率 ，相关系数矩阵为试用主成分分析法求因子分析模型。1x2x3x15/25/15/215/15/15/1120 特征根为： 55. 11 85. 02 6 . 03 6 . 0707. 085. 0331. 055. 1629. 06 . 0707. 085. 0331. 055. 1629.
10、 0085. 0883. 055. 1475. 0A707. 0331. 0629. 0707. 0331. 0629. 00883. 0475. 0U548. 0305. 0783. 0548. 0305. 0783. 00814. 0569. 021 可取前两个因子F1和F2为公共因子，第一公因子F1物价就业因子，对X的贡献为1.55。第一公因子F2为投资因子，对X的贡献为0.85。共同度分别为1，0.706，0.706。211814. 0569. 0FFx3212548. 0305. 0783. 0FFFx3213548. 0305. 0783. 0FFFx22 4 因子旋转（正交变换）
11、 建立了因子分析数学目的不仅仅要找出公共因子以及对变量进行分组，更重要的要知道每个公共因子的意义，以便进行进一步的分析，如果每个公共因子的含义不清，则不便于进行实际背景的解释。由于因子载荷阵是不惟一的，所以应该对因子载荷阵进行旋转。目的是使因子载荷阵的结构简化，使载荷矩阵每列或行的元素平方值向0和1两极分化。有三种主要的正交旋转法。四次方最大法、方差最大法方差最大法和等量最大法。23 百米跑成绩 跳远成绩 铅球成绩 跳高成绩 400米跑成绩 百米跨栏 铁饼成绩 撑杆跳远成绩 标枪成绩 1500米跑成绩 1X2X3X4X5X6X7X8X9X10X奥运会十项全能运动项目奥运会十项全能运动项目得分数
12、据的因子分析得分数据的因子分析 24102. 017. 002. 001. 039. 018. 008. 009. 007. 0124. 034. 018. 013. 017. 044. 021. 011. 0124. 033. 023. 039. 024. 036. 020. 0132. 017. 027. 073. 031. 028. 0134. 046. 036. 052. 040. 0129. 019. 049. 063. 0138. 051. 034. 0142. 035. 0159. 0125变量共同度0.6910.217-0.58-0.2060.840.7890.184-0.19
13、30.0920.70.7020.5350.047-0.1750.80.6740.1340.1390.3960.650.620.551-0.084-0.4190.870.6870.042-0.1610.3450.620.621-0.5210.109-0.2340.720.5380.0870.4110.440.660.434-0.4390.372-0.2350.570.1470.5960.658-0.2790.891F2F3F4F1X2X3X4X5X6X7X8X9X10X 因子载荷矩阵可以看出，除第一因子在所有的变量在公共因子上有较大的正载荷，可以称为一般运动因子。其他的3个因子不太容易解释。似乎
14、是跑和投掷的能力对比，似乎是长跑耐力和短跑速度的对比。于是考虑旋转因子，得下表 2627 通过旋转，因子有了较为明确的含义。 百米跑， 400米跑，需要爆发力的项目在 有较大的载荷， 可以称 为短跑速度因子； 铅球， 铁饼和 标枪在 上有较大的载荷，可以称为爆发性臂力因子； 百米跨栏， 撑杆跳远， 跳远和为 跳高在 上有较大的载荷， 爆发腿力因子； 长跑耐力因子。5X1F1F3X7X9X2F6X8X2X4X3F3F4F1X28变换后因子的共同度变换后因子的共同度设 正交矩阵，做正交变换正交矩阵，做正交变换 AB )()(1mlljilppijabBmjmjmlljilijiabh111222)
15、()(B mjmlmjmlmljttjljitilljilaaa1 11 1 122)(2111222Aimlmjmlilljilhaa变换后因子的共同度没有发生变化！变换后因子的共同度没有发生变化！（二）旋转方法（二）旋转方法29变换后因子贡献变换后因子贡献设 正交矩阵，做正交变换正交矩阵，做正交变换 AB )()(1qlljilppijabBpipiqlljilijjabS111222)()(B piqlpiqlqltttjljitilljilaaa1111 122piqlqlljjljilSa1112222)(A变换后因子的贡献发生了变化变换后因子的贡献发生了变化！30 1、方差最大法
16、方差最大法从简化因子载荷矩阵的每一列出发，使和每个因子方差最大法从简化因子载荷矩阵的每一列出发，使和每个因子有关的载荷的平方的方差最大。当只有少数几个变量在某个因子有关的载荷的平方的方差最大。当只有少数几个变量在某个因子上又较高的载荷时，对因子的解释最简单。上又较高的载荷时，对因子的解释最简单。方差最大的直观意义是希望通过因子旋转后，使每个因子上的载荷尽量拉开距离，一部分的载荷趋于1，另一部分趋于0。2122211211ppaaaaaaA221122212122121111FaFaXFaFaXFaFaXppp31cossinsincosT设旋转矩阵为：cossinsincosAATB则coss
17、insincoscossinsincos112112111211ppppaaaaaaaa*2*1*12*11ppaaaa321,2, ;1,2ijijiadip jh令211(pjijiddp这是列和）max)()(1212 mjpijijddV简化准则为：00V令，则可以解出0000cossinsincosT旋转矩阵为：max(8.4.2)123m即：V +V +V+V33 5 因子得分因子得分 （一）因子得分的概念（一）因子得分的概念 前面我们主要解决了用公共因子的线性组合来表示一组观测变量的有关问题。如果我们要使用这些因子做其他的研究，比如把得到的因子作为自变量来做回归分析，对样本进行分
18、类或评价，这就需要我们对公共因子进行测度，即给出公共因子的值。34 人均要素变量因子分析人均要素变量因子分析。对我国32个省市自治区的要素状况作因子分析。指标体系中有如下指标：X1 ：人口（万人） X2 ：面积（万平方公里）X3 ：GDP（亿元） X4 ：人均水资源（立方米/人）X5：人均生物量（吨/人） X6：万人拥有的大学生数（人）X7：万人拥有科学家、工程师数（人） Rotated Factor Pattern FACTOR1 FACTOR2 FACTOR3 X1 -0.21522 -0.27397 0.89092 X2 0.63973 -0.28739 -0.28755 X3 -0.1
19、5791 0.06334 0.94855 X4 0.95898 -0.01501 -0.07556 X5 0.97224 -0.06778 -0.17535 X6 -0.11416 0.98328 -0.08300 X7 -0.11041 0.97851 -0.0724635 高载荷指标 因子命名 因子1X2；面积（万平方公里）X4:人均水资源（立方米/人）X5:人均生物量（吨/人）自然资源因子 因子2X6：万人拥有的大学生数（人）X7：万人拥有的科学家、工程师数（人） 人力资源因子 因子3 X1;人口（万人）X3:GDP(亿元)经济发展总量因子 X1=-0.21522F1-0.27397F2
20、+0.89092F3 X2=0.63973F1-0.28739F2-0.28755F3 X3=-0.15791F1+0.06334F2+0.94855F3 X4=0.95898F1-0.01501F2-0.07556F3 X5=0.97224F1-0.06778F2-0.17535F3 X6=-0.11416F1+0.98328F2-0.08300F3 X7=-0.11041F1+0.97851F2-0.07246F336 Standardized Scoring Coefficients FACTOR1 FACTOR2 FACTOR3 X1 0.05764 -0.06098 0.50391
21、X2 0.22724 -0.09901 -0.07713 X3 0.14635 0.12957 0.59715 X4 0.47920 0.11228 0.17062 X5 0.45583 0.07419 0.10129 X6 0.05416 0.48629 0.04099 X7 0.05790 0.48562 0.04822F1=0.05764X1+0.22724X2+0.14635X3+0.47920X4+0.45583X5+0.05416X6+0.05790X7F1=0.05764X1+0.22724X2+0.14635X3+0.47920X4+0.45583X5+0.05416X6+0.
22、05790X7F2=-0.06098X1-0.09901X2+0.12957X3+0.11228X4+0.07419X5+0.48629X6+0.48562X7F2=-0.06098X1-0.09901X2+0.12957X3+0.11228X4+0.07419X5+0.48629X6+0.48562X7F3=0.50391X1-0.07713X2+0.59715X3+0.17062X4+0.10129X5+0.04099X6+0.04822X7F3=0.50391X1-0.07713X2+0.59715X3+0.17062X4+0.10129X5+0.04099X6+0.04822X737R
23、EGION FACTOR1FACTOR2FACTOR3beijing-0.081694.23473-0.37983tianjin-0.474221.31789-0.87891hebei-0.22192-0.358020.86263shanxi1-0.48214-0.32643-0.54219neimeng0.54446-0.66668-0.92621liaoning-0.205110.463770.34087jilin-0.214990.10608-0.57431heilongj 0.10839-0.11717-0.02219shanghai-0.200692.38962-0.04259前三个
24、因子得分前三个因子得分38 因子分析的数学模型为： mpmppmmnFFFXXX2121222211121121 原变量被表示为公共因子的线性组合，当载荷矩阵旋转之后，公共因子可以做出解释，通常的情况下，我们还想反过来把公共因子表示为原变量的线性组合。 因子得分函数： pjpjjXXF11mj, 1可见，要求得每个因子的得分，必须求得分函数的系数，而由于pm，所以不能得到精确的得分，只能通过估计。39 回归方法 nmnmnnmmnFFFXXX212121222211121121pjpjjXbXbF11mj, 1mmpmmppbbbbbbbbbbbb212122221112111) 思想40)(
25、jiFxijFXEji)(11pjpjiXbXbXEipjpijbb11jpjjipiibbbrrr2121 则，我们有如下的方程组：41pjjjjpjjppppppaaabbb2121212222111211j=1,2,m矩阵为原始变量的相关系数pppppp21222211121142个因子得分函数的系数为第 jbbbjpjj21列为载荷矩阵的第 jaaapjjj21 注：共需要解注：共需要解m次才能解次才能解出出 所有的得分函数的系数。所有的得分函数的系数。43国民生活质量的因素分析 国家发展的最终目标，是为了全面提高全体国民的生活国家发展的最终目标，是为了全面提高全体国民的生活质量，满足
26、广大国民日益增长的物质和文化的合理需求。在质量，满足广大国民日益增长的物质和文化的合理需求。在可持续发展消费的统一理念下，增加社会财富，创自更多的可持续发展消费的统一理念下，增加社会财富，创自更多的物质文明和精神文明，保持人类的健康延续和生生不息，在物质文明和精神文明，保持人类的健康延续和生生不息，在人类与自然协同进化的基础上，维系人类与自然的平衡，达人类与自然协同进化的基础上，维系人类与自然的平衡，达到完整的代际公平和区际公平到完整的代际公平和区际公平( (即时间过程的最大合理性与空即时间过程的最大合理性与空间分布的最大合理化间分布的最大合理化) )。 从从19901990年开始，联合国开发
27、计划署年开始，联合国开发计划署(UYNP)(UYNP)首次采用首次采用“人文人文发展系数发展系数”指标对于国民生活质量进行测度。人文发展系数指标对于国民生活质量进行测度。人文发展系数利用三类内涵丰富的指标组合，即人的健康状况利用三类内涵丰富的指标组合，即人的健康状况( (使用出生时使用出生时的人均预期寿命表达的人均预期寿命表达) )、人的智力程度、人的智力程度( (使用组合的教育成就使用组合的教育成就表达表达) )、人的福利水平、人的福利水平( (使用人均国民收入或人均使用人均国民收入或人均GDPGDP表达表达) )，并且特别强调三类指标组合的整体表达内涵，去衡量一个国并且特别强调三类指标组合
28、的整体表达内涵，去衡量一个国家或地区的社会发展总体状况以及国民生活质量的总水平。家或地区的社会发展总体状况以及国民生活质量的总水平。44在这个指标体系中有如下的指标：X1预期寿命X2成人识字率X3综合入学率X4人均GDP（美圆）X5预期寿命指数X6教育成就指数X7人均GDP指数45 旋转后的因子结构 Rotated Factor Pattern FACTOR1 FACTOR2 FACTOR3 X1 0.38129 0.41765 0.81714 X2 0.12166 0.84828 0.45981 X3 0.64803 0.61822 0.22398 X4 0.90410 0.20531 0.
29、34100 X5 0.38854 0.43295 0.80848 X6 0.28207 0.85325 0.43289 X7 0.90091 0.20612 0.35052 FACTOR1为经济发展因子 FACTOR2为教育成就因子 FACTOR3为健康水平因子46 被每个因子解释的方差和共同度 Variance explained by each factor FACTOR1 FACTOR2 FACTOR3 2.439700 2.276317 2.009490 Final Communality Estimates: Total = 6.725507 X1 X2 X3 X4 X5 0.987
30、530 0.945796 0.852306 0.975830 0.992050 X6 X7 0.994995 0.976999 47生育率的影响因素分析 生育率受社会、经济、文化、计划生育政策等很多因素影响，但这些因素对生育率的影响并不是完全独立的，而是交织在一起，如果直接用选定的变量对生育率进行多元回归分析，最终结果往往只能保留两三个变量，其他变量的信息就损失了。因此，考虑用因子分析的方法，找出变量间的数据结构，在信息损失最少的情况下用新生成的因子对生育率进行分析。 选择的变量有：多子率、综合节育率、初中以上文化程度比例、城镇人口比例、人均国民收入。下表是1990年中国30个省、自治区、直辖
31、市的数据。48多子率（%）综合节育率（%） 初中以上文化程度比例（%）人均国民收入（元） 城镇人口比例（%）0.9489.8964.51357773.082.5892.3255.41298168.6513.4690.7138.2114819.0812.4690.0445.12112427.688.9490.4641.83108036.122.890.1750.64201150.868.9191.4346.32138342.658.8290.7847.33162847.170.891.4762.36482266.235.9490.3140.85169621.242.692.4235.141717
32、32.817.0787.9729.5193317.914.4488.7129.04131321.3615.2489.4331.0594320.43.1690.2137.85137227.349.0488.7639.7188015.5212.0287.2838.76124828.9111.1589.1336.3397618.2322.4687.7238.38184536.7724.3484.8631.0779815.133.2183.7939.44119324.054.7890.5731.2690320.2521.568622.3865419.9314.0980.8621.4995614.723
33、2.3187.67.786512.5911.1889.7141.0193021.4913.886.3329.6993822.0425.3481.5631.3110027.3520.8481.4534.59102425.8239.664.938.47137431.9149EigenvalueDifferenceProportionCumulative3.249175972.034642910.64980.64981.214533060.962968000.24290.89270.251565070.067433970.05030.94310.184131090.083536290.03680.9
34、7990.100594800.0201 1.0000特征根与各因子的贡献特征根与各因子的贡献50 Factor1Factor2x1-0.760620.55316x20.56898-0.76662x30.891840.25374x40.870660.34618x50.890760.36962没有旋转的因子结构没有旋转的因子结构51Factor1可解释方差Factor2可解释方差2.99754292.1642615各旋转后的共同度各旋转后的共同度0.884540230.911439980.859770610.877894530.9300636952 在这个例子中我们得到了两个因子，第一个因子是社会
35、经济发展水平因子，第二个是计划生育因子。有了因子得分值后，则可以利用因子得分为变量，进行其他的统计分析。 Factor1Factor2x1-0.35310-0.87170 x20.077570.95154x30.891140.25621x40.922040.16655x50.951490.15728 Factor1Factor2x1-0.05897-0.49252x2-0.058050.58056x30.330420.03497x40.35108-0.02506x50.36366-0.03493方差最大旋转后的因子结构方差最大旋转后的因子结构标准化得分函数标准化得分函数53 6 因子分析的步骤
36、、展望和建议 计算所选原始变量的相关系数矩阵计算所选原始变量的相关系数矩阵 相关系数矩阵描述了原始变量之间的相关关系。可以帮助判断原始变量之间是否存在相关关系，这对因子分析是非常重要的，因为如果所选变量之间无关系，做因子分析是不恰当的。并且相关系数矩阵是估计因子结构的基础。 选择分析的变量选择分析的变量 用定性分析和定量分析的方法选择变量，因子分析的前提条件是观测变量间有较强的相关性，因为如果变量之间无相关性或相关性较小的话，他们不会有共享因子,所以原始变量间应该有较强的相关性。一、 因子分析通常包括以下五个步骤54 提取公共因子提取公共因子 这一步要确定因子求解的方法和因子的个数。需要根据研
37、究者的设计方案或有关的经验或知识事先确定。因子个数的确定可以根据因子方差的大小。只取方差大于1(或特征值大于1)的那些因子，因为方差小于1的因子其贡献可能很小；按照因子的累计方差贡献率来确定，一般认为要达到60才能符合要求； 因子旋转因子旋转 通过坐标变换使每个原始变量在尽可能少的因子之间有密切的关系，这样因子解的实际意义更容易解释,并为每个潜在因子赋予有实际意义的名字。 55 计算因子得分计算因子得分 求出各样本的因子得分，有了因子得分值，则可以在许多分析中使用这些因子，例如以因子的得分做聚类分析的变量，做回归分析中的回归因子。 56 因子分析是十分主观的，在许多出版的资料中，因子分析模型都
38、用少数可阐述因子提供了合理解释。实际上，绝大多数因子分析并没有产生如此明确的结果。不幸的是，评价因子分析质量的法则尚未很好量化，质量问题只好依赖一个“哇！”准则 如果在仔细检查因子分析的时候，研究人员能够喊出“哇，我明白这些因子”的时候，就可看着是成功运用了因子分析方法。57因子分析和主成分分析的一些注意事项因子分析和主成分分析的一些注意事项 可以看出，因子分析和主成分分析都依赖于原始变量，可以看出，因子分析和主成分分析都依赖于原始变量，也只能反映原始变量的信息。所以原始变量的选择很重也只能反映原始变量的信息。所以原始变量的选择很重要。要。另外，如果原始变量都本质上独立，那么降维就可能失另外，
39、如果原始变量都本质上独立，那么降维就可能失败，这是因为很难把很多独立变量用少数综合的变量概败，这是因为很难把很多独立变量用少数综合的变量概括。数据越相关，降维效果就越好。括。数据越相关，降维效果就越好。在得到分析的结果时，并不一定会都得到如我们例子那在得到分析的结果时，并不一定会都得到如我们例子那样清楚的结果。这与问题的性质，选取的原始变量以及样清楚的结果。这与问题的性质，选取的原始变量以及数据的质量等都有关系数据的质量等都有关系在用因子得分进行排序时要特别小心，特别是对于敏感在用因子得分进行排序时要特别小心，特别是对于敏感问题。由于原始变量不同，因子的选取不同，排序可以问题。由于原始变量不同
40、，因子的选取不同，排序可以很不一样。很不一样。58SPSSSPSS实现实现( (因子分析与主成分分析因子分析与主成分分析) )拿拿student.savstudent.sav为例，选为例，选AnalyzeAnalyzeData ReductionData ReductionFactorFactor进入进入主对话框；主对话框；把把mathmath、physphys、chemchem、literatliterat、historyhistory、englishenglish选入选入VariablesVariables，然后点击，然后点击ExtractionExtraction，在在MethodMet
41、hod选择一个方法（如果是主成分分析，则选选择一个方法（如果是主成分分析，则选Principal Principal ComponentsComponents），），下面的选项可以随意，比如要画碎石图就选下面的选项可以随意，比如要画碎石图就选ScreeScree plot plot，另外在，另外在ExtractExtract选项可以按照特征值的大小选主成分（或因子），也可选项可以按照特征值的大小选主成分（或因子），也可以选定因子的数目；以选定因子的数目；之后回到主对话框（用之后回到主对话框（用ContinueContinue）。然后点击）。然后点击RotationRotation，再在该，再在
42、该对话框中的对话框中的MethodMethod选择一个旋转方法（如果是主成分分析就选选择一个旋转方法（如果是主成分分析就选NoneNone），），在在DisplayDisplay选选Rotated solutionRotated solution（以输出和旋转有关的结果）和（以输出和旋转有关的结果）和Loading plotLoading plot（以输出载荷图）；之后回到主对话框（用（以输出载荷图）；之后回到主对话框（用ContinueContinue）。）。如果要计算因子得分就要点击如果要计算因子得分就要点击ScoresScores，再选择，再选择Save as Save as varia
43、blesvariables（因子得分就会作为变量存在数据中的附加列上）和（因子得分就会作为变量存在数据中的附加列上）和计算因子得分的方法（比如计算因子得分的方法（比如RegressionRegression）；之后回到主对话框）；之后回到主对话框（用（用ContinueContinue）。这时点）。这时点OKOK即可。即可。59主成分分析：洛衫矶对主成分分析：洛衫矶对12个人口调查区的数据个人口调查区的数据编号编号 总人口总人口 总雇员数总雇员数 中等校中等校 专业服务专业服务 中等房价中等房价 平均校龄平均校龄 项目数项目数 1570012.8250027025000 2100010.960
44、01010000 334008.81000109000 4380013.6170014025000 5400012.8160014025000 682008.326006012000 7120011.44001016000 8910011.533006014000 9990012.534001801800010960013.73600390250001196009.63300801200012940011.440001001300060特征值、累积贡献率特征值、累积贡献率Total Variance Explained2.87357.46657.4662.87357.46657.4661.79
45、735.93393.3991.79735.93393.399.2154.29797.6969.993E-021.99999.6951.526E-02.305100.000Component12345Total% ofVarianceCumulative%Total% ofVarianceCumulative%Initial EigenvaluesExtraction Sums of Squared LoadingsExtraction Method: Principal Component Analysis.61Scree PlotComponent Number54321Eigenvalue
46、3.53.02.52.01.51.0.50.0特征值图特征值图62Component PlotComponent 11.0.50.0-.5-1.0Component 21.0.50.0-.5-1.0中等房价专业服务项目数总雇员数中等校平均校龄总人口二主成分因二主成分因子负荷图子负荷图( ,)ijiijr Y Xa63Component Matrixa.932-.104.791-.558.767-.545.581.806.672.726专业服务项目数中等房价中等校平均校龄总人口总雇员数12ComponentExtraction Method: Principal Component Analys
47、is.2 components extracted.a. 主成分的因子负荷主成分的因子负荷(每列平方和为相应特征值每列平方和为相应特征值, 而每列除以而每列除以相应特征值的平方根为相应的特征向量相应特征值的平方根为相应的特征向量)这是主成分与各这是主成分与各个变量的相关系数个变量的相关系数有的书把它当成特征向量了有的书把它当成特征向量了SPSS没有给出特征向量没有给出特征向量( ,)ijiijr Y Xa64销售人员数据销售人员数据（50个观测值）个观测值）销售增长销售增长 销售利润销售利润 新客户销售额新客户销售额 创造力创造力 机械推理机械推理 抽象推理抽象推理 数学推理数学推理93.00
48、96.0097.809.0012.009.0020.0088.8091.8096.807.0010.0010.0015.0095.00100.3099.008.0012.009.0026.00101.30103.80106.8013.0014.0012.0029.00102.00107.80103.0010.0015.0012.0032.0095.8097.5099.3010.0014.0011.0021.0095.5099.5099.009.0012.009.0025.00110.80122.00115.3018.0020.0015.0051.00102.80108.30103.8010.0
49、017.0013.0031.00106.80120.50102.0014.0018.0011.0039.00103.30109.80104.0012.0017.0012.0032.0099.50111.80100.3010.0018.008.0031.00103.50112.50107.0016.0017.0011.0034.0099.50105.50102.308.0010.0011.0034.0065特征值、累积贡献率特征值、累积贡献率Total Variance Explained5.03571.92371.9235.03571.92371.923.93413.33685.259.934
50、13.33685.259.4987.11392.372.4216.01898.3908.104E-021.15899.5472.034E-02.29199.8381.134E-02.162100.000Component1234567Total% ofVarianceCumulative%Total% ofVarianceCumulative%Initial EigenvaluesExtraction Sums of Squared LoadingsExtraction Method: Principal Component Analysis.66Scree PlotComponent Num
