应用统计分析第三次作业

《应用统计分析第三次作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、应用统计分析第三次作业数学与应用数学一班陈正康2012101010019【摘要】在一些实际问题中，影响结果的变量往往有很多个，本文利用主成分分析法分析2005年至2011年我国不变价国内生产总值的构成对影响我国人均生产总值的因素做主成分分析。由此得出我国国民生产总值的相关经济数据。【关键词】主成分分析法不变价国内生产总值的构成一、问题描述(1)对我国不变价国内生产总值的各个构成进行主成分分析；(2)确定影响我国国内人均生产总值的主要影响因素；数据如下：数据来源于国家统计局网站(http://www.stats・gov.cn/tjsj/nds

2、j/2013/indexch・htm)1不变价国内生产总值草位，亿元国内生产1E份第一产业第二产业»—IXJb工业建筑业気一m按2005年价格计算18493742242008759817723081036737491932083810235410993285871751121534855116200623789282442241142906100170114120599179720072608129257359125579.71101174154622109497420082848448268126138062711973141833131

3、199694200931460252795781549758134175520800313166892010按2010年价格计尊4015128405336187383216072222666101735960201043885304225652066551177403629251518994142011二、问题分析及解决相关变量的说明变量变量说明Y国内生产总值XI第一产业产值X2工业产值X3建筑业产值X4第三产业产值问题的求解2.1数据的输入conomy=read.table(，z3・txt〃,head=TRUE)2.2作主成分分析>con

4、omy.pr<-princomp("xl+x2+x3+x4,data二conomy,cor=T)>summary(conomy.pr,loadings二TRUE)Importanceofcomponents:Comp.4Comp.1Comp.2Comp.3Standard3.515471e-07deviation1.9983660.0657272200.0470316611Proportion3.089633e-14ofVariance0.9983670.0010800170.0005529943Cumulative1.000000e+00

5、Proportion0.9983670.9994470061.0000000000Loadings:Comp.1Comp.2xl-0.5000.381x2-0.5000.465x3-0.500x4-0.500-0.797Comp.3Comp・40.1600.7610.357-0.637-0.856-0.1240.339＞前两个主成分已达到99%的贡献率。入=0.入=02.3主成分回归分析>pre〈-predict(conomy.pr)>conomy$zl<-pre[,1];conomy$z2<-pre[,2]>1m.sol<~1m(y^z1

6、+z2,data=conomy)>summary(lm.sol)Call:Im(formula二y~zl+z2,data二conomy)EstimateStd.(Intercept)291479.6zl-42014.8Error416.7tvaluePr(>111)699.5291.13e-13***20&5-201.5005.71e-ll***z2-45363.76339.5-7.1560.000828***Signif.codes:0'***，10.001，0.010.050.1Residuals:12345678-52.91-156.

7、06227.051684.85-1694.38-132.11-690.37813.93Coefficients:Residualstandeirderror:1179on5degreesoffreedomMultipleR-squared:0.9999,AdjustedR-squared:0.9998F-statistic:2.033e+04on2and5DF,p-value:1.677e~10y=291479.6-42014.8z；—45363.7z；2・4线性回归分析>Im.sol〈-Im(y~xl+x2+x3+x4,data=cono

8、my)>summary(lm.sol)Call:lm(formula二y~xl+x2+x3+x4,data=conomy)Residuals:12345785.135e+022.