第6章 主成分分析实验报告

  • 格式:doc
  • 大小:155.00 KB
  • 文档页数:9

课程实验报告专业年课程名称应用多元统计分析指导教师学生姓名学号实验日期实验地点实验成绩教务处制年月日实验项目名称因子分析的上机实现实验目的及要求SPSS软件中factor analysis的计算机操作及结果分析,使学生能熟练应用计算机软件进行因子分析与结果分析,培养实际应用能力。

实验内容题目:对企业经济效益体系的8项指标建立因子分析模型(数据7-1)。

这8项指标分别为:x1-固定资产利税率,x2-资金利税率,x3-销售收入利税率,x4-资金利润率,x5-固定资产利润率,x6-资金周转天数,x7-万元产值能耗,x8-全员劳动生产率。

在分析过程中,提取因子的方法为“主成分”法,并以数据的“相关阵”为分析矩阵,并且提取3个因子,采用“最大方差旋转法”进行因子旋转。

(1)则这3个因子的累积方差贡献率为多少?(2)请写出原始变量x1和x2的因子表达式;(3)所提取的3个公共因子分别在8个指标中的哪些指标上有较大载荷?并据此说明所提取的公因子概括了企业的何种能力?(4)分别写出因子得分表达式,并计算“大同”企业的综合因子得分。

注:均以因子旋转后的输出结果进行分析。

要求:1、将SPSS软件的分析过程的关键步骤截图说明,需要计算的地方要写出详细计算步骤。

实验步骤1.选择菜单项Analyze→Data Reduction→Factor。

打开Factor Analysis对话框,将原始变量X1到X8移入V ariables列表框中,如图一所示。

如果不想使用全部的样本进行分析,且数据文件中存在一个选择变量的话,将该选择变量移入Selection Variable 框中,并单击右边的Value按钮,在跳出的窗口中输入一个筛选值,这样,只有选择变量的值等于输入的筛选值的case才能参与因子分析。

图一:2.点击Descriptives按钮,展开相应的对话框。

选择Initial solution复选项。

这个选项给出个因子的特征值,各因子特征值占总方差的百分比以及累计百分比。

单击Continue按钮,返回主界面,如图二所示。

图二:3.单击点击Extraction按钮,打开Extraction子对话框,设置有关因子提取的选项,选择“主成分”因子分子方法,如图三所示。

在Method下拉列表中选择因子提取的方法,SPSS提供了七种提取方法可供选择,一般选择默认选项Principal components,即“主成分法”。

在Analyze选项栏中指定用于提取因子的分析矩阵,分别为相关系数矩阵(Correlation matrix)和协方差矩阵(Covariance matrix)。

如果选择相关系数矩阵,则表示首先对原始数据进行标准化,然后再进行因子分析;如果选择协方差矩阵,则表示直接对原始数据进行因子分析。

这里我们选择默认的相关系数矩阵。

在Display选项栏中指定与因子提取有关的输出项,其中,Unrotated factor solutions表示输出旋转前的因子方差贡献表和旋转前的因子载荷阵;Scree Plot表示输出因子碎石图。

因子碎石图其实就是样本协差阵的特征根按大小顺序排列的折线图,可以用来帮助确定提取多少个因子。

典型的碎石图会有一个明显的拐点,拐点之前是较大特征根连接形成的陡峭折线,拐点之后是较小特征根连接形成的平缓折线,一般选择拐点之前的特征根数目为提取因子的数目。

这里我们将两个选项都选中。

在Extract选项栏中指定因子提取的数目,有两种设置方法:一种是在Eigenvalues over后的输入框中设置提取的因子对应的特征值的范围,系统默认值为1,即要求提取那些特征值大于1的因子;第二种设置方法是直接在Number of factors后的输入框中输入要求提取的公因子的数目。

这里选择第二种,提取公因子数为3个。

4.点击Rotation按钮,打开Rotation子对话框,设置有关因子旋转的选项,选择Varimax(最大方差旋转法),如图四所示Method选项栏用于设置因子旋转的方法,可供选择的方法包括方差最大旋转法(Varimax)、直接斜交旋转法(Direct Oblimin)、四次方最大正交旋转法(Quartmax)、平均正交旋转法(Equamax)、斜交旋转法(Promax),如果选择None选项,则不进行旋转。

Display选项栏用于设置与因子旋转有关的输出项。

其中,Rotated factor solutions 表示输出旋转后的因子方差贡献表和旋转后的因子载荷阵;Loading plots表示输出旋转后的因子载荷散点图图,旋转后因子散点图是以因子为坐标轴,以旋转后因子载荷为坐标的散点图,从该散点图中可以直观地观察因子载荷在各因子上的分布状况。

图四:5.点击Scores按钮,打开Factor Scores子对话框,设置有关因子得分的选项。

如图五所示:选中Save as variables复选框,表示将因子得分作为新变量保存在数据文件中。

提取了几个因子则会在数据文件中保存几个因子得分变量,变量名为“fac m_n”,其中,m 表示第m个因子,n表示进行第n次因子分析的结果。

选中Display factor score coefficient matrix复选框,这样在结果输出窗口中会给出因子得分系数矩阵。

图五:6.单击OK按钮,运行因子分析过程。

实验结果与分析结论:CommunalitiesInitial ExtractionX1 1.000 .976X2 1.000 .968X3 1.000 .862X4 1.000 .986X5 1.000 .672X6 1.000 .906X7 1.000 .773 X8 1.000.824Extraction Method: Principal Component Analysis.表1给出了8个原始变量的变量共同度。

变量共同度反映每个变量对提取出的所有公共因子的依赖程度。

表1可以看出,除X3,X5的共同度小于80%,X6共同度为81.7%外,其余的变量共同度都在90%以上,说明提取的因子已经包含了原始变量的大部分信息,因子提取的效果比较理想。

图2:Undefined error #60601 - Cannot open text file "D:\PROGRA~1\Undefined error #60618 - Cannot open text file "D:\PROGRA~1\SPSS\en\wind87654321654321图2给出了因子的碎石图。

图中横坐标为因子的序号,纵坐标为相应特征根的值。

从图中可以看到,前3个因子的特征根普遍较高,连接成了陡峭的折线,而第4个因子之后的特征根普遍较低,连接成了平缓的折线,这进一步说明提取3因子是比较适当的。

表3旋转后的因子载荷Rotated Component Matrix(a)Component1 2 3 X1 .815 .551 -.089 X2 .974 .107 -.092 X3 .675 .636 -.039 X4 .971 .200 -.057 X5.660.225-.431X6 .032 -.376 .873X7 -.608 .138 .620X8 .142 .857 -.264表3给出了旋转后的因子载荷矩阵,根据该表可以写出每个原始变量的因子表达式:X1=0.815F1+0.551F2+(-0.89)F3X2=0.947F1+0.107F2+(-0.092)F3表7-3给出了因子得分系数矩阵,根据表中的因子得分系数和原始变量的标准化值就可以计算每个观测值的各因子的得分。

本例中旋转后的因子得分表达式可以写成:efficiencyFuelcapacityFuelweightCurbLengthW idth W heelbasHorsepowersizeEnginethousandsinicetypeVehicleF107.0012.0070.0105.0011.0177 .0368.0226.0Pr414.0173.01-++-+-+++-=…………由于我们在Factor Scores子对话框中选择了Save as variables复选框,所以,在数据文件中会生成3个因子得分变量,变量名分别为:fac1_1、fac2_1、fac3_1。

这里有两点值得注意的地方:(1)由于我们是以相关系数矩阵为出发点进行因子分析,所以,因子得分表达式中的各变量应该是经过标准化变换后的标准变量,均值为0,标准差为1。

(2)由于因子载荷阵经过了旋转,所以,因子得分不是利用初始的因子载荷阵,而是利用旋转后的因子载荷阵计算得到的。

Component Score Coefficient MatrixComponent1 2 3X1 .159 .259 .160X2 .331 -.173 .068X3 .091 .381 .208X4 .311 -.082 .122X5 .139 -.075 -.248X6 .218 -.123 .695X7 -.187 .428 .515X8 -.195 .654 -.007Undefined error #11401 - Cannot open text file "D:\PROGRA~1\SPSS\en\windows\spss.err": No such file Undefined error #11408 - Cannot open text file"D:\PROGRA~1\SPSS\en\windows\spss.err": No such file Undefined error #11418 - Cannot open text file "D:\PROGRA~1\SPSS\en\windows\spss.err": No such fileComponent Score Covariance MatrixComponent 1 2 31 1.000 -1.222E-16 .0002 -1.222E-16 1.000 .0003 .000 .000 1.000Undefined error #11401 - Cannot open text file "D:\PROGRA~1\SPSS\en\windows\spss.err": No such file注:可根据实际情况加附表x1 x2 x3 x4 x5 x6 x7 x8 琉璃河16.68 26.75 31.84 18.4 53.25 55 28.83 1.75 邯郸19.70 27.56 32.94 19.2 59.82 55 32.92 2.87 大同15.20 23.40 32.98 16.2 46.78 65 41.69 1.53 哈尔滨7.29 8.97 21.30 4.8 34.39 62 39.28 1.63 华新29.45 56.49 40.74 43.7 75.32 69 26.68 2.14 湘乡32.93 42.78 47.98 33.9 66.46 50 32.87 2.60 柳州25.39 37.82 36.76 27.6 68.18 63 35.79 2.43 峨嵋15.05 19.49 27.21 14.2 6.13 76 35.76 1.75 耀县19.82 28.78 33.41 20.2 59.25 71 39.13 1.83 永登21.13 35.20 39.16 26.5 52.47 62 35.08 1.73 工源16.75 28.72 29.62 19.2 55.76 58 30.08 1.52 抚顺15.83 28.03 26.40 17.4 61.19 61 32.75 1.60 大连16.53 29.73 32.49 20.6 50.41 69 37.57 1.31 江南22.24 54.59 31.05 37.0 67.95 63 32.33 1.57 江油12.92 20.82 25.12 12.5 51.07 66 39.18 1.83。