数学建模多元统计分析

1、精选文库精选文库实验报告一、实验名称多元统计分析作业题。二、实验目的（一）了解并掌握主成分分析与因子分析的基本原理和简单解法。（二）学会使用matlab 编写程序进行因子分析，求得特征值、特征向量、载荷矩阵等值。（三）学会使用排序、元胞数组、图像表示最后的结果，使结果更加直观。三、 实验内容与要求精选文库四、实验原理与步骤（一）第一题：1、实验原理：因子分析简介：(1) 1.1 基本因子分析模型设 p 维总体 x=(x1,x2,.,xp)' 的均值为 u=(u1,u2,.,u3)' ，因子分析的一般模型为x1=u1+a11f1+a12f2+.+a1mfm+1x2=u2+a21f

2、1+a22f2+.+a2mfm+ 2.xp=up+ap1f1+fp2f2+.+apmfm+ p其中， f1,f2,.,fm 为 m 个公共因子；i 是变量 xi(i=1,2,.,p) 所独有的特殊因子，他们都是不可观测的隐变量。称aij(i=1,2,.,p;j=1,2,.,m) 为变量 xi 的公共因子fi 上的载荷，它反映了公共因子对变量的重要程度，对解释公共因子具有重要的作用。上式可以写为矩阵形式x=u+Af+ 精选文库其中A=(aij)pxm称为因子载荷矩阵；f=(f1,f2,.,fm)' 为公共因子向量；=( 1, 2,. p)称为特殊因子向量(2) 1.2 共性方差与特殊方差

3、xi 的方差 var(xi) 由两部分组成，一个是公共因子对xi 方差的贡献，称为共性方差；一个是特殊因子对xi 方差的贡献，称为特殊方差。每个原始变量的方差都被分成了共性方差和特殊方差两部分。(3) 1.3 因子旋转因子分析的主要目的是对公共因子给出符合实际意义的合理解释，解释的依据就是因子载荷阵的个列元素的取值。当因子载荷阵某一列上各元素的绝对值差距较大时，并且绝对值大的元素较少时，则该公共因子就易于解释，反之，公共因子的解释就比较困难。此时可以考虑对因子和因子载荷进行旋转（例如正交旋转），使得旋转后的因子载荷阵的各列元素的绝对值尽可能量两极分化，这样就使得因子的解释变得容易。因子旋转方法

4、有正交旋转和斜交旋转两种，这里只介绍一种普遍使用的正交旋转法：最大方差旋转。 这种旋转方法的目的是使因子载荷阵每列上的各元素的绝对值（或平方值） 尽可能地向两极分化，即少数元素的绝对值（或平方值）取尽可能大的值，而其他元素尽量接近于 0.(4) 1.4 因子得分在对公共因子做出合理解释后，有时还需要求出各观测所对应的各个公共因子的得分，就比如我们知道某个女孩是一个美女，可能很多人更关心该给她的脸蛋、身材等各打多少分，常用的求因子得分的方法有加权最小二乘法和回归法。注意：因子载荷矩阵和得分矩阵的区别：因子载荷矩阵是各个原始变量的因子表达式的系数， 表达提取的公因子对原始变量的影响程度。因子得分矩

5、阵表示各项指标变量与提取的公因子之间的关系，在某一公因子上得分高，表明该指标与该公因子之间关系越密切。简单说，通过因子载荷矩阵可以得到原始指标变量的线性组合，如 X1=a11*F1+a12*F2+a13*F3, 其中 X1 为指标变量 1， a11、 a12、a13 分别为与变量 X1 在同一行的因子载荷， F1、 F2、 F3 分别为提取的公因子；通过因子得分矩阵可以得到公因子的线性组合，如F1=a11*X1+a21*X2+a31*X3，字母代表的意义同上。(5) 1.5 因子分析中的 Heywood （海伍德）现象如果 x 的各个分量都已经标准化了，则其方差=1。即共性方差与特殊方差的和为

6、1。也就是说共性方差与特殊方差均大于0，并且小于1。但在实际进行参数估计的时候，共性方差精选文库的估计可能会等于或超过1，如果等于1，就称之为海伍德现象，如果超过1，称之为超海伍德线性。 超海伍德现象意味着某些特殊因子的方差为负，表明肯定存在问题。造成这种现象的可能原因包括：共性方差本身估计的问题；太多的共性因子，出现了过拟合；太少的共性因子，造成拟合不足；数据太少，不能提供稳定的估计；因子模型不适合这些数据。当出现海伍德现象或超海伍德现象时，应对估计结果保持谨慎态度。可以尝试增加数据量，或改变公共因子数目， 让公共因子数目在一个允许的范围内变动， 观察估计结果是否有改观；还可以尝试用其他多元

7、统计方法进行分析，比如主成分析。2、实验步骤：(1) 将原始数据标准化处理 ;(2) 建立相关系数矩阵并计算其特征值和特征向量;将题目所给的相关系数矩阵输入并设为PHO ，利用x,y=eig(PHO) 求得 PHO 的特征值和特征向量。(3) 选择特征值大于等于1 的特征值个数为公共因子数,或者根据特征值累计贡献率大于80%来确定公共因子。从 y 表（特征值表） 中我们可以看出大于等于1 的特征值个数为4，所以公共因子数也为4。(4) 求得正交或斜交因子载荷矩阵;通过上一步，我们得到了因子载荷矩阵lambda。因子分析，公共因子数为4，设置特殊方差的下限为0，使用factoran 函数进行因子

8、旋转。设置表头与变量名，计算贡献率与累计贡献率。将lambda、 Contribut 、 CumCont放在一起，转为元胞数组，并显示最后的结果。(5) 计算公因子得分和综合得分。计算因子得分方法是：用每个共因子的方差贡献率做权数，对每个因子进行加权，然后加总得到每个项目的总因子得分按总得分的多少进行排序，以反映不同免死方面对结果的影响。精选文库根据第3 步我们可以得到因子14 的贡献率分别为：29.112927.330310.76447.2713，即为它们的权重。再由不同面试方面所对应的因子数我们可以得到公因子得分和综合得分。（二）第二题：1、实验原理：主成分分析是由皮尔逊在1901 年首先

9、对非随机变量引入的，后来由霍特林在1933 年推广到随机向量的情形。主成分分析是一种通过降维技术把多个变量化为少数几个主成分(即综合变量 )的多元统计方法，这些主成分能够反映原始变量的大部分信息，且新主成分之间相关性低、重叠少。主成分分析在很多领域都有广泛的应用，一般来说， 当研究的问题涉及多个变量，并且变量间相关性明显，即包含的信息有所重叠时，可以考虑用主成分分析的方法，这样更容易抓住事物的主要矛盾，使问题简化。具体做法是：对原始数据进行标准化处理用 x1 , x2 ,L , xm 表示主成分分析指标的m 个变量， 评价对象有n 个， aij 表示第 i 个评价对象对应于第 j 个指标的取值

10、。将每个指标值aij 转化为标准化指标a%ij ，即a%aijj ,( i1,2,L , n; j 1,2,L , m)ijsj式中： j1 naij ， sj1n(aijj )2n i 1n 1 i1相应地，标准化指标变量为%x jj,( j 1,2,L , m)xjsj计算相关系数矩阵RR (rij )m mn% %akiakjrijk 1,( i , j 1,2,L , m)n1精选文库其中： rii1,rijr ji ， rij 是第 i 个指标和第j 指标之间的相关系数。计算相关系数矩阵的特征值与特征向量解特征方程IR 0，得到特征值(1,2, )；再求出相iiLm12L, m0对应

11、的特征值 i的特征向量 ui (i1,2,L, m) ，其中 u j(u1 j , u2 j ,L, umj )T, 由特征向量组成的 m 个新的指标变量为y1%u11x1u21x2L um1 xmy2%u12x1u22 x2L um 2 xmMym%u1m x1u2 m x2L umm xm其中： y1 为第 1 主成分， y2 为第 1 主成分， ?， ym 为第 m 主成分选择 p（p m） 个主成分，计算综合评价值。（ 1）计算特征值j ( j1,2,L , m) 的信息贡献率和累积贡献率用 bj 表示主成分yi 的信息贡献率，则有bjj( j 1,2,L , m)mkk 1用 ap

12、表示主成分y1, y2 ,L , yp 的累积贡献率，则有pkk 1a pmkk 1若 ap 接近于 1（一般 ap 的范围为 85%95%）时，则用前 p 个指标变量y1 , y2,L , yp作为 p 个主成分，代替原来m 个指标变量，再对p 个主成分进行综合分析。计算综合得分用 b j 表示第 j 个主成分的信息贡献率，则有pZbj y jj1根据综合得分值进行评价。精选文库2、实验步骤：本题使用的是调用 pcacov 函数做主成分分析的方法。(1)首先我们定义相关系数矩阵PHO，并输入题目中所给的数据。(2)设置 pcacov 函数的返回，分别为：主成分表达式的系数矩阵x；相关系数矩阵

13、的特征值向量 y；主成分贡献率向量z。(3) 接下来就可以调用 pcacov 函数根据相关系数矩阵作主成分分析。(4) 为了使结果更加直观，我们选择用元胞数组形式（result1 ）显示结果。在横轴方向分别写出“特征值”、“差值”、“贡献率”和“累积贡献率”。其中，diff函数式用于求导数与差分； comsum函数用于计算一个数组各行的累加值。(5)从 result1我们可以看出， 前两个主成分累积功效率为 97.7680，因此我们用前两个主成分进行后续分析。(6)我们再次使用元胞数组形式（result2 ）显示前 3 个主成分表达式的系数矩阵，即z 矩阵的前 2 列，以使结果更加直观。(7)

14、 通过以上我们得到了 2 个主成分的公式以及主成分综合评价模型。(8)接下来利用公式>> df=PHO*x(:,1:num)和>> tf=df*z(1:num)/100;以及 >>stf,ind=sort(tf,'descend')进行综合得分和各个主成分的得分，再把得分按从高到低的次序排列，即可得到公司的综合排序。四、程序及其运行结果（一）第一题：1、程序1（特征值与特征向量的计算）：>> x,y=eig(PHO)运行结果1：x =Columns 1 through 100.0152-0.07770.11960.08190.17

15、54-0.00310.44030.1378 -0.14500.6108-0.01490.0797-0.22940.2167-0.1166-0.0039-0.05560.0003 -0.0059-0.0478精选文库-0.0616-0.0222-0.0813-0.1235-0.2265-0.2841-0.0169-0.0438-0.31040.23160.3235-0.1245-0.06270.00670.3866-0.2352-0.3941-0.28150.09740.11270.26710.3150-0.3168-0.36250.3504-0.23790.16290.3807-0.11280

16、.0054-0.4623-0.1824-0.35470.31070.0196-0.19790.1087-0.01420.51370.1688-0.1806-0.20840.2690-0.0078-0.01310.02380.45130.3024 -0.1115 -0.31030.12560.17620.62240.2938-0.0938-0.49500.0089-0.1606 -0.0167 -0.09000.0136-0.01650.00460.19870.0755-0.0537-0.30440.60080.1152-0.37320.2237-0.2027-0.32080.41730.007

17、80.25460.1339-0.2384-0.5057-0.2522-0.2654-0.55520.2292-0.33320.13030.1459-0.37890.0808-0.21330.10880.5406-0.12580.0899-0.1343-0.43570.35950.06030.09290.42080.1465-0.29650.57960.19740.03010.24460.5242-0.1141-0.12080.00620.0671-0.20860.2500-0.1729-0.1643-0.5916-0.1372-0.22390.0768-0.24470.0850-0.1151-

18、0.0288-0.0643-0.49610.0186-0.10370.2905-0.42900.1796-0.4213Columns 11 through 150.1181-0.11050.31060.43100.16310.87320.2562-0.0165-0.03500.2128-0.21010.6493-0.40820.23620.0398精选文库-0.11500.32470.4785-0.12640.22520.0003-0.1631-0.2471-0.24800.2913-0.2069-0.0630-0.1520-0.12920.3165-0.07030.40530.2974-0.

19、40530.15820.0344-0.2072-0.2049-0.03930.3221-0.10450.06760.07640.55330.1334-0.2039-0.1514-0.08250.04810.31490.1591-0.1905-0.2110-0.06850.3183-0.07380.0885-0.1110-0.02110.3314-0.12180.1986-0.06050.02260.3333-0.1170-0.21370.4606-0.07990.25930.0745-0.02630.08820.42050.2358y =Columns 1 through 100.034400

20、000000000.057700000000000.076300000000000.093300000000000.150900000000000.205700000000000.259700000000000.310300精选文库000000000.344100000000000.484000000000000000000000000000000000000000000000000000Columns 11 through 15000000000000000000000000000000000000000000000000000.7414000001.2091000001.467700000

21、2.0615000007.5040精选文库2. 程序 2（求因子载荷矩阵、贡献率和累计贡献率）：>>lambda,psi,T=factoran(PHO,4,'xtype','covariance','delta',0,'rotate','none')head = '变量 ', '因子 f1', '因子 f2','因子 f3','因子 f4'varname = ' 申请信的形式','外貌 ',&

22、#39;专业能力 ','讨人喜欢的能力','自信心 ','洞察力 ','诚实 ','推销能力 ','经验 ','驾驶汽车本领','抱负 ','理解能力 ','潜力能力 ','对工作要求强烈程度 ','适应性 ','<贡献率 >','<累积贡献率 >''Contribut = 100*sum(lambda.2)/15;CumCont =

23、 cumsum(Contribut);result1 = num2cell(lambda; Contribut; CumCont);result1 = head; varname, result1运行结果2：lambda =0.47000.06730.5437-0.22290.28000.46270.06090.1335-0.32000.43930.35910.22560.69000.11790.16520.53130.48000.7026-0.3997-0.06850.53000.7141-0.12440.05050.45000.1388-0.23270.63540.55000.7115-0

24、.1032-0.22570.21000.20980.6947-0.24230.61000.60620.0850-0.15850.55000.7099-0.1390-0.18060.55000.72530.08190.11320.54000.73300.20860.21471.0000-0.0000 -0.0000 -0.00000.40000.40780.5997-0.2296psi =0.42930.68600.52480.20040.11150.19120.32030.12970.37000.22810.14170.15190.0815精选文库0.00000.2614T =10000100

25、00100001result1 ='变量 ''因子f1''因子f2''因子f3'' 因子 f4''申请信的形式' 0.47000.06730.5437-0.2229' 外貌' 0.2800 0.4627 0.06090.1335'专业能力'-0.32000.43930.35910.2256'讨人喜欢的能力' 0.69000.11790.16520.5313' 自信心' 0.4800 0.7026 -0.3997-0.0685'

26、洞察力' 0.5300 0.7141 -0.12440.0505' 诚实' 0.4500 0.1388 -0.23270.6354' 推销能力 ' 0.5500 0.7115 -0.1032-0.2257' 经验' 0.2100 0.2098 0.6947-0.2423'驾驶汽车本领' 0.61000.60620.0850-0.1585' 抱负' 0.5500 0.7099 -0.1390精选文库-0.1806' 理解能力 ' 0.5500 0.7253 0.08190.1132'

27、潜力能力 ' 0.5400 0.7330 0.20860.2147' 对工作要求强烈程度 ' 1.0000 -4.5688e-06 -1.7616e-07 -3.2322e-07' 适应性' 0.4000 0.4078 0.5997-0.2296'<贡献率 >'29.112927.330310.76447.2713'<累积贡献率>'29.112956.443167.207574.47883. 运行结果 3：4 个因子的权重分别为：0.3901 0.3670 0.1445 0.0976所以对应的因子得

28、分与综合得分为（最后一列为综合得分）申请信的形式0.18330.02470.0786 -0.0218 0.2650外貌 0.1092 0.1698 0.0088 -0.02180.2660专业能力0.1248 0.1612 0.0519 0.2013 0.5392讨人喜欢的能力 0.26920.04330.02390.02200.3584自信心0.18720.2579-0.05700.05190.440洞察力 0.20680.2612-0.0180-0.00670.4433诚实 0.1755 0.0509 -0.0336 0.00490.1977精选文库推销能力0.21460.2611-0.0

29、180-0.00670.451经验 0.08190.07070.1004-0.02200.2310驾驶汽车本领 0.23800.22480.0123-0.02360.4987抱负0.21460.2605-0.0201-0.01760.4374理解能力0.21460.26620.01180.01100.5036潜力能力0.21070.26900.03010.02100.5308对工作要求强烈程度0.3901-0.3901适应性0.15600.14970.0867-0.02240.3700（二）第二题：1. 程序 1：>> clear allPHO=43.31 7.39 8.73 54

30、.8917.11 12.13 17.29 44.2521.11 6.03 7 89.3729.55 8.62 10.13 7311 8.41 11.83 25.2217.63 13.86 15.41 36.442.73 4.22 17.16 9.9629.11 5.44 6.09 56.2620.29 9.48 12.97 82.233.99 4.64 9.35 13.0422.64 11.13 14.3 50.51精选文库4.43 7.3 14.36 29.045.4 8.9 12.53 65.57.06 2.79 5.24 19.7919.82 10.53 18.55 42.047.26

31、2.99 6.99 22.72;>> PHO=zscore(PHO);>> r=corrcoef(PHO);>> x,y,z=pcacov(r)运行结果1：x =0.5306-0.4122-0.7018-0.23660.59380.4045-0.02290.69520.26060.7207-0.0097-0.64230.5458-0.38350.7119-0.2196y =1.89711.54960.39300.1602z =47.428238.74119.82534.00532.程序 2：>> result1(1,:)=' 特征值 &#

32、39;,' 差值 ',' 贡献率 ',' 累积贡献率 ' result1(2:5,1)=num2cell(y);result1(2:4,2) = num2cell(-diff(y);result1(2:5,3:4) = num2cell(z, cumsum(z)运行结果2：'特征值 ''差值''贡献率''累积贡献率 '1.89710.347547.428247.42821.54961.156638.741186.16930.39300.2328 9.825395.99470.1602

33、 4.0053 100.0000精选文库3. 程序 3：>> s=' 标准化变量 ''x1: 销售净利率 ''x2: 资产净利率 ''x3: 净资产收益率 ''x4: 销售毛利率 'result2(:,1)=s;result2(1, 2:3) = 'Prin1', 'Prin2'result2(2:5, 2:3) = num2cell(x(:,1:2)运行结果3：result2 ='标准化变量 ''Prin1''Prin2'

34、'x1:销售净利率 '0.5306-0.4122'x2:资产净利率 '0.5938 0.4045'x3:净资产收益率 '0.2606 0.7207'x4:销售毛利率 '0.5458-0.38354. 程序 4：>> df=PHO*x(:,1:num);>> tf=df*z(1:num)/100;>> stf,ind=sort(tf,'descend')运行结果4：stf =1.10911.10661.00290.75060.63400.36110.20140.0481-0.0570-0.2014-0.2444-0.3367-0.7301-0.8998-1.2651-1.4791ind =26151191341251738101614五、实验结果分析（一）第一题：精选文库1.从因子