数学建模之因子分析法.docx

《数学建模之因子分析法.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯因子分析因子分析就是一种降维、简化数据的技术。它通过研究众多变量之间的内部依赖关系，探求观测数据中的基本结构，并用少数几个“抽象”的变量来表示其基本的数据结构。这几个抽象的变量被称作“因子”，能反映原来众多变量的主要信息。原始的变量是可观测的显在变量，而因子一般是不可观测的潜在变量。1.因子分析法的应用①汽车行业业绩评价研究（下载文档），②上市公司盈利能力及资本结构实证分析，③生育率影响因素分析。2.步骤①对原始数据进行标准化处理用x1,x2,,xm表示因子分析指标的m个变量，评

2、价对象有n个，aij表示第i个评价对象对应于第j个指标的取值。将每个指标值aij转化为标准化指标aij，即aijaijj,(i1,2,,n;j1,2,,m)sj1式中：jnn1nj)2aij，sj(aiji1n1i1相应地，标准化指标变量为xjxjj,(j1,2,,m)sj②计算相关系数矩阵RR(rij)mmnakiakjk1,(i,j1,2,,m)rijn1式中：rii1,rijrji，rij是第i个指标和第j指标之间的相关系数。1⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯③计算初等载荷矩阵解特征方程IR0，得到特征

3、值i(i1,2,,m)12,m0，再求出相对应的特征值i的特征向量ui(i1,2,,m)，其中uj(u1j,u2j,,umj)T,得到初等载荷矩阵为11u1,,mum④确定主因子的个数k(km)km一般选取使得累计贡献率ii85%的这k个主因子，对k个因子载荷i1i1矩阵作旋转，用1(k)表示1的前k列，T表示正交矩阵，则得矩阵21(k)T，建立因子模型，即x111F11kFk,xmm1F1mkFk.⑥计算因子得分，作出综合评价求出单个因子的得分函数??个样本对第j个因子的得分估Fj，用Fij表示第i计值，Y表示原始数据标准化后的矩阵，则总得分为??)nkYR1kF(

4、Fij例题我国上市公司赢利能力与资本结构的实证分析已知上市公司的数据见表1表1上市公司数据试用因子分析法对上述企业进行综合评价。2⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯模型的建立①对原始数据进行标准化处理用12m表示因子分析指标的m个变量，评价对象有n个，a表示第ix,x,,xij个评价对象对应于第j个指标的取值。将每个指标值aij转化为标准化指标aij，即aijj,(i1,2,,n;j1,2,,m)aijsj式中：1jnn1nj)2aij，sj(aiji1n1i1相应地，标准化指标变量为xjxjj,(j1,2,,

5、m)sj②计算相关系数矩阵RR(rij)mmnakiakjk1,(i,j1,2,,m)rijn1式中：rii1,rijrji，rij是第i个指标和第j指标之间的相关系数。③计算初等载荷矩阵解特征方程IR0，得到特征值i(i1,2,,m)12,m0，再求出相对应的特征值i的特征向量ui(i1,2,,m)，其中uj(u1j,u2j,,umj)T,得到初等载荷矩阵为11u1,,mum④确定主因子的个数k(km)km一般选取使得累计贡献率ii85%的这k个主因子，对k个因子载荷i1i1矩阵作旋转，用1(k)表示1的前k列，T表示正交矩阵，则得矩阵21(k)T，建立因子模型，即

6、3⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯x111F11kFk,xmm1F1mkFk.模型的求解：我们选取两个主因子。利用MATLAB程序计算得旋转后的因子贡献及贡献率见表2，因子载荷阵见表3。表2贡献率数据表3旋转因子分析表计算因子得分，作出综合评价我们用回归方法求单个因子得分函数?bj1x1bjmxm,j1,2,,kFj用F?ij表示第i个样本对第j个因子的得分估计值，则?bj1xi1bjmxim,(i1,2,,n;j1,2,,k)Fij即b11b21bk1b12b22xk2R1Bb1mb2mxkm用Y表示原始数据

7、标准化后的矩阵，则总得分为??YR1BF(Fij)nk计算得出各个因子得分函数为F10.531x10.1615x20.1831x30.5015x4F20.045x10.5151x20.581x30.0199x4总得分为44.49F141.8F2F86.17计算出16家上市公司赢利能力的综合得分见表4。4⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯表416家上市公司赢利能力的综合得分5