多元统计分析课程毕业论文

1、题目：川理工学多元统计分析课程设计报告中国国有控股工业行业的经济效益评价学 生：雷鹏程何君李西京曾学成白俊明专 业：统计学指导教师：柏宏斌四川理工学院理学院二零一四年十二月中国国有控股工业行业的经济效益评价摘要本文主要研究了中国国有控股工业行业的经济效益，对反映行业经济效益的总资产贡 献率、资产负债率、流动资产周转次数、工业成本费用利润率和产品销售率等五个经济指 标进行主成分分析，提取反映行业盈利能力和市场能力的两个综合指标。然后通过因子分 析法分析反映经济效益的各指标的内部结构，表明行业经济效益主要由盈利能力和市场能 力两个公因子决定。根据各行业在盈利能力上的得分和市场能力上的得分将工业行业

2、分为 五类，并对各行业经济效益进行综合评价。然后用聚类分析对综合评价结果进行验证，表 明综合评价较为客观合理。最后，本文给出相应的政策建议。关键字：主成分分析、因子分析、聚类分析。一、引言改革开放以来，工业始终是我国经济发展的主要支柱。作为社会主义国家，我国国有 及国有控股工业行业掌控着国家工业发展命脉，对国民经济、社会协调发展具有巨大推动 作用。因此，考核工业行业的经济效益，对挖掘重点行业和弱势行业，提高整个国有工业 企业的经济效益等具有重大的现实意义。企业或行业的经济效益由众多因素来刻画，目前 反映行业经济效益主要有总资产贡献率、资产负债率、流动资产周转次数、工业成本费用 利润率和产品销售

3、率等五个经济指标10这些众多指标虽然能从多方面对行业的经济效益进 行全面考察，但也在一定程度增加了分析问题的复杂性。在损失少量信息的前提下，设计 一个或少数几个综合指标，并用较少的综合指标对工业经济效益进行分析评价，能够简化 问题。此外，挖掘出反映经济效益的众多指标的内在基本结构，有助于指出各行业经济效 益的主要决定因素及瓶颈，也有助于对各行业经济效益进行综合评价。二、文献综述大量国内文献从灰色系统理论、多元统计分析方法、层次分析法、模糊综合评判法、数据包络分析法等理论与方法，考察了中国各行业、企业或地区经济效益的研究与综合评 价。华中生、梁梁等用模糊聚类方法与数据包络分析分类法考察了合肥工业

4、行业的经济状 况，将各工业行业按经济效益的状况分为高、较高、一般、较差和差等五类1（华中生、梁梁，1995）。王树岭等人利用TOPSIS模型，对吉林省轻工业17个主要行业的经济效益 进行了综合评价与排序，确定出相应的优势行业（王树岭等，1999） o本文以2008年国有及国有控股的主要工业行业为研究对象，通过主成分分析和因子分析法，再次对各工业行 业的经济效益进行分析与评价，并结合聚类分析法来验证综合评价的结果。三、数据来源反映经济效益的指标较多，不同文献中选取的指标不尽相同。本文采用国家统计局最 新公布的五个指标：总资产贡献率、资产负债率、流动资产周转次数、工业成本费用利润 率和产品销售率，

5、分别记为 X1至X5 0总资产贡献率（XJ反映企业全部资产的获利能力。 资产负债率（X2 ）既反映企业经营风险的大小，也反映企业利用债权人提供的资金从事 经营活动的能力。流动资产周转次数（X3 ）反映投入工业企业流动资金的周转速度。成本 费用利润率（X4 ）反映企业投入的生产成本及费用的经济效益。产品销售率（X5 ）反映工业产品已实现销售的程度。选取 39个主要工业行业的数据整理如附录表 1所示。四、模型基本理论建立4.1主成分分析的基本理论设对某一事物的研究涉及p个指标，分别用X1，X2，Xp表示，这p个指标构成 的p维随机向量为X =（XX2,Xp）：设随机向量X的均值为"，协方

6、差矩阵为V 0对X进行线性变换，可以形成新的综合变量，用 丫表示，也就是说，新的综合向量可 以由原来的变量线性表示，即满足下式：国家统计年鉴2009年用这五大指标来反映工业行业的经济效益。由于可以任意地对原始变量进行上述线性变换，由不同的线性变换得到综合变量丫的统计特征也不尽相同。因此为了取得较好的效果，我们总希望 YUi X的方差尽可能大且 各Yi之间相互独立，由于 var(Y) = var(5 X) = 5瓦5 ，面对任意常数c，有FFvar(cui X) = c u ' 比因此对Ui不加限制时，可以使var(YJ任意增大，问题将变得没有意义。我们将线性变 换约束在下面的原则之下：

7、(1) ui ui 1 ( i = 1,2,，P)0(2) Yi 与 Yj 相互无关。(i = j ; i, j =1,2/ ,p 0(30 Yi是Xi，X2，Xp的所有线性组合中方差最大者；丫2是与Yi不相关的Xi， X2，,Xp的所有线性组合中方差最大者；， Yp是与Yi，Y2，Yp不相关的Xi， X2，,Xp的所有线性组合中方差最大者。基于以上这三条原则决定综合变量 Yi，丫2，Yp分别称为原始变量的第一，第二 第p个主成分。其中，各综合变量在总方差所占比重依次递减。在实际研究工作中，通常 指挑选前几个方差较大的主成分，从而达到简化系统结构，抓住问题实质的目的。4.2因子分析的基本理论设

8、有n个样品，每个样品观测p个指标，这些p个指标之间有较强的相关性。为了方 便研究，并消除由观测量纲的差异及数量级不同所造成的影响，将样本的观测数据进行标 准化处理，使标准化的变量均值为0,方差为i。为方便，把原始变量及标准后的变量向量均用X表示，用Fi, F2,Fm (m<p表示标准化的的公因子。如果：(10 X =(Xi,X ,Xp) 是可观测随机变量，且均值向量E(X) = 0 ，协方差矩阵cov (X0 = '： 且协方差矩阵V与相关阵R相等；(20 F -(Fi，F2/ , Fm ) (m<p)是不可观测变量,其均值向量E(F)=0，协方差矩阵(30cov(F)二I

9、 ,即向量F的各分量是相互独立的;=(；i,，；p)与 F相互独立，且E( 0 = 0，；的协方差矩阵a 是对角方阵:cov(；)八；=211J1r2即；的各分量之间也是相互独立的，则模型：称为因子模型，模型的矩阵形式为:其中A称为因子载荷矩阵，而在因子模型中，公共因子的个数少于原始变量的个数，且公 共因子是不可观测的隐变量，载荷矩阵 A不可逆，因而不能直接求得公共因子用原始变量 表示的精确线性组合。解决该问题的一种方法是用建立回归思想求出线性组合系数的估计 值，即建立如下以公共因子为因变量、原始变量为自变量的回归方程：Fj 二Xi2X2Xp j=1,2,m此处因为原始变量与公共因子变量均为标

10、准化量，因此回归方程中没有常数项。在最小二乘意义下，可以得到F估计值：式中，A为因子载荷矩阵；R原始变量的相关阵；X为原始变量向量。这样，在得到一组 样本值之后，就可以带入上面的关系式求出公共因子的估计得分，从而用少数的公共因子 去描述原始变量的数据结构，用公共因子得分去描述原始变量的取值，在估计出公共因子 得分后，可以利用因子得分去进行进一步分析，如样本点的聚类分析，当因子数m较少时,还可以方便地把各样本点在图上表示出来，直观地描述样本分布情况，从而便于把研究工 作引向深入。五、模型的求解与检验5.1工业行业经济效益的综合指标确定由主成分模型的基本原理可得，利用SPSS(19.0)软件对中国

11、国有工业行业的经济效益指标 进行主成分分析得到下表1、2:Total Variance ExplainedComponentInitial EigenvaluesExtraction Sums ofSquared LoadingsRotation Sums of SquaredLoadingsTotal% ofVarianceCumulative %Total% ofVarianceCumulative %Total% ofVarianceCumulative %123452.6831.184.630.376.12653.66923.68512.5937.5272.52553.66977.35

12、489.94897.475100.0002.6831.184.63053.66923.68512.59353.66977.35489.9482.4431.0391.01548.86820.77320.30748.86869.64189.948Extraction Method: Principal Component Analysis.表1解释总方差表Rotated Component Matrix123总资产贡献率X1( %.946.037-.071资产负债率X2( %-.847.085.159流动资产周转次数X3(次/年).050.976-.172成本费用利润率X4( %.892.211-

13、.154产品销售率X5( %-.180-.183.965表2旋转过后的主成分矩阵由上表可得，我们选取的5个经济效益指标被提取出了 3个主成分，提取的3个主成分集中了 5个原始量信息的89.948%。能够很好的反映5个经济效益指标。故三个公因子表示为：其中乙为原始变量的相应标准化变量。第一主成分对原始变量的贡献率为48.868%,第二主成分的贡献率为20.773%，第三个主成分对原始变量的贡献率为 20.307%,其累计 贡献率为89.948%,结果表明前三个个主成分提取了原始变量的绝大部分信息。由上可知,第一主成分中Zi、Z2和乙的系数绝对值较大，第二主成分中 Z3的系数绝对值较大,第三 个成

14、分中Z5的系数绝对值较大。因此，第一主成分主要由总资产贡献率、资产负债率和工 业成本费用利润率组成，第二主成分主要由流动资产周转次数和产品销售率组成。因此， 第一个公因子主要由“总资产贡献率”、“资产负债率”、“成本费用利用率”组成，该三个 指标主要反映一个行业的盈利能力，所以第一个综合指标可以表示为“盈利能力”，第二个公因子可得主要反映行业的运行能力，第二个综合指标可以表示为“行业运行能力”，第三个综合指标可以表示为“市场能力”。我们把5个影响行业的经济效益的自变量最终 用3个综合指标进行表示，接下来我们利用这三个综合指标对每个行业进行综合评价。5.2基于因子分析的经济效益的综合评价利用SP

15、S( 19.0)软件对中国国有工业行业的经济效益指标进行因子分析得到每个行业的 经济效益的最终评分，得到下表 KM创球形Bartlett检验结果：KMO and Bartlett's TestKaiser-Meyer-Olkin Measure of Sampling Adequacy.684Bartlett's Test ofApprox. Chi-Square83.527Sphericitydf10Sig.000表 3KMO口 Bartlett 检验首先巴特利特球度检验表明：巴特利特球度统计量值为83.527，显着性水平少于0.001, 即拒绝相关系数矩阵是一个单位阵的原假

16、设。KMO统计量为0.684，表明简单相关系数平 方和远大于偏相关系数平方和，比较适合做因子分析2。KMO值越接近1,则越适合做因子分析，反之亦反。Kaiser认为0.9以上就非常适合，0.8-09很适合，0.7-0.8适合,为了能够计算出各行业的因子得分，我们需要将公共因子表示成原始变量的线性组 合。公共因子对原始变量基于最小二乘法的最优线性估计为：F = atraz（ 2）其中A为利用“最大方差法”旋转后的因子载荷矩阵，R为样本相关系数矩阵，Z为标准化原始指标向量。按照（2）式，估计出三个个公因子的因子得分表达式为：Fi =0.946Zi 0.847Z2 0.05Z3 0.892Z 0.I

17、8Z5讦2 =0.037乙 +0.085Z2 +0.976Z3+0.211Z4-0.183Z5 （ 3）f3 = -0.070 +0.159Z2 -0.172Z3 -0.154Z4 + 0.965Z5各行业经济效益的综合指标得分：Soccer=48.868/89.948* F1 +20.713/89.948* F2 +20.307/89.948* F3（4）按照计算的综合指标得分公式，算出了如下的最终综合评价得分表：工业行业市场能力F1行业运行能力F2市场能力F3总得分（Socer）石油和天然气开采业烟草制品业其他采矿业黑色金属矿采选业饮料制造业有色金属矿采选业家具制造业石油加工、炼焦及核燃料

18、加工业电力、热力的生产和供应业煤炭开采和洗选业皮革、毛皮、羽毛（绒）及其制品业印刷业和记录媒介的复制金属制品业黑色金属冶炼及压延加工业有色金属冶炼及压延加工业非金属矿采选业交通运输设备制造业医药制造业燃气生产和供应业-0.006052878化学原料及化学制品制造业文教体育用品制造业-0.04804615仪器仪表、文化、办公机械制造业农副食品加工业食品制造业塑料制品业0.002260768非金属矿物制品业电气机械及器材制造业纺织服装、鞋、帽制造业橡胶制品业废弃资源和废旧材料回收加工业化学纤维制造业通信设备、计算机及其他电子设备制造 业工艺品及其他制造业通用设备制造业纺织业造纸及纸制品业专用设备制

19、造业水的生产和供应业木材加工及木、竹、草等制品业表4工业行业的因子得分及综合得分根据上述公式（2）、（ 3）和（4）,计算出工业主要行业的因子得分及排名（由高到低 进行排列）见上表4所示，其中可知石油和天然气开采业、烟草制品业和黑色金属矿采选 业等行业的经济效益最高；而其他采矿业、橡胶制品业和纺织业等行业的经济效益最低。 根据上表4，我们可以对39个工业行业的盈利能力、行业的运行能力及市场能力进行分类 5.3行业的聚类分析我们已经对各工业行业进行了综合评价，为了验证评价的客观性和可靠性，本文采用 聚类分析加以验证，其主要思想和依据：如果经济效益高的行业自然会成为一类，经济效 益差的行业会成为一

20、类，如果聚类分析能够验证这个结论，则说明我们的综合评价具有一 定的可靠性。本文用欧式距离来度量样本之间的距离，用瓦尔德法测度来类与类之间的距 离，将39个行业分为四类：第一类：石油和天然气开采业、烟草制造业。第二类：黑色金属矿采选业、有色金属矿采选业、非金属矿采选业、饮料制造业、皮 革毛皮羽毛（绒）及其制品业。第三类：印刷业和记录媒介的复制、文教体育用品制造业、医药制造业、仪器仪表及 文化与办公用机械制造业、水的生产和供应业。第四类为其它行业。结合表4中行业的综合经济效益可知：聚类分析中第一类为经济效益最好的两个行业， 第二类中除水的生产和供应业和文教体育用品制造业的综合排名较后，分别为26和

21、30名,而它们在盈利能力的排名分别是13名和14名。仪器仪表及文化与办公用机械制造业排在第 13名，其它行业的经济效益均排在前十名。 即聚类分析中的第二类为经济效益较好的行业， 第三类为经济效益较差的行业。 因此，聚类分析基本验证了基于因子分析的综合评价模型， 说明表4中所示的综合评价比较合理。图 1 聚类树状图六 问题的思考及政策建议6.1 问题思考一个企业或行业的经济效益由众多因素来刻画，目前反映行业经济效益主要有总资产 贡献率、资产负债率、流动资产周转次数、工业成本费用利润率和产品销售率等五个经济 指标。这些众多指标虽然能从多方面对行业的经济效益进行全面考察，但也在一定程度增 加了分析问

22、题的复杂性。在损失少量信息的前提下，设计一个或少数几个综合指标，并用 较少的综合指标对工业经济效益进行分析评价，能够简化问题。此外，挖掘出反映经济效 益的众多指标的内在基本结构，有助于指出各行业经济效益的主要决定因素及瓶颈，也有 助于对各行业经济效益进行综合评价。所以如何对一个企业的经济效益的评价对当下的企 业和全国的经济增长有着直接的作用，如何去发展一个企业的经济也是至关重要的，对当 下企业的经济效益有一定的益处。6.2 政策建议 本文以2008年中国国有及国有控股行业进行较为深入的分析与研究，据此提出以下几 点建议：一. 在我国国有及国有控股工业行业中，各行业的经济效益及其潜在的优势和劣势

23、不 尽相同。政府部门在扶持行业的发展时，应该同时考虑现阶段行业的优势和劣势，做到强 者更强，弱者不弱。二尽管分析表明，行业经济效益的主要因子是行业的盈利能力（占比 53.65%），但 政府部门在决定支持行业的发展时不能只着重行业的盈利能力。比如烟草制造品业的市场 盈利能力为全行业最高，但其市场能力为全行业倒数第二。要想尽快提高工业行业的经济 效益，应该首先重点扶持盈利能力和市场能力均较高的行业，如石油和天然气开采业、黑 色金属矿采选业、非金属矿采选业等。三. 合理引导各行业的投资，不能只将投资挤向盈利能力较高的行业，而要适度考虑 到市场能力较强的行业。七 参考文献1 华中生，梁梁 . 地区工业

24、行业经济状况的综合评价与分析 J. 管理工程学报 1995，6.2 王树岭等 . 吉林省轻工业行业经济效益综合评价 J. 吉林工业大学学报自然科学版，1999, 1.3 何晓群.多元统计分析M.北京：中国人民大学出版社，20044 林秀梅.多元统计方法M.长春：吉林人民出版社，19962013.3 张文彤，董伟.SPSS统计分析高级教程（第二版）M.北京：高等教育出版社,附录表1： 2008年国有及国有控股工业行业主要经济效益指标的原始数据行业序号总资产贡献率人（%资产负债率X（ %流动资产周转次数人（次/年）成本费用利润率X4（ %产品销售率 X5（ %煤炭开采和洗选业116.2961.04

25、1.9217.3198.2石油和天然气开采业249.4739.064.6583.4789.16黑色金属矿采选业319.7140.691.932.1996.26有色金属矿采选业417.8547.532.0317.298.88非金属矿采选业513.1252.511.7911.5496.26其他采矿业60.6374.370.370.51118.57农副食品加工业78.9565.542.82.6696.5食品制造业88.7863.22.253.0897.33饮料制造业921.4342.611.6117.6799.07烟草制品业107023.351.4941.9498.94纺织业114.2767.871

26、.610.298.1纺织服装、 鞋、帽制造业128.2964.651.375.6798.82皮革、毛皮、 羽毛（绒）及 其制品业1316.6743.221.2211.9298.27木材加工及 木、竹、草等 制品业148.2564.041.773.3293.61家具制造业1515.72582.559.2299.67造纸及纸制品业166.5359.821.65.0694.99印刷业和记 录媒介的复制1712.941.571.6611.9897.56文教体育用185.3750.581.382.7199.94品制造业石油加工、炼焦及核燃料 加工业19-6.9758.885.54-7.6499.15化学原料及化学制品制造业206.255