1、第三讲第三讲 因子分析因子分析Factor Analysis目目 录录1 1 引言引言2 2 因子分析模型因子分析模型 3 3 因子载荷矩阵的估计方法因子载荷矩阵的估计方法4 4 因子旋转(正交变换)因子旋转(正交变换)5 5 因子得分因子得分 6 6 因子分析的因子分析的SPSSSPSS操作操作 因子分析因子分析(factor analysis)(factor analysis)是一种数据简化的技是一种数据简化的技术。它通过研究众多变量之间的内部依赖关系,探求术。它通过研究众多变量之间的内部依赖关系,探求观测数据中的基本结构,并用少数几个假想变量来表观测数据中的基本结构,并用少数几个假想变量
2、来表示其基本的数据结构。这几个假想变量能够反映原来示其基本的数据结构。这几个假想变量能够反映原来众多变量的主要信息。原始的变量是可观测的显变量众多变量的主要信息。原始的变量是可观测的显变量(Observed Variable),而假想变量是不可观测的潜变,而假想变量是不可观测的潜变量量(Latent Variable),称为因子。,称为因子。例如,在企业形象或品牌形象的研究中,消费者例如,在企业形象或品牌形象的研究中,消费者可以通过一个有可以通过一个有2424个指标构成的评价体系,评价百货个指标构成的评价体系,评价百货商场的商场的2424个方面的优劣。个方面的优劣。1 引引 言言 但消费者主要
3、关心的是三个方面,即商店的环境但消费者主要关心的是三个方面,即商店的环境、商店的服务和商品的价格。因子分析方法可以通过、商店的服务和商品的价格。因子分析方法可以通过2424个变量,找出反映商店环境、商店服务水平和商品价格个变量,找出反映商店环境、商店服务水平和商品价格的三个潜在的因子,对商店进行综合评价。的三个潜在的因子,对商店进行综合评价。而这三个公而这三个公共因子可以表示为:共因子可以表示为:iiiiiiFFFx33221124,1i 称称 是不可观测的潜在因子。是不可观测的潜在因子。2424个变量个变量共享这三个因子,但是每个变量又有自己的个性,共享这三个因子,但是每个变量又有自己的个性
4、,即不被即不被包含的部分包含的部分 ,称为特殊因子。,称为特殊因子。321FFF、i注:注:因子分析与回归分析不同,因子分析中的因因子分析与回归分析不同,因子分析中的因子是一个比较抽象的概念,而回归因子有非常明子是一个比较抽象的概念,而回归因子有非常明确的实际意义;确的实际意义;主成分分析分析与因子分析也有不同,主成主成分分析分析与因子分析也有不同,主成分分析仅仅是变量变换,而因子分析需要构造因分分析仅仅是变量变换,而因子分析需要构造因子模型。子模型。主成分分析主成分分析:原始变量的线性组合表示新的原始变量的线性组合表示新的综合变量,即主成分;综合变量,即主成分;因子分析:潜在的假想变量和随机
5、影响变因子分析:潜在的假想变量和随机影响变量的线性组合表示原始变量。量的线性组合表示原始变量。2 因子分析模型 一、数学模型一、数学模型 设设 个变量,如果表示为个变量,如果表示为iX),2,1(pip11iiiimmiXa Fa F)(pm 11111211122212222212mmpppppmpmXFXFXF或XAF或 称为称为 公共因子,是不可观测的变量,公共因子,是不可观测的变量,他们的系数称为因子载荷。他们的系数称为因子载荷。是特殊因子,是不能被是特殊因子,是不能被前前m个公共因子包含的部分。并且满足:个公共因子包含的部分。并且满足:mFFF,21iIFD111)(cov(,)0,
6、F,F不相关;不相关;mFFF,21即即 互不相关,且各自的方差为互不相关,且各自的方差为1 1。22221)(pD即各即各 之间互不相关,且各自的方差不一定相等之间互不相关,且各自的方差不一定相等。),0(2iiNi用矩阵的表达方式:用矩阵的表达方式:X-=AF+()EF0()E0()VarFI22212()(,)pVardiag1 11212 122212()()()()()()cov()()()()()ppppppE FE FE FE FE FE FEE FE FE FF,F0二、因子分析模型的性质二、因子分析模型的性质 1、原始变量、原始变量X的协方差矩阵的分解的协方差矩阵的分解X-=
7、AF+()()()VarVarVarX-=AF A+x=AA+DA是因子模型的系数22212()(,)pVardiagD D的主对角线上的元素值越小,则公共因子共享的成的主对角线上的元素值越小,则公共因子共享的成分越多。分越多。2、因子载荷不是惟一的、因子载荷不是惟一的 设设T为一个为一个pp的正交矩阵,令的正交矩阵,令A*=AT,F*=TF,则模型可以表示为则模型可以表示为*X+A F+()ET F0()E0*()()()VarVarVarFTFTF TI22212()(,)pVardiag*cov()()EF,F 0且满足条件因子模型的条件且满足条件因子模型的条件 三、三、因子载荷矩阵中的
8、几个统计特征因子载荷矩阵中的几个统计特征 1 1、因子载荷、因子载荷aij的统计意义的统计意义 因子载荷因子载荷 是第是第i个个标准化变量标准化变量与第与第j个公共因子的相关系数个公共因子的相关系数 ija模型为模型为 imimiiFaFaX11 在上式的左右两边乘以在上式的左右两边乘以 jF,再求数学期望再求数学期望)()()()()(11jijmimjjijjijiFEFFEaFFEFFEaFXE 根据公共因子的模型性质,有根据公共因子的模型性质,有ijFxji (载荷矩阵中第(载荷矩阵中第i行,第行,第j列的元素)反映了第列的元素)反映了第i个变量与第个变量与第j个公共因子的相关重要性。
9、绝对值越大,相个公共因子的相关重要性。绝对值越大,相关的密切程度越高。注意关的密切程度越高。注意标准化变量的方差为标准化变量的方差为1 1。2 2、变量共同度的统计意义、变量共同度的统计意义定义:变量定义:变量 的共同度是因子载荷矩阵的第的共同度是因子载荷矩阵的第i行的元行的元素的平方和。记为素的平方和。记为iX统计意义统计意义:imimiiFaFaX11两边求方差:两边求方差:)()()()(2112imimiiVarFVaraFVaraXVarmjiija1221 所有的公共因子和特殊因子对变量所有的公共因子和特殊因子对变量 的贡献为的贡献为1 1。如果。如果 非常非常靠近靠近1 1,非常
10、小,则因子分析的效果好,从原变量空间到公共因非常小,则因子分析的效果好,从原变量空间到公共因子空间的转化性质好。子空间的转化性质好。iXmjija122imjija12。mjijiah122iX(注意注意 为标准化变量为标准化变量)3 3、公共因子、公共因子 方差贡献的统计意义方差贡献的统计意义jF因子载荷矩阵中各列元素的平方和因子载荷矩阵中各列元素的平方和 称为第称为第j个公共因子个公共因子 对对 X 的所有分量的所有分量 的的方差贡献和,它衡量了第方差贡献和,它衡量了第j个公共因子个公共因子 在全体公共因子中的在全体公共因子中的相对重要性。相对重要性。piijjaS12jF),1(piiX
11、jF3 3 因子载荷矩阵的估计方法因子载荷矩阵的估计方法 设随机向量设随机向量 的均值为的均值为,协方差为,协方差为,为为 的特征根,的特征根,为对应的为对应的标准化特征向量,则标准化特征向量,则pxxx,21x021pp21u,u,u12p=UUAA+D(一)主成分分析法(一)主成分分析法 上式给出的上式给出的表达式是精确的,然而,它实际上是毫表达式是精确的,然而,它实际上是毫无价值的,因为我们的目的是寻求用少数几个公共因子无价值的,因为我们的目的是寻求用少数几个公共因子解释,故略去后面的解释,故略去后面的p-m项的贡献,有项的贡献,有21111mmmmmmp1122ppu uu uu uu
12、uu up2uuuuuuppp21122111100p212ppuuuuuu 12 mmm1122AA+Du uu uu uD1121122 mmp mpmm p2uuuuuDAADu22212(,)pdiagD其中221miiiijjsampijmmax2211)(uuuA),1(pi(二)主因子法(二)主因子法 主因子方法是对主成分方法的修正,假定我们首先主因子方法是对主成分方法的修正,假定我们首先对变量进行标准化变换。设对变量进行标准化变换。设 R=AA+D R*=AA=R-D称称R*为约相关矩阵。为约相关矩阵。R*对角线上的元素是对角线上的元素是 ,而不,而不是是1。2ih211212
13、2122212ppppphrrrhrRrrhR-D直接求直接求R*的前的前p个特征根和对应的正交特征向量。得如下个特征根和对应的正交特征向量。得如下的矩阵:的矩阵:*1122ppAuuu*10pR特征根:*12,pu uu正交特征向量:21222pRR 当特殊因子当特殊因子 的方差不大且的方差不大且已知的,问题非常好解决。已知的,问题非常好解决。i*11*221122*ppppuuuuuu*1122mmAuuu2121100phhD 在实际的应用中,特殊因子方差矩阵一般都是未知在实际的应用中,特殊因子方差矩阵一般都是未知的,可以通过一组样本来估计。的,可以通过一组样本来估计。估计的估计的方法有
14、如下几方法有如下几种种:首先,求首先,求 的初始估计值,构造出的初始估计值,构造出 2ih*R 1)取取 ,在这个情况下主因子解与主成分解等,在这个情况下主因子解与主成分解等价;价;2 2)取)取 ,为为xi与其他所有的原始变量与其他所有的原始变量xj的复的复相关系数的平方,即相关系数的平方,即xi对其余的对其余的p-1-1个个xj的回归方程的的回归方程的判定系数,这是因为判定系数,这是因为xi 与公共因子的关系是通过其余与公共因子的关系是通过其余的的p-1-1个个xj 的线性组合联系起来的;的线性组合联系起来的;12ih22iiRh 2iR 3 3)取)取 ,这意味着取,这意味着取xi与其余
15、的与其余的xj的的简单相关系数的绝对值最大者;简单相关系数的绝对值最大者;)(|max2ijrhiji 4 4)取)取 ,其中要求该值为正数。,其中要求该值为正数。pjijijirph,1211 5 5)取)取 ,其中,其中 是是 的对角元素。的对角元素。iiirh/112iir1R例例 假定某地固定资产投资率假定某地固定资产投资率 ,通货膨胀率,通货膨胀率 ,失,失业率业率 ,相关系数矩阵为,相关系数矩阵为试用主成分分析法求因子分析模型。试用主成分分析法求因子分析模型。1x2x3x15/25/15/215/15/15/11 特征根为特征根为:55.11 85.02 6.03 6.0707.0
16、85.0331.055.1629.06.0707.085.0331.055.1629.0085.0883.055.1475.0A707.0331.0629.0707.0331.0629.00883.0475.0U548.0305.0783.0548.0305.0783.00814.0569.0 可取前两个因子可取前两个因子F1和和F2为公共因子,第一公因为公共因子,第一公因子子F1物价就业因子,对物价就业因子,对X的贡献为的贡献为1.55。第二公因子。第二公因子F2为投资因子,对为投资因子,对X的贡献为的贡献为0.85。共同度分别为。共同度分别为1,0.706,0.706。211814.056
17、9.0FFx3212548.0305.0783.0FFFx3213548.0305.0783.0FFFx 4 4 因子旋转(正交变换)因子旋转(正交变换)建立了因子分析数学目的不仅仅要找出公共因子建立了因子分析数学目的不仅仅要找出公共因子以及对变量进行分组,更重要的要知道每个公共因子以及对变量进行分组,更重要的要知道每个公共因子的意义,以便进行进一步的分析,如果每个公共因子的意义,以便进行进一步的分析,如果每个公共因子的含义不清,则不便于进行实际背景的解释。由于因的含义不清,则不便于进行实际背景的解释。由于因子载荷阵是不惟一的,所以应该对因子载荷阵进行旋子载荷阵是不惟一的,所以应该对因子载荷阵
18、进行旋转。目的是使因子载荷阵的结构简化,使载荷矩阵每转。目的是使因子载荷阵的结构简化,使载荷矩阵每列或行的元素平方值向列或行的元素平方值向0 0和和1 1两极分化。有三种主要的两极分化。有三种主要的正交旋转法。四次方最大法、正交旋转法。四次方最大法、方差最大法方差最大法和等量最大和等量最大法。法。(一)为什么要旋转因子(一)为什么要旋转因子 百米跑成绩百米跑成绩 跳远成绩跳远成绩 铅球成绩铅球成绩 跳高成绩跳高成绩 400米跑成绩米跑成绩 百米跨栏百米跨栏 铁饼成绩铁饼成绩 撑杆跳远成绩撑杆跳远成绩 标枪成绩标枪成绩 1500米跑成绩米跑成绩 1X2X3X4X5X6X7X8X9X10X奥运会十
19、项全能运动项目奥运会十项全能运动项目得分数据的因子分析得分数据的因子分析 102.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.01变量共同度0.6910.217-0.58-0.2060.840.7890.184-0.1930.0920.70.7020.5350.047-0.1
20、750.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 因子载荷矩阵可以看出,除第一因子在所有的变量在公共因子因子载荷矩阵可以看出,除第一因子在所有的变量在公共因子上有较大的正载荷,可以称为一般运动因子。其他的上有较大的正载荷,可以称为一
21、般运动因子。其他的3个因子不太个因子不太容易解释。似乎是跑和投掷的能力对比,似乎是长跑耐力和短跑速容易解释。似乎是跑和投掷的能力对比,似乎是长跑耐力和短跑速度的对比。于是考虑旋转因子,得下表度的对比。于是考虑旋转因子,得下表 变量共同度0.844*0.1360.156-0.1130.840.631*0.1940.515*-0.0060.70.2430.825*0.223-0.1480.810.2390.150.750*0.0760.650.797*0.0750.1020.4680.870.4040.1530.635*-0.170.620.1860.814*0.147-0.0790.72-0.0
22、360.1760.762*0.2170.66-0.0480.735*0.110.1410.570.045-0.0410.1120.934*0.891F2F3F4F1X2X3X4X5X6X7X8X9X10X 通过旋转,因子有了较为明确的含义。通过旋转,因子有了较为明确的含义。百米跑,百米跑,跳远和跳远和 400米跑,需要爆发力的项目在米跑,需要爆发力的项目在 有较大的有较大的载荷,载荷,可以称为短跑速度因子;可以称为短跑速度因子;铅球,铅球,铁饼和铁饼和 标枪在标枪在 上有较大的载荷,可以上有较大的载荷,可以称为爆发性臂力因子;称为爆发性臂力因子;百米跨栏,百米跨栏,撑杆跳远,撑杆跳远,跳远和为
23、跳远和为 跳高在跳高在 上上有较大的载荷,有较大的载荷,爆发腿力因子;爆发腿力因子;长跑耐力因子。长跑耐力因子。2X5X1F1F3X7X9X2F6X8X2X4X3F3F4F1X 方差最大法方差最大法 方差最大法从简化因子载荷矩阵的每一列出发,使和每个因方差最大法从简化因子载荷矩阵的每一列出发,使和每个因子有关的载荷的平方的方差最大。当只有少数几个变量在某个子有关的载荷的平方的方差最大。当只有少数几个变量在某个因子上又较高的载荷时,对因子的解释最简单。因子上又较高的载荷时,对因子的解释最简单。方差最大的直方差最大的直观意义是希望通过因子旋转后,使每个因子上的载荷尽量拉开观意义是希望通过因子旋转后
24、,使每个因子上的载荷尽量拉开距离,一部分的载荷趋于距离,一部分的载荷趋于 1 1,另一部分趋于,另一部分趋于0 0。2122211211ppaaaaaaA221122212122121111FaFaXFaFaXFaFaXppp(二)旋转方法(二)旋转方法cossinsincosT设旋转矩阵为:cossinsincosAATB则cossinsincoscossinsincos112112111211ppppaaaaaaaa*2*1*12*11ppaaaa1,2,;1,2ijijiadip jh令211(pjijiddp这是列和)max)()(1212 mjpijijddV简化准则为:00V令,则
25、可以解出0000cossinsincosT旋转矩阵为:max(8.4.2)123m即:V+V+V+V1000cossin0sincosT1000cossin0sincos T111 TT变换后因子的共同度变换后因子的共同度设设 正交矩阵,做正交变换正交矩阵,做正交变换 AB)()(1mlljilppijabBmjmjmlljilijiabh111222)()(B mjmlmjmlmljttjljitilljilaaa1 11 1 122)(2111222Aimlmjmlilljilhaa变换后因子的共同度没有发生变化!变换后因子的共同度没有发生变化!(三)旋转结果(三)旋转结果变换后因子贡献变
26、换后因子贡献设设 正交矩阵,做正交变换正交矩阵,做正交变换 AB)()(1qlljilppijabBpipiqlljilijjabS111222)()(B piqlpiqlqltttjljitilljilaaa1111 122piqlqlljjljilSa1112222)(A变换后因子的贡献发生了变化!变换后因子的贡献发生了变化!5 5 因子得分因子得分(一)因子得分的概念(一)因子得分的概念 前面我们主要解决了用公共因子的线性组合来表示前面我们主要解决了用公共因子的线性组合来表示一组观测变量的有关问题。如果我们要使用这些因子做一组观测变量的有关问题。如果我们要使用这些因子做其他的研究,比如把
27、得到的因子作为自变量来做回归分其他的研究,比如把得到的因子作为自变量来做回归分析,对模型进行诊断析,对模型进行诊断,进一步分析原始数据进一步分析原始数据(如如对样本进对样本进行分类或评价行分类或评价),这就需要我们对公共因子进行测度,即这就需要我们对公共因子进行测度,即给出公共因子的值。给出公共因子的值。例例1 1 人均要素变量因子分析人均要素变量因子分析。对我国32个省市自治区的要素状况作因子分析。指标体系中有如下指标:X1:人口(万人)X2:面积(万平方公里)X3:GDP(亿元)X4:人均水资源(立方米/人)X5:人均生物量(吨/人)X6:万人拥有的大学生数(人)X7:万人拥有科学家、工程
28、师数(人)Rotated Factor Pattern FACTOR1 FACTOR2 FACTOR3 X1 -0.21522 -0.27397 0.89092 X2 0.63973 -0.28739 -0.28755 X3 -0.15791 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.07246 高载荷指标 因子命名 因子1X2;面积(万平方公里)X4:人均水资源(立方米/人)X5:人均
29、生物量(吨/人)自然资源因子 因子2X6:万人拥有的大学生数(人)X7:万人拥有的科学家、工程师数(人)人力资源因子 因子3 X1;人口(万人)X3:GDP(亿元)经济发展总量因子 X1=-0.21522F1-0.27397F2+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.11
30、041F1+0.97851F2-0.07246F3 Standardized Scoring Coefficients FACTOR1 FACTOR2 FACTOR3 X1 0.05764 -0.06098 0.50391 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.
31、14635X3+0.47920X4+0.45583X5+0.05416X6+0.05790X7F1=0.05764X1+0.22724X2+0.14635X3+0.47920X4+0.45583X5+0.05416X6+0.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.
32、17062X4+0.10129X5+0.04099X6+0.04822X7F3=0.50391X1-0.07713X2+0.59715X3+0.17062X4+0.10129X5+0.04099X6+0.04822X7REGION 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.205
33、110.463770.34087jilin-0.214990.10608-0.57431heilongj 0.10839-0.11717-0.02219shanghai-0.201992.38962-0.04259前三个因子得分前三个因子得分例例2 国民生活质量的因素分析国民生活质量的因素分析 国家发展的最终目标,是为了全面提高全体国民的生国家发展的最终目标,是为了全面提高全体国民的生活质量,满足广大国民日益增长的物质和文化的合理需求。活质量,满足广大国民日益增长的物质和文化的合理需求。在可持续发展消费的统一理念下,增加社会财富,创自更在可持续发展消费的统一理念下,增加社会财富,创自更多的物质
34、文明和精神文明,保持人类的健康延续和生生不多的物质文明和精神文明,保持人类的健康延续和生生不息,在人类与自然协同进化的基础上,维系人类与自然的息,在人类与自然协同进化的基础上,维系人类与自然的平衡,达到完整的代际公平和区际公平平衡,达到完整的代际公平和区际公平(即时间过程的最即时间过程的最大合理性与空间分布的最大合理化大合理性与空间分布的最大合理化)。从从19901990年开始,联合国开发计划署年开始,联合国开发计划署(UYNP)(UYNP)首次采用首次采用“人文发展系数人文发展系数”指标对于国民生活质量进行测度。人文指标对于国民生活质量进行测度。人文发展系数利用三类内涵丰富的指标组合,即人的
35、健康状况发展系数利用三类内涵丰富的指标组合,即人的健康状况(使用出生时的人均预期寿命表达使用出生时的人均预期寿命表达)、人的智力程度、人的智力程度(使用使用组合的教育成就表达组合的教育成就表达)、人的福利水平、人的福利水平(使用人均国民收入使用人均国民收入或人均或人均GDPGDP表达表达),并且特别强调三类指标组合的整体表达,并且特别强调三类指标组合的整体表达内涵,去衡量一个国家或地区的社会发展总体状况以及国内涵,去衡量一个国家或地区的社会发展总体状况以及国民生活质量的总水平。民生活质量的总水平。在这个指标体系中有如下的指标:在这个指标体系中有如下的指标:X X1 1预期寿命预期寿命X X2
36、2成人识字率成人识字率X X3 3综合入学率综合入学率X X4 4人均人均GDPGDP(美圆)(美圆)X X5 5预期寿命指数预期寿命指数X X6 6教育成就指数教育成就指数X X7 7人均人均GDPGDP指数指数 旋转后的因子结构旋转后的因子结构 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.34100 X5 0.38854 0.43295 0
37、.80848 X6 0.28207 0.85325 0.43289 X7 0.90091 0.20612 0.35052 FACTOR1为经济发展因子为经济发展因子 FACTOR2为教育成就因子为教育成就因子 FACTOR3为健康水平因子为健康水平因子 被每个因子解释的方差和共同度被每个因子解释的方差和共同度 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
38、.987530 0.945796 0.852306 0.975830 0.992050 X6 X7 0.994995 0.976999 Standardized Scoring Coefficients标准化得分系数标准化得分系数 FACTOR1 FACTOR2 FACTOR3 X1 -0.18875 -0.34397 0.85077 X2 -0.24109 0.60335 -0.10234 X3 0.35462 0.50232 -0.59895 X4 0.53990 -0.17336 -0.10355 X5 -0.17918 -0.31604 0.81490 X6 -0.09230 0.62
39、258 -0.24876*6*5*4*3*2*109230.017918.05399.035462.024109.018875.01xxxxxxf*6*5*4*3*2*162258.031604.017336.050232.060335.034397.02xxxxxxf*6*5*4*3*2*124876.081490.010335.059895.010234.085077.03xxxxxxf例例3 生育率的影响因素分析生育率的影响因素分析 生育率受社会、经济、文化、计划生育政策等很多生育率受社会、经济、文化、计划生育政策等很多因素影响,但这些因素对生育率的影响并不是完全独立因素影响,但这些因素
40、对生育率的影响并不是完全独立的,而是交织在一起,如果直接用选定的变量对生育率的,而是交织在一起,如果直接用选定的变量对生育率进行多元回归分析,最终结果往往只能保留两三个变量,进行多元回归分析,最终结果往往只能保留两三个变量,其他变量的信息就损失了。因此,考虑用因子分析的方其他变量的信息就损失了。因此,考虑用因子分析的方法,找出变量间的数据结构,在信息损失最少的情况下法,找出变量间的数据结构,在信息损失最少的情况下用新生成的因子对生育率进行分析。用新生成的因子对生育率进行分析。选择的变量有:多子率、综合节育率、初中以上文化选择的变量有:多子率、综合节育率、初中以上文化程度比例、城镇人口比例、人均
41、国民收入。下表是程度比例、城镇人口比例、人均国民收入。下表是1990年中国年中国30个省、自治区、直辖市的数据。个省、自治区、直辖市的数据。多子率(%)综合节育率(%)初中以上文化程度比例(%)人均国民收入(元)城镇人口比例(%)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
42、.891.4762.36482266.235.9490.3140.85169621.242.692.4235.14171732.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.789
43、0.5731.2690320.2521.568622.3865419.9314.0980.8621.4995614.7232.3187.67.786512.5911.1889.7141.0193021.4913.886.3329.6993822.0425.3481.5631.3110027.3520.8481.4534.59102425.8239.664.938.47137431.91EigenvalueDifferenceProportionCumulative3.249175972.034642910.64980.64981.214533060.962968000.24290.89270.
44、251565070.067433970.05030.94310.184131090.083536290.03680.97990.100594800.0201 1.0000特征根与各因子的贡献特征根与各因子的贡献 Factor1Factor2x1-0.760620.55316x20.56898-0.76662x30.891840.25374x40.870660.34618x50.890760.36962没有旋转的因子结构没有旋转的因子结构各原始变量在旋转后的共同度各原始变量在旋转后的共同度0.884540230.911439980.859770610.877894530.93006369 在这个
45、例子中我们得到了两个因子,第一个因子是社会经济在这个例子中我们得到了两个因子,第一个因子是社会经济发展水平因子,第二个是计划生育因子。有了因子得分值后,则发展水平因子,第二个是计划生育因子。有了因子得分值后,则可以利用因子得分为变量,进行其他的统计分析。可以利用因子得分为变量,进行其他的统计分析。Factor1Factor2x1-0.35310-0.87170 x20.077570.95154x30.891140.25621x40.922040.16655x50.951490.15728 Factor1Factor2x1-0.05897-0.49252x2-0.058050.58056x30.
46、330420.03497x40.35108-0.02506x50.36366-0.03493方差最大旋转后的因子结构方差最大旋转后的因子结构标准化得分函数标准化得分函数6 6 因子分析的因子分析的SPSSSPSS操作操作 原始数据:原始数据:2020口盐泉的化学成分。口盐泉的化学成分。建立建立Excel数据表文件数据表文件(数据存放在数据存放在Sheet 1)打开打开SPSS 16.0;在;在SPSS开始向导窗口,选定开始向导窗口,选定“Create new query using Database Wizard”,点击,点击“OK”按钮,进入按钮,进入“Database Wizard”向导窗
47、口向导窗口在在“Database Wizard”向导窗口,选定向导窗口,选定“Excel Files”,点击,点击“Next”按钮按钮在在“ODBC Driver Login”向导窗口,使用向导窗口,使用“Browse”按钮,搜索选定按钮,搜索选定Excel数据数据源文件,点击源文件,点击“Open”按钮,点击按钮,点击“OK”按钮按钮在在“Database Wizard”向导窗口,在其左边的窗口中选定向导窗口,在其左边的窗口中选定Excel数据表(本例为数据表(本例为Sheet 1),点击位于中间的),点击位于中间的“箭头箭头”按钮,将该数据表(本例为按钮,将该数据表(本例为Sheet 1)
48、移入到)移入到右边的空白窗口之中,点击右边的空白窗口之中,点击“Next”按钮,第二次点击按钮,第二次点击“Next”按钮,第三次点击按钮,第三次点击“Next”按钮按钮在在“Database Wizard”向导窗口,使用向导窗口,使用“Browse”按钮指定文件夹,命名按钮指定文件夹,命名SPSS数据转换文件(数据转换文件(*.spq)(建议采用与)(建议采用与Excel数据文件相同的名称),点击数据文件相同的名称),点击“Save”按按钮钮点击点击“Next”按钮,点击按钮,点击“Finish”按钮按钮自动回到自动回到SPSS 16.0视窗,显示为未命名的视窗,显示为未命名的SPSS数据集
49、编辑窗口数据集编辑窗口在在SPSS数据集编辑窗口的底部,点击数据集编辑窗口的底部,点击“Variable View”,定义各变量的数据类型、,定义各变量的数据类型、宽度、小数位数等宽度、小数位数等在在SPSS数据集编辑窗口的底部,点击数据集编辑窗口的底部,点击“Data View”,显示出最终的数据集,显示出最终的数据集利用该数据集,在利用该数据集,在SPSS主菜单,选定主菜单,选定“AnalyzeData ReductionFactor”点击,点击,得到因子分析引导窗口得到因子分析引导窗口“Factor Analysis”在因子分析引导窗口在因子分析引导窗口“Factor Analysis”
50、,将待分析的原始变量用箭头输入到右,将待分析的原始变量用箭头输入到右边的边的“Variables”窗口窗口在因子分析引导窗口在因子分析引导窗口“Factor Analysis”,点击,点击“Descriptives”按钮,进入按钮,进入“Descriptives”对话框,从中选取需要的统计量,选完后点击对话框,从中选取需要的统计量,选完后点击“Continue”在因子分析引导窗口在因子分析引导窗口“Factor Analysis”,点击,点击“Extraction”按钮,进入按钮,进入“Extraction”对话框,从中选择提取因子的方法,选完后点击对话框,从中选择提取因子的方法,选完后点击“