应用统计分析第三次作业.docxVIP

应用统计分析第三次作业 数学与应用数学一班 陈正康 2012101010019 【摘要】在一些实际问题中，影响结果的变量往往有很多个，本文利用主成分分析法分析2005年至2011年我国不变价国内生产总值的构成对影响我国人均生产总值的因素做主成分分析。由此得出我国国民生产总值的相关经济数据。 【关键词】主成分分析法 不变价国内生产总值的构成 问题描述 （1）对我国不变价国内生产总值的各个构成进行主成分分析； （2）确定影响我国国内人均生产总值的主要影响因素； 数据如下:数据来源于国家统计局网站(/tjsj/ndsj/2013/indexch.htm) 二、问题分析及解决 相关变量的说明 变量变量说明 Y国内生产总值 X1 第一产业产值 X2 工业产值 X3 建筑业产值 X4 第三产业产值 问题的求解 2.1数据的输入 conomy=read.table(3.txt,head=TRUE) 2.2 作主成分分析 conomy.pr-princomp(~x1+x2+x3+x4,data=conomy,cor=T) summary(conomy.pr,loadings=TRUE) Importance of components: Comp.1 Comp.2 Comp.3 Comp.4 Standard deviation 1.998366 0.065727220 0.0470316611 3.515471e-07 Proportion of Variance 0.998367 0.001080017 0.0005529943 3.089633e-14 Cumulative Proportion 0.998367 0.999447006 1.0000000000 1.000000e+00 Loadings: Comp.1 Comp.2 Comp.3 Comp.4 x1 -0.500 0.381 0.160 0.761 x2 -0.500 0.465 0.357 -0.637 x3 -0.500 -0.856 -0.124 x4 -0.500 -0.797 0.339 前两个主成分已达到99%的贡献率。 2.3 主成分回归分析 pre-predict(conomy.pr) conomy$z1-pre[,1];conomy$z2-pre[,2] lm.sol-lm(y~z1+z2,data=conomy) summary(lm.sol) Call: lm(formula = y ~ z1 + z2, data = conomy) Residuals: 1 2 3 4 5 6 7 8 -52.91 -156.06 227.05 1684.85 -1694.38 -132.11 -690.37 813.93 Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 291479.6 416.7 699.529 1.13e-13 *** z1 -42014.8 208.5 -201.500 5.71e-11 *** z2 -45363.7 6339.5 -7.156 0.000828 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 1179 on 5 degrees of freedom Multiple R-squared: 0.9999, Adjusted R-squared: 0.9998 F-statistic: 2.033e+04 on 2 and 5 DF, p-value: 1.677e-10 2.4线性回归分析 lm.sol-lm(y~x1+x2+x3+x4,data=conomy) summary(lm.sol) Call: lm(formula = y ~ x1 + x2