多元統(tǒng)計(jì)正交因子分析實(shí)驗(yàn)報(bào)告

WordWord文檔正交因子分析［役計(jì)性賣(mài)驗(yàn)丿(Orthogonalfactoranalysis)實(shí)菠點(diǎn)理：因子分析是主成分分析的推/?和發(fā)展，其目的是用*數(shù)兒個(gè)不可觀測(cè)的隱變量，即因子，來(lái)解粹原始變量之間的和關(guān)關(guān)糸，它&是厲于多元分析中處理吟維的一種統(tǒng)計(jì)方凍。因子分析的基本思?想是通過(guò)變量間的協(xié)方差矩陣(戎相關(guān)糸數(shù)矩陣丿部結(jié)構(gòu)的研?死，尋找能技制所有變量的少救幾個(gè)因子去描述.多個(gè)變量之間的和關(guān)關(guān)糸。因子分析中最常用的數(shù)學(xué)樓更是正交因子棧型，其特點(diǎn)是棋型中的因子相互之間正交。下表中給岀了二戰(zhàn)以來(lái)典運(yùn)會(huì)運(yùn)動(dòng)員十頊運(yùn)動(dòng)成績(jī)的相關(guān)糸數(shù)矩陣：(E9a6)100卷1.00■?跳遺0.591.00??鉛球0.350.421.00?跳壽0.340.510.381.00400耒0.630.490.190.29110卷跨欄0.400.520.360.46鐵餅0.280.310.730.27擇竿跳壽0.200.360.240.39標(biāo)槍0.110.210.440.171500*-0.070.09-0.080」8

1.00 TOC\o"1-5"\h\z0.341.00 . ^ ^ .0」70.321.00 . ^ .0.230.330.241.00 . .0.130.180340.241.00 .0.390.00-0.020.17-0.001.00賣(mài)絵要求:(1J試由相關(guān)糸數(shù)矩陣作因子分析；covmat(2)試根據(jù)因子戟椅，并結(jié)合題目背景知識(shí)，對(duì)公共因子進(jìn)行命名。賣(mài)檢題目二：下在中給出了不同國(guó)家及地區(qū)的女子徑躱記錄：(t1a7)100m200m400m800m1500m3000mMarathonCountry(s)(min)(min)(min)(min)argentin11.6122.9454.52」54.439.79178.52australi11.222.3551.081.989.08152.37austria11.4323.0950.621.994.229.34159.37belgium11.4123.045224」48.88157.85bermuda11.4623.055332」64.589.81169.98brazil11.3123.1752.82」4.499.77168.75burma12.1424.47552」84.459.51191.02Canada1122.2550.0624.068.81149.45chile1224.5254.92.054.239.37171.38china11.9524.4154.972.084.339.31168.48Columbia11.62453.262」14.359.46165.42cookis12.927」60.42.34.8411」233.22costa11.9624.658.252.214.6810.43171.8czech11.0921.9747.991.894」48.92158.85denmark11.4223.5253.62.034」88.71151.75domrep11.7924.0556.052.244.749.89203.88finland223950.142.034」8.92154.23france11」522.5951.7324」48.98155.27gdr10.8121.7148.161.933.968.75157.68他11.0122.3949.751.954.038.59148.53gbni1122.1350.461.984.038.62149.72greece11.7924.0854.932.074.359.87182.2guatemal11.8424.5456.092.284.8610.54215.08hungary11.4523.0651.52.014」48.98156.37india11.9524.2853.62」4.329.98188.03indonesi11.8524.2455.342.224.6110.02201.28Ireland11.4323.5153.242.054.118.89149.38israel11.4523.5754.92」4.259.37160.48italy11.292352.011.963.988.63151.82japan11.732453.732.094.359.2150.5kenya11.7323.8852.724」59.2181.05korea11.9624.4955.72」54.429.62164.65dprkorea12.2525.7851.21.974.259.35179.17luxembou12.0324.9656.12.074389.64174.68malaysia12.2324.2155.092」94.6910.46182.17mauritiu11.7625.0858.12.274.7910.9261.13mexico11.8923.6253.762.044.259.59158.53

netherla11.2522.8152.381.994.069.01152.48nz11.5523.1351.62.024.188.76145.48norway11.5823.3153.122.034.018.53145.48png12.2525.0756.962.244.8410.69233philippi11.7623.5454.62」94.610.16200.37poland11.1322.2149.291.953.998.97160.82Portugal11.8124.2254.32.094」68.84151.2rumania11.4423.4651.21.923.968.53165.45singapor12.32555.082」24.529.94182.77spain11.823.9853.592.054」49.02162.6Sweden11.1622.8251.792.024」28.84154.48switzerl11.4523.3153.112.024.078.77153.42taipei11.2222.6252.52」4.389.63177.87thailand11.7524.4655.82.24.7210.28168.45turkey11.9824.4456.452」54.37938201.08usa10.7921.8350.621.963.958.5142.72ussr11.0622.1949.191.893.878.45151.22wsamoa12.7425.8558.732.335.8113.04306（數(shù)據(jù)來(lái)源：1984年洛杉機(jī)典運(yùn)會(huì)IAAF/AFT徳球與田球統(tǒng)計(jì)手冊(cè)丿ussr11.0622.1949.191.893.878.45151.22rumania11.4423.4651.21.923.968.53165.45賣(mài)絵要求:

(1丿很據(jù)以上數(shù)據(jù)對(duì)女子後賽頊目作因子分析；對(duì)亦共因孑進(jìn)行解釋?zhuān)挥?jì)算各個(gè)國(guó)家的笫一因子得分并進(jìn)行排名。要求列出排名前10的國(guó)家或地區(qū)，并給岀中國(guó)的名次。賣(mài)蠢題目一分析掖告：R?4：record<-read.table(Hdata4.txf\head=F) #導(dǎo)入救據(jù)record<-record[,-1] #刪除笫一列record<-as.matrix(record) #將原救據(jù)矩陣化option$(digit$=2) #保Q兩住小教pca.datal<-princomp(covmat=record)#以和關(guān)糸數(shù)矩陣作為基礎(chǔ)，建立主成分分析summary(pca.datal) #輸出主成分分析報(bào)表factl.$t<-factanal(covmat=record,factor$=5,rotation=”none”)#作因子分析,不淡轉(zhuǎn)factl.stfactl.rofactl.ro<-factan8l(covm8t=recorcLfdctors=5、rotdtion=\8rimax”)#作因子分析,淡攜factl.stfactl.ro#輸出不炎轉(zhuǎn)的結(jié)果#輸出淡轉(zhuǎn)的結(jié)黑#計(jì)算共同度apply((factl.ro$loadings)八2J』um)#計(jì)算共同度f(wàn)act2.ro<-factan8l(covm8t=recorcLfdctors=4、rot8tion=”vdrimdx”)#作因孑分析,淡轉(zhuǎn)fact2.ro#輸出淡轉(zhuǎn)的結(jié)果#計(jì)算共同度apply((fact2.ro$loadings)^2J,sum)#計(jì)算共同度輸岀姑系及分析：(1丿試由柑關(guān)糸數(shù)矩陣作因子分析；record<-read.table(ndata4.txt\head=F)#導(dǎo)入數(shù)據(jù)record<-record[,-1]#刪徐第一列record<-as.matrix(record)#將療教據(jù)矩陣化option$(digit$=2)#保倒兩住小數(shù)pca.datal<-princomp(covmat=record)#以和關(guān)糸數(shù)矩陣作為基礎(chǔ)，建立主成分分析summary(pca.datal) #輸出主成分分析圾表為了確主因子分析中因子的數(shù)目,我們先對(duì)柑關(guān)糸數(shù)矩陣做主成分分析叔1主成分分析掖在Comp.1Comp.2Comp.3Comp.4Comp.5Comp.6Comp.7Comp.8Comp.9Comp」0Standarddeviation1.951.231.060.9560.8490.7710.7260.6190.4850.456ProportionofVariance0.380.150.110.0910.0720.0590.0530.0380.0240.021CumulativeProportion0.380.530.640.7330.8050.8650.9170.9560.9791.000

由方差累計(jì)賈伙率得刊，柱笫五主成分，累積賈故率達(dá)到了80%以上，并趨于稔定。我們確定因子分析中因子數(shù)目為5.factl.st<-factanal(covmat=record,factor$=5,rotation=MnoneM)#作因子分析,不淡轉(zhuǎn)factl.stfactl.rofactl.ro<-factanal(covmat=record,factor$=5,rotation=HvarimaxM)#factl.stfactl.ro#輸出不炎轉(zhuǎn)的結(jié)果#輸出淡轉(zhuǎn)的結(jié)果apply((factl.ro$loadings)^2J,sum) #計(jì)算共同度做因子分析，得到未淡轉(zhuǎn)的因子栽持以及淡轉(zhuǎn)的因子栽希表2未炎轉(zhuǎn)的因孑載荷FactorFactorlFactor2Factor3Factor4Factor5100耒0.2080.7910.301-0.167跳迄0.3780.5950.2460.242鉛球0.6440.761跳需0.41503440.1570.471-0.139400*0.4460.688-0.113-0.2030.116110耒跨欄0.2650.4350.2610.343鐵餅0.5030.534撐竿跳需0.3070.2400.4020.214標(biāo)槍0.3130.引40.3781500*0.707-0.704

累積貢伙率0.20.380.550.6160.640表3炎轉(zhuǎn)的因孑我椅FactorFactorlFactor2Factor3Factor4Factor5Communalities100*0.1710.8150.276-0.1410.79跳遺0.2230.4800.5800.62鉛球0.9550.1390.2411.00跳需0.2110.1520.6870.1170.56400*0.7600.19303260.1260.74110耒跨欄0.1870.2780.5650.45鐵餅0.6930.1250.1940.1110.55撐竿跳壽0.1120.5010.1190.2820.36標(biāo)槍0.4080.1400.4010.351500卷0.9891.00累積貢故率0」70.340.500.610.640觀疼表格中枚標(biāo)注為綠色的兩個(gè)因子我満(標(biāo)槍項(xiàng)目一行丿,在Factor!中的因子栽'荷為0.408,A.Factor5中的因子我持％0.401,比較兩個(gè)因子我椅，0.408>0.401,因此我們最終選取0.408。這樣一來(lái)，我們做因子分肘，只謝要4個(gè)因子即可。因此.我們下面冉做4個(gè)因子的簸桔因子分析。fact2.rov?factanal(covmat=record,factor$=4,rotation=HvarimaxM)#作因孑分析,淡轉(zhuǎn)

fact2.ro#輸出淡轉(zhuǎn)的結(jié)果fact2.roapply((fact2.ro$loadings)八2J,sum) #計(jì)算共同度在4<44的因孑我椅FactorFactorlFactor2Factor3Factor4Communalities100耒0.1670.8570.246-0.1380.84跳遠(yuǎn)0.2390.4760.5810.62鉛球0.9630.1530.2011.00跳需0.2420.1720.6320.1130.50400耒0.7100.23603310.67110卷跨欄0.2050.2610.5880.46鐵餅0.6990.1330.1790.54撐華跳壽0.1380.5120.1170.30標(biāo)槍0.4180.1750.211500耒0.1130.9881.00累積賈秋率0」80340.500.61(2)試根據(jù)因子戟椅，并結(jié)今題目背景知幟，對(duì)公共因子進(jìn)行令名由淡轉(zhuǎn)后的栽持可發(fā)現(xiàn)，第一因子中，鉛球、鐵餅和標(biāo)槍的栽持較大，可令名為儀擲因子；笫二因子中，100卷和400耒的我椅較大，可命名為短跑因子；笫三因子中，跳迄、跳需、110卷跨欄、擇竿跳壽較丸，可令名為彈跳因子；第四因子中，1500耒的我荷較丸，可令名為長(zhǎng)跑因子。

賣(mài)絵題目二分析掖告:R程序:bv-read.csvfdata42.csvv)#導(dǎo)入數(shù)據(jù)b1<-b[,-l]#刪除笫一列pc.b1<-princomp(b1,cor=T)#做主成分分析summary(pc.bl)#主成分分析結(jié)果fact.b1<-factanal(bl,factor=2,method=umleM,rotation=HnoneH)#未淡轉(zhuǎn)的因子分析fact.b1$loadings#輸出不炎轉(zhuǎn)的結(jié)果fact.b2<-factanal(bljactor=2,method=vmleM,rotation=Mvarimax,\scores=Hregression11)#淡轉(zhuǎn)的因子分析fact.b2$loadings#輸出炎轉(zhuǎn)的結(jié)果apply((fact.b2$loadings)2J,sum)#計(jì)算共同度shapiro.test(fact.b2$scores)#檢驗(yàn)正態(tài)性fact.b3<-factanal(bLfactor=2,method=,,mle,\rotation=Mvarimax,\scores=MBartlett1')b[order(fact.b3$scores[J],decreasing=F)J] #排名輸岀結(jié)果及分析：(1J根據(jù)以上數(shù)據(jù)對(duì)女子徑瘵項(xiàng)目作因子分析；bv-read.csvfdata42.c$vu) #4■入數(shù)據(jù)bl<-b[,-1]pc.b1<?princomp(b1,cor=T)summary(pc.bl)裹4主成分分析結(jié)果CompCompCompCompCompCompComp

.1.2.3.4.5.6.7Standard0.2320」970.1492.410.8080.5480.354deviation068Proportionof0.0070.0050.0030.830.0930.0430.018Variance762Cumulative0.9910.9961.0000.830.9230.9660.984Proportion280根據(jù)主成分分析的結(jié)果可以看出，在笫2個(gè)特征根處，累計(jì)賈伙率就已經(jīng)達(dá)到了92.3%。因此，我們選用2個(gè)因孑進(jìn)行因孑分析。fact.b1<-factanal(b1,factor=2,method=nmleM,rotation=”none”)fact.bl$loadingsfact.b2<-factanal(blJactor=2,method=nmleH,rotation=MvarimaxK,scores=nregression")fact.b2$loadingsapply((fact.b2$loadings)2J,sum)表5未炎轉(zhuǎn)的因孑我椅Factor!Factor2X100.m..s.0.95-0.13X200.m..$.0.97-0.22X400.m..$.0.900

X800.m..min.0.830.38X1500.m..min.0.840.53X3000.m..min.0.840.49Marathon..min?0.800.40表6<44的因孑我椅Factor!Factor2CommunalitiesX100.m..