资源描述:
《北京工业大学薛毅老师工程数据建模实验6多元分析》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、6多元分析实验学号:G201212xxx姓名:xxxemail:xxx@163.com6.1回归分析法国经济数据分析,考虑进口总额Y与三个自变量:国内总产值X1,存储量X2,总消费量X3(单位为10亿法郎)之间的关系。现收集了1949年至1959年共11年的数据,如表6.1所示,试对此数据进行分析。表6.1法国经济分析数据序号X1X2X3Y1149.34.2108.115.92161.24.1114.816.43171.53.1123.219.04175.53.1126.919.15180.81.1132
2、.118.86190.72.2137.720.47202.12.1146.022.78212.45.6154.126.59226.15.0162.328.110231.95.1164.327.611239.00.7167.626.3(1)求出Y关于X1,X2和X3的线性回归方程,并对方程作显著性检验;(2)分析所得到的回归方程是否合理,对变量作逐步回归;(3)假设某年的国内总产值(X1)、存储量(X2)和总消费量(X3)分别为240、4.5和170(单位为10亿法郎),给出该年进口总额(Y)的预测值、相应
3、的置信区间和预测区间(α=0.05)。解:(1)输入程序:X1<-c(149.3,161.2,171.5,175.5,180.8,190.7,202.1,212.4,226.1,231.9,239)X2<-c(4.2,4.1,3.1,3.1,1.1,2.2,2.1,5.6,5.0,5.1,0.7)X3<-c(108.1,114.8,123.2,126.9,132.1,137.7,146.0,154.1,162.3,164.3,167.6)Y<-c(15.9,16.4,19,19.1,18.8,20.4,2
4、2.7,26.5,28.1,27.6,26.3)lm.sol<-lm(Y~X1+X2+X3)summary(lm.sol)运行结果:Call:lm(formula=Y~X1+X2+X3)Residuals:Min1QMedian3QMax-0.52367-0.389530.054240.226440.78313Coefficients:EstimateStd.ErrortvaluePr(>
5、t
6、)(Intercept)-10.127991.21216-8.3556.9e-05***X1-0.051400.
7、07028-0.7310.488344X20.586950.094626.2030.000444***X30.286850.102212.8070.026277*---Signif.codes:0‘***’0.001‘**’0.01‘*’0.05‘.’0.1‘’1Residualstandarderror:0.4889on7degreesoffreedomMultipleR-squared:0.9919,AdjustedR-squared:0.9884F-statistic:285.6on3and7DF,
8、p-value:1.112e-07结果分析:可知Y关于的线性回归方程:y=-10.12799-0.05140x1+0.58695x2+0.28685x3其中:x1无标记说明不显著,x2标记“***”说明极为显著,x3标记“*”显著。(2)由(1)中的数据可以得知新的分析函数anova(lm.sol)R程序输入程序:X1<-c(149.3,161.2,171.5,175.5,180.8,190.7,202.1,212.4,226.1,231.9,239)X2<-c(4.2,4.1,3.1,3.1,1.1,2
9、.2,2.1,5.6,5.0,5.1,0.7)X3<-c(108.1,114.8,123.2,126.9,132.1,137.7,146.0,154.1,162.3,164.3,167.6)Y<-c(15.9,16.4,19,19.1,18.8,20.4,22.7,26.5,28.1,27.6,26.3)lm.sol<-lm(Y~X1+X2+X3,data=blood)summary(lm.sol)anova(lm.sol)运行结果:Call:lm(formula=Y~X1+X2+X3)Residuals
10、:Min1QMedian3QMax-0.52367-0.389530.054240.226440.78313Coefficients:EstimateStd.ErrortvaluePr(>
11、t
12、)(Intercept)-10.127991.21216-8.3556.9e-05***X1-0.051400.07028-0.7310.488344X20.586950.094626.2030.000444***X30.286850.