统计建模与R软件第六章习题.doc

第六章 1． a． x - c(5.1, 3.5, 7.1, 6.2, 8.8, 7.8, 4.5, 5.6, 8.0, 6.4) y - c(1907, 1287, 2700, 2373, 3260, 3000, 1947, 2273, 3113,2493) plot(x,y) X与Y线性相关 b． x - c(5.1, 3.5, 7.1, 6.2, 8.8, 7.8, 4.5, 5.6, 8.0, 6.4) y - c(1907, 1287, 2700, 2373, 3260, 3000, 1947, 2273, 3113,2493) lm.sol-lm(y~1+x) summary(lm.sol) Call: lm(formula = y ~ 1 + x) Residuals: Min 1Q Median 3Q Max -128.591 -70.978 -3.727 49.263 167.228 Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 140.95 125.11 1.127 0.293 x 364.18 19.26 18.908 6.33e-08 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 96.42 on 8 degrees of freedom Multiple R-Squared: 0.9781, Adjusted R-squared: 0.9754 F-statistic: 357.5 on 1 and 8 DF, p-value: 6.33e-08 回归方程为Y=140.95+364.18X，极为显著 d． new - data.frame(x=7) lm.pred - predict(lm.sol,new,interval=prediction,level=0.95) lm.pred fit lwr upr [1,] 2690.227 2454.971 2925.484 Y(7)= 2690.227, [2454.971,2925.484] 2． out-data.frame( + x1 - c(0.4,0.4,3.1,0.6,4.7,1.7,9.4,10.1,11.6,12.6,10.9,23.1,23.1,21.6,23.1,1.9,26.8,29.9), + x2 - c(52,34,19,34,24,65,44,31,29,58,37,46,50,44,56,36,58,51), + x3 - c(158,163,37,157,59,123,46,117,173,112,111,114,134,73,168,143,202,124), + y - c(64,60,71,61,54,77,81,93,93,51,76,96,77,93,95,54,168,99) + ) lm.sol-lm(y~x1+x2+x3,data=out) summary(lm.sol) Call: lm(formula = y ~ x1 + x2 + x3, data = out) Residuals: Min 1Q Median 3Q Max -27.575 -11.160 -2.799 11.574 48.808 Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 44.9290 18.3408 2.450 0.02806 * x1 1.8033 0.5290 3.409 0.00424 ** x2 -0.1337 0.4440 -0.301 0.76771 x3 0.1668 0.1141 1.462 0.16573 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 19.93 on 14 degrees of freedom