数学模型实验8.docx

数学建模实验八 多元分析实验1 回归分析解：(1)R程序如下：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)summary(lm.sol, correlation = T, symbolic.cor = F)得到的结果如下：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.293x 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 freedomMultiple R-squared: 0.9781, Adjusted R-squared: 0.9754 F-statistic: 357.5 on 1 and 8 DF, p-value: 6.33e-08 由计算结果知：系数不显著，极为显著。相关系数和方程通过检验。得到的一元线性回归模型为：(2)R程序如下：new - data.frame(x = 7)predict(lm.sol, new, interval=prediction, level=0.95)得到的结果如下： fit lwr upr1 2690.227 2454.971 2925.484预测值为2690.227,区间为?(2454.971,2925.484)(3)输入程序如下：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)new - data.frame(x = seq(3.5, 8.8, by=0.1))pp-predict(lm.sol, new, interval=prediction)pc-predict(lm.sol, new, interval=confidence)matplot(new$x, cbind(pp, pc[,-1]), type=l, xlab=X, ylab=Y, lty=c(1,5,5,2,2), col=c(blue, red, red, brown, brown), lwd=2)points(x,y, cex=1.4, pch=21, col=red, bg=orange)legend(4.1, 3000, c(Points, Fitted, Prediction, Confidence), pch=c(19, NA, NA, NA), lty=c(NA, 1,5,2), col=c(orange, blue, red, brown))savePlot(predict, type=eps)得到的结果如下：分析：从图中看出所有的点都落在预测区间内。只有第7个点不在置信区间内，其余都在置信区间内。可以看出拟合情况较好。2 回归分析解：输入程序如下：intellect-data.frame( 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, data=intellect)summary(lm.sol)influence.measures(lm.sol)op - par(mfrow=c(2,2), mar=0.4+c(4,4,1,1), oma= c(0,0,2,0))plot(lm.sol, 1:4)par(op)savePlot(diagnoses1, type=eps)运行结果如下：Call:lm(formula = y ~ 1 + x, data = intellect)Residuals: Min