因子分析报告和主成分分析报告Word文档格式.docx

1、9.7713.744.788.594.991.3410.027.510.162.136.226.144.529.842.172.731.097.577.287.0712.661.792.100.829.037.082.5911.764.541.285.513.981.306.925.337.302.403.270.620.443.367.638.848.398.747.013.3111.683.534.761.129.649.491.0313.5713.1318.522.359.731.331.009.8711.063.702.981.179.177.859.912.627.125.493.6

2、89.722.643.431.194.693.015.982.763.552.011.275.814.575.381.661.611.572.801.782.093.725.905.761.555.407.501.979.271.5113.69.0212.671.754.922.5410.053.965.241.434.944.386.686.499.062.817.232.301.777.794.395.372.279.467.311.0412.0011.5816.182.429.555.354.2511.742.773.511.054.508.071.298.213.089.103.754

3、.661.729.416.445.1112.52.453.100.91实验步骤：1建立数据文件。定义变量名：分别为X1、X2、X3、X4、X5、X6、X7，按顺序输入相应数值，建立数据文件,保存为“生化检验”。2选择菜单“分析降维因子分析” ，弹出“因子分析”对话框。在对话框左侧的变量列表中选变量X1至X7，进入“变量”框，如图1。3单击“描述”按钮，弹出“因子分析：描述统计”对话框，在“统计量”中选“单变量描述性”项，输出各变量的均数与标准差，“在相关矩阵”栏选“系数”，计算相关系数矩阵，并选“KMO 和Bartlett的球型度检验”项，对相关系数矩阵进行统计学检验，如图2。图1图24单击“

4、抽取”按钮，弹出“因子分析：抽取”对话框，选用“主成分”方法提取因子，如图3。图35单击“旋转”按钮，弹出“因子分析：旋转”对话框，在“方法”栏中选择“ 最大方差”进行因子正交旋转，如图4。 6单击“得分”按钮，弹出“因子分析：得分”对话框，选择“回归”项估计因子得分系数，如图5。7单击“确定”，得到输出结果。图4图5实验结果（1）考察原有变量是否适合进行因子分析表1表1显示原有变量的相关系数矩阵，可以看出大部分的相关系数都比较高，各变量呈较强的线性关系，能够从中提取公共因子，适合进行因子分析。表2由表1-2可知，巴特利特球度检验统计量观测值为326.191，相应的概率p值接近0，如果显著性水

5、平为0.05，由于概率p值小于显著性水平，应拒绝零假设，可以认为相关系数矩阵与单位阵有显著差异，同时KMO值为0.321，根据Kaiser给出的KMO度量标准可知KMO值偏小，意味着因子分析的结果可能不能接受。（2）提取因子表3表4表3和表4为因子方差表，提取因子后因子方差的值均很高，表明提取的因子能很好的描述这7个指标。方差分解表也表明前两个因子能够解释7个指标的88.593%。综合以上，提取前两个因子最好了。（3）因子的命名解释表5表6由表5、表6可知，由旋转成份矩阵可以看出，经旋转后，因子便于命名和解释。因子1主要解释的是X1、X2、X4、X7，命名为FAC1_1；而因子2主要解释的是其

6、余三个指标，X3、X5、X6。命名为FAC2_1。九、课外作业：1. 数据文件“development.sav”是某年我国各省发展状况的一些指标，包括人均GDP、人力资源指数CAPITAL、人均收入INCOME、人均净收入NETINC、教育指数、健康指数，试用主成分法或者因子分析法寻找这些指标主要代表了发展状况的哪些特征，以及各省市的发展程度排序。2. 对某市15个大中型工业企业经济效益进行分析。经研究，从有关经济效益指标中选择7个指标作分析，即：固定资产产值率、固定资产利税率、资金利润率、资金利税率、流动资金周转天数、销售收入利税率和全员劳动生产率，数据文件为“某市工业企业效益指标.sav”

7、，试研究该市大中型工业企业经济效益的状况及差异。1、实验步骤：1打开数据文件“development.sav”。在对话框左侧的变量列表中选变量gdp，capital，income，netinc，education，healcare，进入“变量”框，如图1-1。描述统计”对话框，在“统计量”中选“单变量描述性”项，输出各变量的均数与标准差，“在相关矩阵”栏选“系数”，计算相关系数矩阵，并选“KMO 和Bartlett的球型度检验”项，对相关系数矩阵进行统计学检验，如图1-2。图1-1图1-2抽取”对话框，选用“主成分”方法提取因子，如图1-3。图1-3旋转”对话框，在“方法”栏中选择“ 最大方差

8、”进行因子正交旋转，如图1-4。得分”对话框，选择“回归”项估计因子得分系数，如图1-5。7单击“确定”钮，得到输出结果。图1-4图1-52、实验结果表1-1表1-1显示原有变量的相关系数矩阵，可以看出大部分的相关系数都比较高，各变量呈较强的线性关系，能够从中提取公共因子，适合进行因子分析。表1-2由表1-2可知，巴特利特球度检验统计量观测值为148.715，相应的概率p值接近0，如果显著性水平为0.05，由于概率p值小于显著性水平，应拒绝零假设，可以认为相关系数矩阵与单位阵有显著差异，同时KMO值为0.635，根据Kaiser给出的KMO度量标准可知原有变量适合进行因子分析。表1-3表1-4

9、表1-3和表1-4为因子方差表，提取因子后因子方差的值均很高，表明提取的因子能很好的描述这6个指标。方差分解表也表明前两个因子能够解释6个指标的85.255%。表1-5表1-6由表1-5、表1-6可知，由旋转成份矩阵可以看出，经旋转后，因子便于命名和解释。因子1主要解释的是人均GDP、人均收入和人均净收入，可以命名为经济因子；而因子2主要解释的是其余三个指标，人力资源指数、EDUCATION和健康指数。可以命名为民生因子。因子分析要求，最后得到的因子之间没有相关性，而因子转换矩阵显示，两个因子相关。可见，对因子进行旋转是完全有必要的。（4）计算因子得分表1-7根据表1-7可写出以下因子得分函数

10、：F1=0.362人均GDP+0.038人力资源指数+0.332人均收入+0.350人均净收入-0.030EDUCATION-0.151健康指数F2=-0.074人均GDP+0.324人力资源指数-0.059人均收入-0.025人均净收入+0.397EDUCATION+0.445健康指数（5）计算综合得分。单击“转换计算变量”，新建输出变量，命名为“score分”并输入其计算公式，如图1-6所示。然后，新建一个变量“rank”，一次输入1到29。很清楚地可以对29个省进行评价。图1-6图1-7综合以上分析，可以认为指标人均GDP、人均收入、人均净收入代表了发展状况中的经济特征，而指标人力资源指

11、数、健康指数、教育则代表了发展状况中的民生特征。各省市的排序如图1-7所示，即为广西蒙2、实验步骤：1打开数据文件“某市工业企业效益指标.sav”。在对话框左侧的变量列表中选变量x1至x7，进入“变量”框，如图2-1。描述统计”对话框，在“统计量”中选“单变量描述性”项，输出各变量的均数与标准差，“在相关矩阵”栏选“系数”，计算相关系数矩阵，并选“KMO 和Bartlett的球型度检验”项，对相关系数矩阵进行统计学检验，如图2-2。图2-1图2-2抽取”对话框，选用“主成分”方法提取因子，如图2-3。图2-3旋转”对话框，在“方法”栏中选择“ 最大方差”进行因子正交旋转，如图2-4。得分”对话

12、框，选择“回归”项估计因子得分系数，如图2-5。图2-4图2-5表2-1表2-1显示原有变量的相关系数矩阵，可以看出大部分的相关系数都比较高，各变量呈较强的线性关系，能够从中提取公共因子，适合进行因子分析。表2-2由表2-2可知，巴特利特球度检验统计量观测值为136.426，相应的概率p值接近0，如果显著性水平为0.05，由于概率p值小于显著性水平，应拒绝零假设，可以认为相关系数矩阵与单位阵有显著差异，同时KMO值为0.635，根据Kaiser给出的KMO度量标准可知原有变量适合进行因子分析。表2-3表2-4表2-3和表2-4为因子方差表，提取因子后因子方差的值均很高，表明提取的因子能很好的描

13、述这7个指标。方差分解表也表明前两个因子能够解释7个指标的84.619%。表2-5表2-6由表2-5、表2-6可知，由旋转成份矩阵可以看出，经旋转后，因子便于命名和解释。因子1主要解释的是固定资产产值率、固定资产利税率、资金利润率、资金利税率和销售收入利税率，可以命名为盈利能力因子；而因子2主要解释的是其余两个指标，流动资金周转天数、和全员劳动生产率。可以命名为资金和人力因子。表2-7根据表2-7可写出以下因子得分函数：F1=0.238固定资产产值率+0.191固定资产利税率+0.265资金利润率+0.270资金利税率+0.173流动资金周转天数+0.135销售收入利税率-0.046全员劳动生

14、产率F2=-0. 087固定资产产值率+0.096固定资产利税率-0.126资金利润率-0.160资金利税率-0.627流动资金周转天数+0.167销售收入利税率+0.477全员劳动生产率单击“转换计算变量”，新建输出变量，命名为“score分”并输入其计算公式FAC1_1 * 0.66259 + FAC2_1 * 0.18360。将FAC1_1、FAC2_1和score原始变量标准化后，保存为新变量ZFAC1_1、ZFAC2_1和zscore然后，分别按ZFAC1_1、ZFAC2_1和zscore排序。得到因子得分值及排序表如表2-8所示。企业 盈利能力资金和人力综合得分ZFAC1_1名次Z

15、FAC2_1Zscore数源科技2.130561-0.64627111.88062中华电子2.211071.641632二轻1.23755-0.8872130.95573南方制药0.6955940.481850.799五羊自行0.2747650.07570.285石化-0.37374101.896590.146286白云制药0.19232-0.7109812-0.004527广发卷烟-0.381660.057-0.352588康佳电子-0.561410.69404-0.355699华南冰箱-0.09552-1.0999714-0.38577岭信-0.34094-0.24467-0.3939华空空调-0.58518-0.33462-0.65328中国长城-0.36788-1.2678915-0.69309稀土高科-0.87028-0.29465-0.91736三星集团-2.044980.07001-1.95203表2-8