实验报告-因子分析(多元统计)精选全文
精选全文整版可完编修改辑实验告报成分主分析综合(实性验()riPnicaplcmopnoetannlaysis)验原理实主成分分析:利用指标间的相之关性将多,个标指化转少数几个综为指标,合从达到降维和而数结构简据化目的。的些这合综标指映了原始反标指绝大的分信部,息常通表为原始指标示的种线性组某合,综且合标指间不相关利用。矩阵代的知数识可解主成分求。实验题目一:将彩胶色卷在显影液处下理在不同后情形曝光,下然后过通、红绿、三种滤蓝色并片高在中、、低种三度密下进测行,量每胶卷个有高、红高绿高、、蓝红中、…低、等蓝九指标个(分别记为X1-X九个变量)。试9了验10个胶卷8,由据数已算如得下协差阵:S2a1()17717959693532-74-3-419245113811271-2134026010941244111581204243217369456128228343133393948实验要:(1求试)从差阵出协进发主行分分成;析(2)算计差方积贡累献率(;)作Scree图3,并结(2合)的果结确主成分定的数;个(4)试对果结进行释。解1实验目题二下表:给出了不中国家同及地的区子径男赛记录:(ta68)Cuontry100m200m40m008m01500m500m10,0000Mamrahotn(s)(s)(s)(imn)mi(n)m(ni)(min)(mni)sArgneina10t3.920.816.8441.8137.140.429.3617.327Austarila10.1302.6440.841.743.5713.8227.6612.83uAsriat104.20.84146.281.793.6132.67272.135.9elBiumg0.34120.864.541.7303.613222..471592.95Berudam1.028205845..911.8375.1.46380.55146.62Brazil012.22.0344.5121.33.7616.63228.2613.313Burma10.4621.52483.18.38.54.154302.8319.59Caadna1.1702.2204.6851.67.63313.55820.9130.51hCle1i.0423.8046.1279.3.711.61329.13340.3Chia1n0.512104.473.18.137.3139.291.333.135oCumlbai1.4021305.6.1418.237.143492.7881.133.5CooIslkadns12.183.225.9422.024.24167.5.3836147.CsoaticRa01.42919.4.66818.7.38414.038.82113.685Czecohsloavika10.5302.564.5461.763.5813.428.12193.423eDmnrak1.0562.52045.81.7893.6131.52.811130.87DoinmiancRpeulicb01.41206.56.841.83.28214.191.43551412.Filnand10.3420.649.4591.73461.312727.5.21308.7Frnaec1.01120.385.2481.33.577313427.9.732.31GemranD.R.()10.120.33244.87173.3.5613.1727.2419.922Grema(nF.R).01.620.17344.51.733.353121.2.61172.332reatGBri.tN&I.elandr101120..2144.3917..51133.0127.1152.931Grecee1022.20.7146.651.873.641.4598245.31.64Guatealam0.9182.8124.84.891.8341.1630.1113.93Hu3gnary0.2162.6024.602.717362.3.1498244.132.58Idnia106.2.1244.753.761.731337.2878.1131.892IdonnesaIirleandsIrealItaylJaapnenyKKaroaDe.PRKo.raeuxemLburgoMalasiyaMuriaitsuMeicoxNetherlandseNwZelaadn
NorawyapuaNPewGinuaePilhippnesPilaodnProtuaglRmauinainSgaoperpainSSewendSwtizelrndaTapeiihTalaidnurTeykUASUSSR10.591.2447.89.814.3291.7341.6102096.643..7193.651.3230.172114.81.7773.271.3661.0011.729452.61.73.36132.31.0432.81045.861.7936.41341.104.620.664.4921.73.3553.1110.340.29846.91.937.77139.60.1921.149473..813577.14.1310.352.07774.1.483.2761364.1.4002.2946.318.3.821.64411.9122.5474.17883.83.15.0601.224.314.611.83.513646.015.220.595.1417.4.3261.3630.1512088.46.1.713.4451321.1.5502.11646711.7.636.21.3430.9162.1847.7199.401.14.27017.281.644624..1813.38147.410.160.2245.436.763.611.29310.3521.746.711.97.3623.13101.124.08945.817.67.64313.25103.21828.7.4148.83.895.1110.1224.7075.9841.7635.53131.0.2120.5146.561.737.6133.121903.72.4604.7581.87355.31.2210.952.29146.817.3.791740710.392..1094.7911.8338.4512.103.1271.434.761.97.3673.16959319.7.543.81.67335.313.2010.27044.6.15735913..230.7927.1288.397.5222772.27.382.93292.762908.1.0133.71277.952.6127.72776.3193.36.6402789.7.3822.77631.322.732779.47.92103.073.5262.85827.437253.14.83812.335137.5531.01188.623129.7531.65230.87141.172154.15212.13922.219.0212.89813.18418.42214.25131758.18.2561235.157.7731.57113.63031.1239.271149.931.5128.22130.5513WetsernamSo1a.8202.18649.2024.2416.832.417(据来源:数918年4洛杉奥运会I机AFAAFT/赛与田径赛统手计)册实验求要:1()试求成主,并对分结进行解释果;2(试)用差方积贡累率和S献reec图定确主成分个数;的(3)计算各国第主一成分得分并排的。名11.683实题验一分目析告报(:1>)jscov<.-atmri(c(1x7,177,959,6,5933,2,-,74,-3-,71,914,925,141,311,8127-2,,,14,5,9254,30260,,109,412,,4,41,96,113,106,185102,,244,3,2,53,,81110,,912,017396,,4,56,,2,317,122,44,96,2182,22,8,,7,-2,-4,,4,4234,3,1,33,4,-14,,35,,,23,139,9,-3,34,1,1,6,8,32,393,4),89,9>)js.orc<co-(sjr.ovc)#关相系数矩阵>rond(ujs.oc,2r),1[],2[[,]]3[,4],[5],6][,7][[,8][9],1[,]1.000.840.620.78.60064.-0.84-0.4800-80.[,2]0.841.00.8600.27.901.70-0.897-.707-0.73[,]30.620.681.0004.30.80.7190.68--.608-059.[4,].087.7200.43100.07.0942.-080.0.80--085.[,]50.660.910.78.09710.0.850-.79-007.90-.77[6],0.440.780.910.42.801.05-0.070-07.1-.62[70,-]08.40-.970-68.-.080-0.79-.701.000.980.096[8,]0-80-0.77-..086-0.80-079-.0.17.9801.000.899,][-080-0.73.0.-95-08.5-0.7-0.620796..091.800:答据根相系数关阵矩可看以:出关相较性大,可以做成分主分.析(2)>jsc.vo#协方矩差[,1阵[,]][2,3][4,],[]5[,]6[7,[,]8],9[[1]],77171999563532-74--34[2,]179419425311181217-214[,]95324503206109142
4411[4],9631610518120244235,]5318[110901127369456[,6]322171422961284228[,]7-7-24442343133[8,-4]14352313399[,9-]3411263383948>pac.sj-p<irconm(pcvmotas=.jov)c#成主分析分>pca.j$sloaidngs载#荷Ladonigs:omp.1C1[],03.50[2,0.]65[33,]0.48[34],0.612[5]0,.34[26],02.1[7,][8,7[]9,]Comp2..4860.01510.-5880.491-0.373Comp3.04.2011.83.0263-0.547-04.59-.069-0.22650-.26-0.2628Comp.40.82701.51-.0108-0.82-0.20430438.0.055055.8Comp5.0.3850.-640.84205.6410-.195.1045Comp6.046.4-.105-0.30090.3-0506.96.0130Com.70p.41-60130.0-8.810.41705.61Cmp.o80.263-.2070.01620.3460.674--0.331.2075-0296.Comp9.-0.210.21950-1.85.0522-07.160.843>sumamy(pra.scj)mpIoratncefoocpmoenns:toCpm.1Stnadrdadeviaiton2.93682950roPporitnofOarVanic0e.0962366CumluaitvPeorortponi0.0623966Cop.m5Satdnardedvaitoin90.1442695ProoprtoinfOVaiarcne.056053043CumultaieProvportino0962.15699Cop.2m1.40034420213.988850.7522554Com.p6.615120406.006227446.098768454omC.3p1.1324944010.089211580.84433075omCp.27.416843230.0048329790.199603243Co4m.p10.174004681.0017269107.90616936omC.p82.388282519.000359331709.975975475Copm9S.andarddteivtiano18.57879645ProoptroinfVOarinca0e002.44206Cumu2ltiavPeorporiont10.00000000答由结:果得可三个前本主样成提取的分信息达83.44%(3>)crseepol(trincomp.pap,ypte"=lneis")#Srcee图prncoip.amp3000520500Varaicnes1500000500Comp1.Cop.m3oCmp5.omp.C7Com.9p答由:Sreec图以看可出:坡度从第个三主分成开逐始渐缓变由;(2的结果)得:从可三个主成第后分始累开计贡献率为83.4%480>%结合两者可以得出,主成分的数为3个个。(4答:)第一主成表达分式:y=10.35x100+.536x2+0.834x3+02.1x6+4.023x45+.07126x第二主成分达表:式2y0.=8461+0x1.1x2-5.0885x+0.43194-x.495x0-0.373x566第三主分成表达:式3=0y412x.+0.18132+0.x326x-0.34754-0.x20x5-08.2696x0.25-x67-.0266x80-28.2x9第主成分一映反了低密的度绿蓝滤红片色在卷胶光中应曝较少用;二第主成为对分比分,成映了反中密度的红色滤色片对高色彩卷胶曝光响较影,运大在风景用照建或摄筑影,可时以作出制风“暴”似的云彩,光及树木月,造出超写实的特营效这,种色片是绝滤必对要的,中而高密的蓝色度色滤强调山片中的谷岚及穿山透水与面烟雾光的等方式,线达来凸显到天层次空的效果对彩色胶卷;三第主分成综为合分,成反了映同样是绿蓝红种三颜色滤色,只片在高密有度才下对卷曝胶起到光决作定用,密越度高,明对说绿红蓝种三色吸收颜越,多而影从响照到片的彩。色7实验题目分析报告:二所有的格表经过xcEel表格处理及留保小数三点。位(
1)>ATFre=d.cas("v:DFA.Tcsv,"eha=d)T导入数据#>ATF=FT[-A1,]去#第一行缺掉失值>AFT=A1F[,T-]1#去掉一列名第>字Simgaap.va=rAF(T1#算出)数据方差>的Pa.=por(cFAT1)#算出据的数方差协P.>paX10m0.s..2X0m.0s..X40m0.s..X008m..mi.nX150m0.mi.n.X5000m..im.X1n.000m0..min.Maarhtn..mios.nX001..s.mX02m0..s.X04m..0sX800m...mi.n1.00009.32.0438.05780.9231.0000.5380.0880.8340.5381.0000870..7058.8080.87001.000.00710.770.58360.19906.12069.70.778086.30.3460.9860.7780.69805.120.597.7050.08071X050.m.im.n07.010775.0.863.91091.000.09290.9530.668X005m.0.in.m0.216.069077780.863..0929100.009.57.039X2010.0m..mi0n.0.364.0698.0877.0680.99305975.1.0000.493Mraahton.mi.sn0.5210..597.7005.0807.0866.0320.9431.9000>priconmpap.=ripncom(ApTF1cor,T=)用相关系#数主成做分析分>sumarm(pyrincmop.pa#)总统计汇Com.p1omC.p2Com.3pComp4.omCp.5Comp.6Com.7Comp.p8Standarddveaitoin2574.0.396039.703520..82032.61.014201.50ropoPtrinoofVaiarcne0.8820.1010.200.0050.10010.0090.0060.003uCulatimvePorprtiono.82080.938.907509.3709.30.99180.9719.000>ripncop.mp$aloadngis#主分分成析X100..sm.X020.ms..4X00..ms.8X00..mmn.Xi5100..mim.nX0500.m.imn.X1.000m0.m.i.nMaatrhno.min..sCmop10..18303.370.5630.6930.370.3346.36700.423omp.C205.66.4600.512-.54-0.0141-.310-2030.7-.044omC.p03.360.335-70.5440.-74-017.01.9510185.0265.omC.4p0.120-2.609.6560.359-.39050.2-55-.1033029.6omp.50C2.3-60156.-0217.-0.422-0.874014.022.90.1499Comp.60.953-0.56-60.561.1580-0.592-.1027-03.16Comp.07.1390-1.60.61001.4540-.9770.93C7mpo8..00170.1420.16SlSoaingsdrPooptirnVoaCruulamtivVeraComp.110.125.0215Copm21.0125..025Cmpo3.1.0215.0375Cmop4.1.125005.Comp5.10.1250625.Comp6.10.125.750Comp7.1.12050.857Comp.81.12051oCpm1.为综合力Comp能.2为限极力能即(间爆发力和瞬力)耐82()su>mmrya(prncimop.ap)汇#总计统omC.1pomCp2Comp..3oCmp.Co4p.5moCmp.C6mo.7pComp.8Satnardddeiavtino257.40.9603397.0.5320.283.0612.01420.51Pr0poortoinoVfriaance0.2880.100.12000.100.01050.090.0006.0030CmulutaievrPpoortoni0.288.93080.9570.973.90830991..997010.00后最列为方差一计累献贡,可率取Comp以.1和Cmp.o2为主分成。>csreelptopr(nicmp.oa,type="piles")n画出了scr#e图e根据可图以得,取出omCp.1和Copm2.为成分9主(3)y=>pricnompap.$cores[s,]#计1主算成得分分>y在Excel中理并整序:排CuontryoPitNetnhrealndCoskoIlsnad1s.0536BraizleWtersnSaoam7.298NeweaZlanMdauriitus.2984weSdnePapuNaeGuwniae.3945wStierlzanSindapgor3e1.15FilnndaTahlaind27.78CaadnaInondesai.7732PolndGaateuamla.2976eBliumgoCstaRcia231.8KenyaPhliippnei2.s900
FrncaeBrma1.990uAsurtlaaiDomiicnnaRpebuicl.7131eGmar(n.R.F)alMayias1.247Gemrna(DR..)D.P.RoKea1r6.99USSRaipTie.9509ItalyBrmedau0.746GreatBri.&Nt.IerandlsIrela0.394SAUChian04.13Trueky.0296Argentian.2604uxemboLurg0.223orKea.0029Idian0-1.6Gree7ec-0.338Chile-0.358Colubmai-0.394exicMo-0685.usAtrai-.8051orNway0.-819Irelnad0-8.92Portgaul-092.5eDnamkr1.12-3uRanmia-.2180Hnuargy1.-12Ja6pan-.1249Cezchsloovakai1-385.Spain-1.449-1.507-1.75-31.61-4161.-8.6145-1.708-1.762-2019-.2.006-2.881-2.1922-4.96-25.67-2.6142-6.15-.72253-.520-3.42601