多元统计分析方法概述.docx

多元统计分析方法概述目录引言第四页多元线性回归方法原理简介第四页多元线性回归案例叙述分析第四页多元线性回归分析方法在社会的应用第八页聚类分析方法原理简介第八页聚类分析案例叙述分析第八页聚类分析方法在社会的应用第十页主成份分析方法原理简介第十页主成份分析案例叙述分析第十一页主成份分析方法在社会的应用第十四页因子分析方法原理简述第十四页因子分析案例叙述分析第十四页因子分析方法在社会的应用第十七页偏最小二乘回归分析方法原理简介第十八页偏最小二乘回归分析案例叙述分析第十九页偏最小二乘回归分析方法在社会的应用第二十一页总结第二十一页参考文献第二十二页谢辞第二十三页摘要本文主要概述了多元统计分析的各个方法，然后在后面介绍了多元统计分析方法在社会生活等方面的实际案例以及分析。并由案例分析找出各个统计分析方法的对应使用领域。关键词多元统计分方法回归分析聚类分析因子分析主成份分析偏最小二乘回归分析因素股市模型财务SummaryTheChemometricsincludeschemicalexperimentaldesignandoptimization(suchasorthogonaldesign,simplexmethodandvarianceanalysis),chemicalpatternrecognition(suchasclusters,PCA,k-nearestneighbouranalysis,SIMCAandANN),multi-variancecalibration(suchasM1.R,C1.S,PCRandP1.S)andspectrumanalysis(SuchasITTFA,EFAandFSWEFA),signalprocessing(suchasfiltering,smoothing,derivationandconvolution).keywordmultivariatestatisticalanalysismethodregressionanalysisclusteranalysisfactoranalysisprincipalcomponentanalysislinearleastsquaresestimatecomplicationequitymarketmodelfinance一、引言多元统计分析的基本方法。多元统计分析是数理统计学30多年来迅速发展起来的一个分支随着计算机的普及，各种统计软件不断推出，多元统计分析方法已广泛应用到教育管理的各个领域。多元统计分析是运用数理统计的方法来研究多变量问题的理论和方法，它是一元统计学的推广在许多教育问题中，教育现象涉及到的变量不是一个，而是多个变量，且这些变量间又存在一定的联系，需要处理多个变量的观测数据，如果用一元统计方法就要对多方面分别进行分析，而一次分析一个方面，同时忽视了各方面之间存在的相关性，这样会丢失很多信息,分析的结果不能客观全面地反映情况。多元统计分析方法主要包括线性回归分析方法、判别分析方法、聚类分析方法、主成份分析方法、因子分析方法、对应分析方法、典型相关分析方法以及片最小二乘回归分析方法等。在这里我们主要针对了最常用的回归分析方法、聚类分析方法、主成份分析方法、因子分析方法以及偏最小二乘回归分析方法结合实际案例进行分析总结。最终总结出各个方法所使用的领域。二、多元线性回归分析方法（一）多元线性回归分析方法原理简介根据相关性的大小把变量分组，使得同组内的变量之间的相关性较高，但不同组的变量相关性较低各个因子间互不相关，所有变量都可以表示成公因子的线性组合。因子分析的目的就是减少变量的数目，用少数因子代替所有变量去分析整个问题。（二）多元线性回归案例叙述分析公路客货运输量多元线性回归预测方法探讨1.背景公路客、货运输量的定量预测，近几年来在我国公路运输领域大面积广泛地开展起来，并有效的促进了公路运输经营决策的科学化和现代化。关于公路客、货运输量的定量预测方法很多，本文主要介绍多元线性回归方法在公路客货运输量预测中的具体操作。根据笔者先后参加的部、省、市的科研课题的实践，证明了多元线性回归方法是对公路客、货运输量预测的一种置信度较高的有效方法。2 .多元线性回归预测线性回归分析法是以相关性原理为基础的.相关性原理是预测学中的基本原理之一。由于公路客、货运输量受社会经济有关因素的综合影响。所以，多元线性回归预测首先是建立公路客、货运输量与其有关影响因素之间线性关系的数学模型。然后通过对各影响因素未来值的预测推算出公路客货运输量的预测值。3 .公路客、货运输量多元线性回归预测方法的实施步骤：影响因素的确定影响公路客货运输量的因素很多，主要包括以下一些因素：（1）客运量影响因素人口增长量裤保有量、国民生产总值、国民收入工农业总产值，基本建设投资额城乡居民储蓄额铁路和水运客运量等。货运量影响因素人口货车保有量（包括拖拉机），国民生产总值，国民收入、工农业总产值，基本建设投资额，主要工农业产品产量，社会商品购买力，社会商品零售总额.铁路和水运货运量等。上述影响因素仅是对一般而言，在针对具体研究对象时会有所增减。因此，在建立模型时只须列入重要的影响因素，对于非重要因素可不列入模型中。若疏漏了某些重要的影响因素，则会造成预测结果的失真。另外，影响因素太少会造成模型的敏感性太强.反之，若将非重要影响因素列入模型，则会增加计算工作量，使模型的建立复杂化并增大随机误差。影响因素的选择是建立预测模型首要的关键环节，可采取定性和定量相结合的方法进行.影响因素的确定可以通过专家调查法，其目的是为了充分发挥专家的聪明才智和经验。具体做法就是通过对长期从事该地区公路运输企业和运输管理部门的领导干部、专家、工作人员和行家进行调查。可通过组织召开座谈会.也可以通过采访，填写调查表等方法进行，从中选出主要影响因素为了避免影响因素确定的随意性，提高回归模型的精度和减少预测工作量，可通过查阅有关统计资料后，再对各影响因素进行相关度（或关联度）和共线性分析，从而再次筛选出最主要的影响因素.所谓相关度分析就是将各影响因素的时间序列与公路客货运量的时间序列做相关分杯事先确定一个相关系数，对相关系数小于的影响因素进行淘汰.关联度是灰色系统理论中反映事物发展变化过程中各因素之间的关联程度，可通过建空公路客、货运量与各影响影响因素之间关联系数矩阵，按一定的标准系数舍去关联度小的影响因素.所谓共线性是指某些影响因素之问存在着线性关系或接近于线性关系.由于公路运输经济自身的特点，影响公路客，货运输量的诸多因素之问总是存在着一定的相关性，持别是与国民经济有关的一些价值型指标。4 .建立经验线性回归方程利用最小二乘法原理寻求使误差平方和达到撮小的经验线性回归方程：y=+QlCl+Q2N2+.+anxny一一预测的客、货运量g一一各主要影响因数5 .数据整理对收集的历年客、货运输量和各主要影响因素的统计资料进行审核和加工整理是为了保证预测工作的质量。资料整理主要包括下列内容：（1）资料的补缺和推算。（2）对不可靠资料加以核实调整.对查明原因的异常值加以修正。（3）对时间序列中不可比的资料加以调整和规范化；对按当年价格计算的价值指标应折算成按统。6 .多元线性回归模型的参数估计在经验线性回归模型中，源是要估计的参数，可通过数理统计理论建立模型来确定。在实际预测中，可利用多元线性回归复相关分析的计算机程序来实现对模型参数的估计值进行检验。此项工作的目的在于判定估计值是否满意、可靠。一般检验工作须从以下几方面来进行。（1）经济意义检验（2）统计检验(3)拟合度检验(4)回归方程的显著性检验(5)参数估计值的标准差检验应当强调指出.统计检验相对于经济意义检验来说是第二位的。如果经济意义检验不合理，那么即使统计检验可以达到很高的置信度，也应当抛弃这种估计结果，因为用这样的结果来进行经济预测是没有意义的。7 .最优回归方程的确定经过上述的经济意义和统计检验后，挑选出的线性回归方程往往是好几个、为了从中优选出用于进行实际预测的方程，我们可以采用定性和定量相结合的办法。从数理统计的原理来讲，应挑选方程的剩余均方和S-E较小为好.但作为经济预'删还必须尽量考虑到方程中的影响因素更切合实际和其未来值更易把握的原则来综合考虑。当然、有时也可以从中挑选出好几个较优的回归方程.通过预测后，分别作为不同的高、中、低方案以供决策人员选择。8 .模型的实际预测检验在获得模型参数估计值后，又经过了上述一系列检验而选出的最优(或较优)回归方程，还必须对模型的预测能力加以检验。不难理解、最优回归方程对于样本期间来说是正确的，但是对用于实际预测是否合适呢?为此，还必须研究参数估计值的稳定性及相对于样本容量变化时的灵敏度，也必须研究确定估计出来的模型是否可以用于样本观察值以外的范国，其具休做法是：(1)采用把增大样本容量以后模型估计的结果与原来的估计结果进行比较，并检验其差异的显著性。(2)把估计出来的模型用于样本以外某一时间的实际预测，并将这个预测值与实际的观察值作一比较，然后检验其差异的显著性。9 .模型的应用公路客、货运输量多元线性回归预测模型的研究目的主要有以下几个方面。(1)进行结构分析，研究影响该地区的公路客、货运输量的主耍因素和各影响因素影响程度的大小，进一步探讨该地区公路运输经济理论。(2)预测该地区今后年份的公路客、货运输量的变化，以便为公路运输市场、公路运输政策及公路运辅建设项目投资作出正确决策提供理论依据。另外，还可以通过公路客.货运输量与公路交通量作相关分析来对公路的饱和度发展趋势进行预测。从而为公路的新建、扩建项目的投资提供决策分析。(3)模拟各种经济政策下的经济效果，以便对有关政策进行评价。四、经调查分析，影响某地区旅客运输量的因素为。川国民收入X2工农业总产值%3社会总产值X人口Xs客车保有量及一一城乡居民储蓄存款经计算得下列相关系数表：X】Xi否X吊在Y0.94390.92870.90430.99140.96700.7021Z0.97360.96140.93260.86450.93210.6678Y客运盈Z旅客周转量若令=0.85,则可以舍去及这个影响因索，也就是认为“城乡居民储蓄存款”不能作为响旅客运输量的主要因素。2.经调查分析、影响某地区旅客运输量的因素为：X国民收入Xi工农业总产值%3一社会总产值Xx人口XB一一客车保有量X6国民生产总值Xi公路通车里程经计算得客运量和旅客周转量的经验线性回归方程如下：Y=0+QH+5X5R2=0.9997V=1+a；©+*3+Q4R2=0.9962Z=Bo+B4X4+B5吊+1X1R2=0.9983Z=+RN4+限5+房叫R=O.9990Y客运盈Z旅客周转量各自变量问的相关系数表如下:rU5×23%与x0.98360.96210.97100.93420.90360.9344叼0.98360.95230.96070.94660.98250.9823×30.9621O.8S230.89240.92050.97730.9114×40.97100.96070.89240.93480.97660.89280.93420.94660.92050.93480.98170.9036%090360.9825097730.97660.98170.9715X70.9344098230.91140.89280.90360.9715由上述计算可知，四个方程中均未出现rj>的情况.因此可以认为各自方程中的影响因素之间不存在严重共线性问题。3.经调查分析，影响某地区货运周转量的因素为：X国民收入%工农业总产值%3基建投资额X、原煤产量后一一钢铁、化肥、水泥、粮食总产量%6国民总产值用一一社会商品零售总额不一一相邻地、市工农业总产值的平均值V=4+国+%Xg+aix1（1）其中：4=0.9875F=206.33SE=1673.24f4=-2.8321备=3.1407t7=2.7431/=+b1x2+b、X（2）其中：#=0.9764F=164.59SE=1044.27（三）多元线性回归分析方法在社会的应用由上述案例分析可知多元线性回归分析方法在通过线性约束等条件将一个事物进行多元的分析处理，最终能筛选出影响这个事物发展的因素。这样就能通过多元线性回归分析将一个事物进行资源最优化配置，在交通、航运的等领域都比较实用。三、聚类分析方法（一）聚类分析方法原理简介聚类分析是研究事物分类的一种方法，是将一批样本或变量按照它们在性质上的亲疏程度加以分类。实质是按照距离的远近将数据分为若干个类别，以使得类别内数据的差异尽可能小，类别间的“差异”尽可能大。聚类分析中包括:层次聚类法，费层次聚类法，智能聚类法等多种详细的方法。（二）聚类分析案例叙述分析1 .问题的提出随着我国市场经济建设的高速发展，人们的金融意识和投资意识日益增强，而作为市场经济的组成部分一股票市场，正逐步走向成熟与规范，越来越多的投资者把眼光投向了股票，历史已经证明股票是一种不仅在过去已提供了投资者可观的长期利益，并且在将来也将提供良好机遇的投资媒体。然而，股价涨跌无常,股市变幻莫测，投资者要想在股市投资中赢取丰厚的投资回报，成为一个成功的投资者，就得认真研究上市公司的历史、业绩和发展前景，详细分析上市公司的财务状况，树立以基本分析为主，技术分析为辅的投资理念，找出真正具有投资价值的股票，进行长期投资。股票投资的基本分析分为宏观分析、中观分析、微观分析三大部分，宏观分析指对国家的国民经济以及政治、文化的分析，微观分析指公司分析，而中观分析主要指行业分析和地区分析等，板块分析主要归属于中观分析，兼有微观分析。中国股市从无到有，发展至今己颇具规模，前些年，在中国股市发展的初生阶段,由于市场规模小，上市公司数量不多，加上股民的投资思维和操作方法不太成熟,因此，投机性特强，这时用不上多少板块分析。但是，随着股市发展、投资手法和证券监管方法的成熟，以及上市公司数量的不断增多，如果再和以往一样，面对上千种股票胡乱抓一气，碰运气，甚至受各种股评和谣言所左右，则很难走向理性化，进而难以最终取得投资成功。因此，在成熟股市中，一个股民若想成功,必须学会板块分析，习惯理性操作，树立板块投资理念。2 .聚类分析在股市板块分析中的应用系统聚类分析的基本思想是首先将每个样本当作一类，然后根据样本之间的相似程度并类，并计算新类与其它类之间的距离，再选择相近者并类，每合并一次减少一类，继续这一过程，直到所有样本都并成一类为止。在聚类过程中，我们选用欧氏距离来度量类与类之间的相似程度，聚类方法采用类平均法。我们以高科技板块中的31个上市公司为研究对象，分析中选取了这31个上市公司1997年的每股收益、每股净资产、股东权益率、净资产收益率、净利润率等五个反映上市公司综合盈利能力的指标，数据取自4（略）。应用SAS软件中的系统聚类过程C1.USTER对31个样本进行聚类，得到表1所示的聚类过程。表1中NC1.为聚类数；ClustersJoined为每次聚成一个新类的2个样品（标有OB）或旧类（标有C1.）；FREQ为新类中所含有的样品数；SPRSQ为半偏R2,它表示每一次合并对信息的损失程度，看这一列的数值可知：从4类合并成3类时信息损失（为0.1042）最多，此统计量表明聚成4类较合适；CCC在NC1.=4时达到唯一的峰值-2.49,它支持分4类；PSF为伪F统计量，PSF出现峰值时所对应的分类数较合适，从这一列的数值可知PSF在NC1.=I7、NC1.=I2、NC1.=9和NC1.=4时4次达到峰值，但在NC1.=4时峰更陡些；PST2为t2统计量，PST2出现峰值的前一行所对应的分类数较合适，从这一列的数值可知NC1.=3时出现峰值9.8,它也支持分4类。综合这四个统计量可知：将31个样本分成4类较合适。根据表1,作出图1所示的聚类谱系图。当我们取分类阀值为1.0时，31个样本被分成了各类包含的样本如下：表1系统聚类过程NC1.Cluster三JoinedFREQSPRSQCCCPSFFST2DIST300B20B1920.1819.00.2330290B160B3020.2616.10.2792260B90B2320.002615.70.2804270B250B27202715.702839ZA0B60B2820.4114.30.351625C1.3O0B2030005612.63.10.3751240B3C1.2730.5811.72.20.388023OBl0B26206311.20.434022C1.25C1.250one9.53.40.449221OBlT0B3120.709.50.459920C1.240B540.00879.22.10.465419OBlO0B1520759404744180B40B820829.50495217OBll0B21200082970496516C1.26C1.2940.0133944.00499815C1.200B750Olll9.S1.90.522814CI3C1.140.01439.5200.5680130B14C1.Zl30.0129971.80.5857120B130B2420.0141990651211C1.22C1.l690.04278.57.20.6546100B12C1.1340.01858.81.90.6848i9C1.l90B2230.0224,9.13.007485IiC1.llC1.l5140.08687.09.20.794370B180B2920.02137.80.79896C1.l4C1.lT60.0492-3.12785.30.85015C1.8C1.12160.0573HS8.03.80.9007i4C1.9C1.lO70.0703-2.4984S.10.9533C1.6C1.5220.1745-3.55629.81.0341'2C1.3C1.4290.1656-3.494.8671.0618C1.2C1.T310.1428-0.004.81.3179第一类:0B18（天津磁卡），0B29（燕化高新）；第二类：0B22（深科技），0B15（实达电脑）,OBlO（清华同方）,0B12（东大阿派），0B14（长城电脑）,0B31（风华高科），0B17（东方通信）；第三类：OBl（华光科技），0B26（佛山照明），0B4（国脉通信），0B8（工大高新）,OB11（振华科技），0B21（彩虹股份）；第四类：0B2（冰箱压缩）,0B19（同济科技），0B20（华东电脑），0B9（长安信息），0B23（中科健）,0B16（湘计算机）,0B30（倍特高新），0B6（南华西中0B28（深圳华强），0B7（厦门信达），0B25（华意压缩），0B27（粤TC1.）,0B3（复华实业）,0B5（南洋实业）,0B13（中国高新），0B24（深华源）。第一、二类公司在经营规模、经营实力、技术水平等方面具有一定优势，竟争能力强，经营业绩优良，综合财务状况良好，属高科技板块的绩优龙头股，颇具发展潜力和长期投资价值，是高科技板块中投资者的首选投资对象。其中第一类的天津磁卡、燕化高新1997年净利润率分别为45.86%和44.32%,净利润增长率分别是93.6%和96.95%,其获利能力远远大于其它公司。第四类公司业绩一般,投资者应谨慎介入，可多加观望。（二）聚类分析方法在社会的应用由上述案例可以看出聚类分析能综合多项财务指标来反映上市公司的盈利能力和水平，所得聚类结果与公司的实际财务状况和经营状况相吻合。我们还可以对所选出的各个板块的龙头潜力股再进行聚类分析，找出最具实力的板块龙头股。因此聚类分析方法适用于分析社会上的一些公司的盈利能力和水平，在经济类比较突出。四、主成份分析方法（一）主成份分析方法原理简介主成分分析是将多指标化为少数几个综合指标的一种统计方法主成分分析是从原始变量中导出少数几个主分量，使他们尽可能多地保留原始变量的信息，且彼此互不相关主成分分析的应用目的是数据的压缩、数据的解释，它常被用来寻找判断某种事物或现象的综合指标，并且给综合指标所包含的信息以适当的解释，从而更加深刻地揭示事物的内在规律。（二）主成份分析案例叙述分析1.中学生身体四项指标的主成分分析在某中学随机抽取某年级30名学生，测量起身高（XI）,体重（X2）,胸围（X3）和坐高（X4）,数据如下表。试对这30名中学生身体四项指标数据做主成分分析。IXlII1148417278161523573792139347176171494782793160497786181453570774149366779191604774875159458086201564478856142316676211514273827153437683221473873788150437779231573968809151427780241473065751013931687425157488088111402964742615136748012161477884271443668761315849788328141306776141403367772913932687315137316673301483870782.对数据的相关阵作主成分分析，有> pr.stud<-princomp(student,cor=TRUE)> summary(pr.stud,Ioadings=TRUE)Importanceofcomponents:Comp.1Comp.2Comp.3Comp.4Standarddeviation1.88178050.559806360.281795940.25711844ProportionofVariance0.88527450.078345790.019852240.01652747CumulativeProportion0.88527450.963620290.983472531.000000001.oadings:Comp.1Comp.2Comp.3Comp.4Xl-0.4970.543-0.4500.506X2-0.515-0.210-0.462-0.6910.1750.4610.744-0.232X3-0.481-0.725X4-0.5070.368其中StandarddeViatiOn为主成分的标准差，即方差的开方，也就是相应的特征值的开方。ProportionofVariane表示方差的贡献率，而CUmUlatiVeProportion表示方差的累计贡献率。1.OadingS=FA1.SE或缺省就不列出loadings。3.分析：从主成分分析结果可看出前两个主成分的累计贡献率高达96%,选择两个主成分。第一个主成分对应系数的符号都相同，其值在0.5左右，反映了中学生身材的魁梧程度，身材高大的学生，他的四个部分的尺寸都比较大，因此第一主成分的值就较小。而身材矮小的同学他的四部分都比较小,第一主成分的值较大。第一主成分为大小因子。第二主成分是高度和围度之差，比较大表明该学生细高，比较小为“矮胖”，称第二因子为形体因子。看一下各样本的主成份值画第一个主成分的散点图，可看出10,11,15,29值较大,说明学生比较瘦小,而3,5,25值较小,说明学生比较高大.predict(pr.stud)->scoreplot(1:30,score,1)plot(1:30,score,2)1:301:30从这个图很容易看出，那些学生属于高大魁梧型，比如25号学生,3、5号学生,那些学生属于高瘦型比如23、19、4等等。（）主成份分析方法在社会的应用根据主成分分析的定义及性质，我们已大体上能看出主成分分析的一些应用。概括起来说，主成分分析主要有以下几方面的应用。1.主成分分析能降低所研究的数据空间的维数。即用研究m维的Y空间代替P维的X空间（InVP）,而低维的Y空间代替高维的X空间所损失的信息很少。即：使只有一个主成分Yl（即11=l）时,这个Yl仍是使用全部X变量（P个）得到的。例如要计算YI的均值也得使用全部X的均值。在所选的前m个主成分中，如果某个Xi的系数全部近似于零的话，就可以把这个Xi删除，这也是一种删除多余变量的方法。2 .有时可通过因子负荷aij的结构，弄清X变量间的某些关系。3 .多维数据的一种图形表示方法。我们知道当维数大于3时便不能画出几何图形，多元统计研究的问题大都多于3个变量。要把研究的问题用图形表示出来是不可能的。然而，经过主成分分析后，我们可以选取前两个主成分或其中某两个主成分，根据主成分的得分，画出n个样品在二维平面上的分布况，由图形可直观地看出各样品在主分量中的地位。4 .由主成分分析法构造回归模型。即把各主成分作为新自变量代替原来自变量X做回归分析。5 .用主成分分析筛选回归变量。回归变量的选择有着重的实际意义，为了使模型本身易于做结构分析、控制和预报，好从原始变量所构成的子集合中选择最佳变量，构成最佳变量集合。用主成分分析筛选变量，可以用较少的计算量来选择量,获得选择最佳变量子集合的效果。五、因子分析方法（一）因子分析方法原理简述因子分析在某种程度上可以被看成是主成分分析的推广和发展，它对问题的研究更加深入，研究相关阵或协方差阵的内部依赖关系，它将多个变量综合为少数几个因子，以再现原始变量与因子之间的相关关系，也是多元统计分析中降维的一种方法。因子分析是通过研究多个变量间相关系数矩阵的内部依赖关系，找出能综合所有变量的少数几个随机变量，这几个随机变量是不可测量的，通常称为因子。然它们在性质上的亲疏程度加以分类.实质是按照距离的远近将数据分为若干个类别，以使得类别内数据的差异尽可能小，类别间的“差异”尽可能大。（二）因子分析案例叙述分析众所周知，建立套利定价模型的关键在于因素的筛选，计算量很大。而因子分析能将为数众多的原始指标变量经过分析综合为少数几个公共因子变量，从而大大减少计算的复杂度。本文利用因子分析的方法对11个因素进行筛选，确定四个能够很好地反映所有因素包含的信息但又互不相关的公共因子变量，并建立套利定价模型，实证检验说明，通过该方法进行因素筛选建立的套利定价模型具有较好的定价效果。1 .问题的提出1976年，StePhenRoSS提出了著名的资产定价模型套利定价理论（ArbitragePricingTheory,APT）o该理论假设任何风险证券的收益率受K个因素的影响，由一个K因素线性模型给出：ri=ai+kk=lbikfk+i,i=l,2,n(l)其中：E(i)=E(fk)=E(ij)=E(ifk)=0；E(2i)=s2i<S2；ri为第i种风险证券的收益率；ai表示所有影响风险证券收益率的因素都为零时风险证券i的平均收益率;fk表示第k个因素的值;bik表示风险证券i对第k个因素的敏感性；i为随机扰动项。当不存在渐进套利机会时，由K因素线性模型可以得到如下的近似定价模型套利定价模型(APT)：E(ri)=ai0+Kk=lbikk(2)其中，k称为风险证券i对第k个因素的风险溢价。如果将误差记为Vi三ai-O-Kk=Ibikk,则当不存在渐进套利机会时，有limnlnni=lv2i=0.建立套利定价模型的关键在于因素的筛选。然而，一种风险证券的收益率受多方面因素的影响，同时我们也不知道究竟需要多少个因素来构造APT.假设有n个因素对证券的收益率有影响，则可能nm=lCmn种因素的组合。要从如此众多的因素组合中筛选出最优的因素组合，其计算量可想而知。一般来说，因子的辨识和确定有两种基本的方法：统计方法和推理方法。统计方法涉及从一个全面的资产收益集(通常远超过用来估计和检验的样本资产收益)来确定因子，采用这些收益的样本数据来构造表示因子的资产组合，如ConnorKorajczyk(1988)>1.ehmannftModest(1988),前者使用因子分析方法，后者利用主元分析方法。推理方法是基于捕捉经济的系统风险原则来辨识因子的，例如Fama和FrenCh(1988,1996)使用公司特征来构建因子资产组合。在这类研究中，在将股票分组后，对每一组股票首先采用因子分析方法来估计影响股票收益率的因子数目，并估计每只股票的因子载荷；然后，利用股票收益率数据和已估计出的因子载荷做横截面回归，估计因子的风险溢价，进而检验多因子模型的适用性。此外，由于APT认为股票收益率的风险可以分为可分散风险和不可分散风险,其中可分散风险部分的均值为零，在大样本中可忽略不计,而不可分散风险部分由K个共同因子决定，并通过K个因子系数反映股票收益率与每个非零风险溢价之间的关系。但是，现实中可能某一变量本身与不可分散风险不相关(即不应当作为因子出现)，但在APT模型中却被不恰当定价，成为一个共同因子。虽然由实际数据生成的因子模型通过了显著性检验，但却无法肯定这些因子就是不可分散风险的溢价，也无法排除可分散风险成为共同因子的可能。鉴于此，我们有必要对APT进行“自方差”检验这里用“自方差”只是一种强调性说法，其实质就是该项资产收益率的方差从长期来看，证券收益率的自方差与收益率均值之间总是保持很高的相关性，而自方差又是每一种证券所特有的，属于可分散风险。如果APT有效，那么单个证券的自方差就不应当对期望收益率起作用，因为APT认为只有不可分散的风险才对定价起作用，才可以成为定价因子。“自方差”检验就是要证明单个股票收益率的自方差是否为共同因子，可否用于定价，要接受还是否定APT.鉴于此，他们也利用“自方差”检验来对多因子模型做了补充研究。到目前为止，我国在套利定价理论因素确定方面的研究并不多，主要是利用多元线性回归构造套利定价模型，这一方法的计算量大，其包含的因素要么过多要么不全面，而且因素之间的关联程度较高。而因子分析是一种常用的统计降维技术，能够利用原始指标变量中某些指标之间的相关性对多变量的面板数据进行最佳综合和简化，将为数众多的指标综合为少数几个公共因子，以较少的几个公共因子变量反映原始指标变量的大部分信息，从而大大降低了分析问题的难度。2 .用因子分析法确定APT中的因素组合在已有的研究中，一般认为APT中至少包含有三类不同的因素：反映总体经济活动的指标、通货膨胀率以及某些类型的利率因素。鉴于此，本文将国民生产总值、工业生产总值、第二产业生产总值、第三产业生产总值、全国居民消费水平、通货膨胀率、全社会固定资产投资增长速度、社会消费品零售总额、货币供应总量、年净出口贸易总额、利率期限结构等11个因素作为原始指标变量，利用我国1980年到200=3年统计年鉴中的数据进行因子分析。（1）对原始指标变量进行相关性分析因子分析从众多的原始指标变量中构造出少数几个具有代表意义的公共因子变量，它要求原始指标变量之间要具有比较强的相关性，否则就无法从中综合出能反映某些变量共同特性的少数公共因子变量来，原始指标变量就不适于进行因子分析。因此，在因子分析之前需要对原始指标变量进行相关性分析。本文采用的是KMO（Kaiser-Meyer-Olkin）检验和Bartlett球度检验。统计量KMO的值为0.771,大于0.6,根据统计学家KaiSer给出的标准，原始指标变量适合做因子分析；BartIett球度检验给出的自由度为55的卡方近似值为780.924,相伴概率为0.000,小于显著水平0.05,适合进行因子分析。由KMC）检验和BartIett球度检验结果可知，原始指标变量适合做因子分析。与此同时，本文还对11个原始指标变量进行了反映像相关矩阵检验，在反映像相关矩阵中，所有偏相关系数的绝对值均小于0.05,说明所有的原始指标变量都适于进行因子分析。（2）构造公共因子变量构造公共因子变量是因子分析的一个核心问题。因子分析中确定公共因子变量的方法很多，本文采取的是主元分析法。确定保留公共因子变量的数目根据公共因子变量与其特征值的散点图（图1）可以看出，前面4个公共因子变量的特征值变化非常明显，从8.744到0.126,而从第5个公共因子变量开始，特征值的变化趋于平稳。这说明提取前4个公共因子变量对原始指标变量的信息描述有显著作用。为了能够得到更