毛宇豪_13062304_上机实习作业2..docx

2.1、研究货运总量y（万吨）与工业总产值x1（亿元）、农业总产值x2（亿元）、居民非商品支出x3（亿元）的关系。数据见表2.1表2.1编号货运总量y(万吨)工业总产值x1（亿元）农业总产值x2（亿元）居民非商品支出x3（亿元）116070351.0226075402.4321065402.0426574423.0524072381.2622068451.5727578424.0816066362.0927570443.21025065423.0计算出y，x1，x2，x3的相关系数矩阵。求出y与x1，x2，x3的三元回归方程。对所求的方程作方程的显著性检验。对每一个回归系数作显著性检验如果有的回归系数没有通过显著性检验，将其剔除，重新建立回归方程，并作回归方程的显著性检验和回归系数的显著性检验。求出每一个回归系数的置信水平为95%置信区间求标准化回归方程求当x01=75，x02=42，x03=3.1时的y的预测值，给定置信水平95%，软件计算精确置信区间。作异常值检验结合回归方程对问题做一些简单分析。解：R程序：X-data.frame(y=c(160,260,210,265,240,220,275,160,275, 250),x1=c(70,75,65,74,72,68,78,66,70,65),x2=c(35,40,40,42,38,45,42,36,44,42),x3=c(1.0,2.4,2.0,3.0,1.2,1.5,4.0,2.0,3.2,3.0))cor(X)输出结果： y x1 x2 x3y 1.0000000 0.5556527 0.7306199 0.7235354x1 0.5556527 1.0000000 0.1129513 0.3983870x2 0.7306199 0.1129513 1.0000000 0.5474739x3 0.7235354 0.3983870 0.5474739 1.0000000由以上输出结果可知：第一题相关矩阵为 R程序：lm.sol-lm(y~x1+x2+x3,data=X)summary(lm.sol)输出结果：Call:lm(formula = y ~ x1 + x2 + x3, data = X)Residuals: Min 1Q Median 3Q Max -25.198 -17.035 2.627 11.677 33.225 Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) -348.280 176.459 -1.974 0.0959 .x1 3.754 1.933 1.942 0.1002 x2 7.101 2.880 2.465 0.0488 *x3 12.447 10.569 1.178 0.2835 ---Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1Residual standard error: 23.44 on 6 degrees of freedomMultiple R-squared: 0.8055, Adjusted R-squared: 0.7083 F-statistic: 8.283 on 3 and 6 DF, p-value: 0.01487由以上结果知：回归方程为：根据模型摘要，F统计量服从自由度为(3,6)的F分布，给定显著性水平，查表得，由模型摘要得F值=8.2834.76，p值=0.015，拒绝原假设。由模型摘要得到F=8.283,P=0.0150.05，说明在置信水平95%下，回归模型显著。做t检验:设原假设为，统计量服从自由度为n-p-1=6的t分布，给定显著性水平0.05，查得单侧检验临界值为1.943，由模型摘要，的t值=1.9421.943，处在否定域边缘。的t值=2.4651.943。拒绝原假设。根据模型摘要，在显著性水平时，只有的P值0.05，通过检验，即只有的回归系数为显著；其余的资料变量的P值均大于0.05，即，的系数均不显著。R程序：lm.new-lm(y~x1+x2,data=X)summary(lm.new)输出结果：Call:lm(formula = y ~ x1 + x2, data = X)Residuals: Min 1Q Median 3