多元统计分析课程论文

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

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

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

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

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

6、，用表示，也就是说，新的综合向量可以由原来的变量线性表示，即满足下式：由于可以任意地对原始变量进行上述线性变换，由不同的线性变换得到综合变量的统计特征也不尽相同。因此为了取得较好的效果，我们总希望的方差尽可能大且各之间相互独立，由于 ，面对任意常数，有 因此对不加限制时，可以使任意增大，问题将变得没有意义。我们将线性变换约束在下面的原则之下： （1）（）。 （2）与相互无关。（；）（3）是，, 的所有线性组合中方差最大者；是与不相关的，, 的所有线性组合中方差最大者；，是与，不相关的，, 的所有线性组合中方差最大者。基于以上这三条原则决定综合变量，分别称为原始变量的第一，第二第p个主成分。其中

7、，各综合变量在总方差所占比重依次递减。在实际研究工作中，通常指挑选前几个方差较大的主成分，从而达到简化系统结构，抓住问题实质的目的。4.2因子分析的基本理论设有n个样品，每个样品观测p个指标，这些p个指标之间有较强的相关性。为了方便研究，并消除由观测量纲的差异及数量级不同所造成的影响，将样本的观测数据进行标准化处理，使标准化的变量均值为0，方差为1。为方便，把原始变量及标准后的变量向量均用表示，用, （m<p）表示标准化的的公因子。如果：（1） 是可观测随机变量，且均值向量 ，协方差矩阵 = ，且协方差矩阵与相关阵相等；（2） (m<p)是不可观测变量,其均值向量，协方差矩阵,即向

8、量的各分量是相互独立的；（3） 与相互独立，且，的协方差矩阵是对角方阵：=即的各分量之间也是相互独立的，则模型：称为因子模型，模型的矩阵形式为：其中称为因子载荷矩阵，而在因子模型中，公共因子的个数少于原始变量的个数，且公共因子是不可观测的隐变量，载荷矩阵不可逆，因而不能直接求得公共因子用原始变量表示的精确线性组合。解决该问题的一种方法是用建立回归思想求出线性组合系数的估计值，即建立如下以公共因子为因变量、原始变量为自变量的回归方程： j=1,2,m此处因为原始变量与公共因子变量均为标准化量，因此回归方程中没有常数项。在最小二乘意义下，可以得到估计值：式中，为因子载荷矩阵；原始变量的相关阵；为原

9、始变量向量。这样，在得到一组样本值之后，就可以带入上面的关系式求出公共因子的估计得分，从而用少数的公共因子去描述原始变量的数据结构，用公共因子得分去描述原始变量的取值，在估计出公共因子得分后，可以利用因子得分去进行进一步分析，如样本点的聚类分析，当因子数m较少时，还可以方便地把各样本点在图上表示出来，直观地描述样本分布情况，从而便于把研究工作引向深入。五、模型的求解与检验5.1工业行业经济效益的综合指标确定由主成分模型的基本原理可得，利用spss(19.0)软件对中国国有工业行业的经济效益指标进行主成分分析得到下表1、2：total variance explainedcomponentini

10、tial eigenvaluesextraction sums of squared loadingsrotation sums of squared loadingstotal% of variancecumulative %total% of variancecumulative %total% of variancecumulative %12.68353.66953.6692.68353.66953.6692.44348.86848.86821.18423.68577.3541.18423.68577.3541.03920.77369.6413.63012.59389.948.6301

11、2.59389.9481.01520.30789.9484.3767.52797.4755.1262.525100.000extraction method: principal component analysis.表1解释总方差表rotated component matrixacomponent123总资产贡献率x1（%）.946.037-.071资产负债率x2（%）-.847.085.159流动资产周转次数x3（次/年）.050.976-.172成本费用利润率x4（%）.892.211-.154产品销售率x5（%）-.180-.183.965表2旋转过后的主成分矩阵由上表可得，我们选取

12、的5个经济效益指标被提取出了3个主成分，提取的3个主成分集中了5个原始量信息的89.948%。能够很好的反映5个经济效益指标。故三个公因子表示为：其中为原始变量的相应标准化变量。第一主成分对原始变量的贡献率为48.868%，第二主成分的贡献率为20.773% ，第三个主成分对原始变量的贡献率为20.307%，其累计贡献率为89.948%，结果表明前三个个主成分提取了原始变量的绝大部分信息。由上可知，第一主成分中、和的系数绝对值较大，第二主成分中的系数绝对值较大,第三个成分中的系数绝对值较大。因此，第一主成分主要由总资产贡献率、资产负债率和工业成本费用利润率组成，第二主成分主要由流动资产周转次数

13、和产品销售率组成。因此，第一个公因子主要由“总资产贡献率”、“资产负债率”、“成本费用利用率”组成，该三个指标主要反映一个行业的盈利能力，所以第一个综合指标可以表示为“盈利能力”，第二个公因子可得主要反映行业的运行能力，第二个综合指标可以表示为“行业运行能力”，第三个综合指标可以表示为“市场能力”。我们把5个影响行业的经济效益的自变量最终用3个综合指标进行表示，接下来我们利用这三个综合指标对每个行业进行综合评价。5.2基于因子分析的经济效益的综合评价利用spss（19.0）软件对中国国有工业行业的经济效益指标进行因子分析得到每个行业的经济效益的最终评分，得到下表kmo和球形bartlett检验

14、结果：kmo and bartlett's testkaiser-meyer-olkin measure of sampling adequacy.684bartlett's test of sphericityapprox. chi-square83.527df10sig.000表3kmo和bartlett检验首先巴特利特球度检验表明：巴特利特球度统计量值为83.527，显著性水平少于0.001,即拒绝相关系数矩阵是一个单位阵的原假设。kmo统计量为0.684，表明简单相关系数平方和远大于偏相关系数平方和，比较适合做因子分析 kmo值越接近1，则越适合做因子分析，反之亦反。k

15、aiser认为0.9以上就非常适合，0.8-09很适合，0.7-0.8适合，0.6-0.7比较适合，0.5-0.6勉强，0.5以下不适合。为了能够计算出各行业的因子得分，我们需要将公共因子表示成原始变量的线性组合。公共因子对原始变量基于最小二乘法的最优线性估计为： （2）其中a为利用“最大方差法”旋转后的因子载荷矩阵，r为样本相关系数矩阵，z为标准化原始指标向量。按照（2）式，估计出三个个公因子的因子得分表达式为：（3）各行业经济效益的综合指标得分：soccer=48.868/89.948*+20.713/89.948*+20.307/89.948* （4）按照计算的综合指标得分公式，算出了如

16、下的最终综合评价得分表：工业行业市场能力行业运行能力市场能力总得分（socer）石油和天然气开采业 3.2366072352.847904035-1.1085059192.165869204烟草制品业 3.880133782-0.7179905560.9018566262.145834953其他采矿业 -0.338689942-0.5795507495.1968800280.85541574黑色金属矿采选业 1.316706228-0.193679604-0.3145251920.599617895饮料制造业 1.065204838-0.4906709430.3344165740.540898

17、295有色金属矿采选业 0.751264250.0178027970.3327810470.48739698家具制造业 0.1487616410.6869083020.604894760.376022012石油加工、炼焦及核燃料加工业 -1.1727125093.6538721240.649468680.35334339电力、热力的生产和供应业 -0.6139139662.0479883340.5792198380.270204243煤炭开采和洗选业 0.238700680.0928350570.1574905120.186679496皮革、毛皮、羽毛(绒)及其制品业 0.736082724-

18、0.979664216-0.0407083740.164469036印刷业和记录媒介的复制 0.625904514-0.573691999-0.2019406840.161966785金属制品业 -0.1580977580.1746264710.2555480420.012129325黑色金属冶炼及压延加工业 -0.4708672750.6927576830.34212949-0.018589219有色金属冶炼及压延加工业 -0.2718005140.3981640940.020683685-0.051043505非金属矿采选业 0.210968339-0.322698036-0.495336

19、648-0.071737079交通运输设备制造业 -0.201924849-0.2016348490.352617116-0.076662388医药制造业 0.442582731-0.830304619-0.633719878-0.094373799燃气生产和供应业 -0.2470875870.066344551-0.006052878-0.120285239化学原料及化学制品制造业 -0.3088391930.261045862-0.073210663-0.124030961文教体育用品制造业 -0.04804615-0.7058428030.285857441-0.124577364仪器仪

20、表、文化、办公机械制造业 0.225149559-1.011153889-0.207111795-0.157956934农副食品加工业 -0.6036182430.868965399-0.304020917-0.195894859食品制造业 -0.4703772590.295955823-0.184548633-0.228867064塑料制品业 -0.345529774-0.2575129140.002260768-0.246684255非金属矿物制品业 -0.297512988-0.223117102-0.209806691-0.260530759电气机械及器材制造业 -0.32157923

21、4-0.487846883-0.086675011-0.306944977纺织服装、鞋、帽制造业 -0.382512265-0.5103711260.078181656-0.308032573橡胶制品业 -0.894252540.3144376430.188327868-0.370704695废弃资源和废旧材料回收加工业 -0.4654374050.796076505-1.419027766-0.389383808化学纤维制造业 -0.664949631-0.253473642-0.213264699-0.467947401通信设备、计算机及其他电子设备制造业-0.41880571-0.713

22、185265-0.428322213-0.488939544工艺品及其他制造业 -0.781779398-0.180031692-0.189109901-0.509005745通用设备制造业 -0.587485138-0.663378194-0.224067297-0.522965642纺织业 -0.788943911-0.29145315-0.142298644-0.528062057造纸及纸制品业 -0.440990493-0.520960211-0.964197908-0.577580344专用设备制造业 -0.633569735-0.680842747-0.766897175-0.67

23、4587685水的生产和供应业 -0.32667013-1.472493298-0.782511578-0.694204227木材加工及木、竹、草等制品业 -0.622072925-0.354136195-1.286753666-0.710255232表4工业行业的因子得分及综合得分根据上述公式（2）、（3）和（4），计算出工业主要行业的因子得分及排名（由高到低进行排列）见上表4所示，其中可知石油和天然气开采业、烟草制品业和黑色金属矿采选业等行业的经济效益最高；而其他采矿业、橡胶制品业和纺织业等行业的经济效益最低。根据上表4，我们可以对39个工业行业的盈利能力、行业的运行能力及市场能力进行分类

24、5.3 行业的聚类分析我们已经对各工业行业进行了综合评价，为了验证评价的客观性和可靠性，本文采用聚类分析加以验证，其主要思想和依据：如果经济效益高的行业自然会成为一类，经济效益差的行业会成为一类，如果聚类分析能够验证这个结论，则说明我们的综合评价具有一定的可靠性。本文用欧式距离来度量样本之间的距离，用瓦尔德法测度来类与类之间的距离，将39个行业分为四类：第一类：石油和天然气开采业、烟草制造业。第二类：黑色金属矿采选业、有色金属矿采选业、非金属矿采选业、饮料制造业、皮革毛皮羽毛(绒)及其制品业。第三类：印刷业和记录媒介的复制、文教体育用品制造业、医药制造业、仪器仪表及文化与办公用机械制造业、水的

25、生产和供应业。第四类为其它行业。结合表4中行业的综合经济效益可知：聚类分析中第一类为经济效益最好的两个行业，第二类中除水的生产和供应业和文教体育用品制造业的综合排名较后，分别为26和30名，而它们在盈利能力的排名分别是13名和14名。仪器仪表及文化与办公用机械制造业排在第13名，其它行业的经济效益均排在前十名。即聚类分析中的第二类为经济效益较好的行业，第三类为经济效益较差的行业。因此，聚类分析基本验证了基于因子分析的综合评价模型，说明表4中所示的综合评价比较合理。图1聚类树状图六 问题的思考及政策建议6.1问题思考一个企业或行业的经济效益由众多因素来刻画，目前反映行业经济效益主要有总资产贡献率

26、、资产负债率、流动资产周转次数、工业成本费用利润率和产品销售率等五个经济指标。这些众多指标虽然能从多方面对行业的经济效益进行全面考察，但也在一定程度增加了分析问题的复杂性。在损失少量信息的前提下，设计一个或少数几个综合指标，并用较少的综合指标对工业经济效益进行分析评价，能够简化问题。此外，挖掘出反映经济效益的众多指标的内在基本结构，有助于指出各行业经济效益的主要决定因素及瓶颈，也有助于对各行业经济效益进行综合评价。所以如何对一个企业的经济效益的评价对当下的企业和全国的经济增长有着直接的作用，如何去发展一个企业的经济也是至关重要的，对当下企业的经济效益有一定的益处。6.2政策建议本文以2008年

27、中国国有及国有控股行业进行较为深入的分析与研究，据此提出以下几点建议：一.在我国国有及国有控股工业行业中，各行业的经济效益及其潜在的优势和劣势不尽相同。政府部门在扶持行业的发展时，应该同时考虑现阶段行业的优势和劣势，做到强者更强，弱者不弱。二尽管分析表明，行业经济效益的主要因子是行业的盈利能力（占比53.65%），但政府部门在决定支持行业的发展时不能只着重行业的盈利能力。比如烟草制造品业的市场盈利能力为全行业最高，但其市场能力为全行业倒数第二。要想尽快提高工业行业的经济效益，应该首先重点扶持盈利能力和市场能力均较高的行业，如石油和天然气开采业、黑色金属矿采选业、非金属矿采选业等。三.合理引导各

28、行业的投资，不能只将投资挤向盈利能力较高的行业，而要适度考虑到市场能力较强的行业。七 参考文献1 华中生，梁梁. 地区工业行业经济状况的综合评价与分析j.管理工程学报1995，6.2 王树岭等.吉林省轻工业行业经济效益综合评价j.吉林工业大学学报自然科学版，1999，1.3 何晓群.多元统计分析m.北京：中国人民大学出版社，20044 林秀梅.多元统计方法m.长春：吉林人民出版社，19965张文彤，董伟.spss统计分析高级教程（第二版）m.北京：高等教育出版社，2013.3附录表1：2008年国有及国有控股工业行业主要经济效益指标的原始数据行业序号总资产贡献率x1（%）资产负债率x2（%）流

29、动资产周转次数x3（次/年）成本费用利润率x4（%）产品销售率x5（%）煤炭开采和洗选业116.2961.041.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.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.9