北京工业大学数学建模工程数学实验6.docx

下载文档 降价啦

1
0
约1.77万字
约 36页
2017-10-27 发布于重庆
举报
版权申诉
保障服务

北京工业大学数学建模工程数学实验6.docx

1、本文档共36页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

北京工业大学数学建模工程数学实验6

(1)输入：?import-data.frame(?X1=c(149.3,161.2,171.5,175.5,180.8,190.7,202.1,212.4,226.1,231.9,239.0),?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.0,19.1,18.8,20.4,22.7,26.5,28.1,27.6,26.3)?)?lm.sol-lm(Y?~?X1+X2+X3,?data=import)?summary(lm.sol)?结果为：?Call:lm(formula = Y ~ X1 + X2 + X3)Residuals: Min 1Q Median 3Q Max -0.52367 -0.38953 0.05424 0.22644 0.78313 Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) -10.12799 1.21216 -8.355 6.9e-05 ***X1 -0.05140 0.07028 -0.731 0.488344 X2 0.58695 0.09462 6.203 0.000444 ***X3 0.28685 0.10221 2.807 0.026277 * ---Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.4889 on 7 degrees of freedomMultiple R-squared: 0.9919, Adjusted R-squared: 0.9884 F-statistic: 285.6 on 3 and 7 DF, p-value: 1.112e-07得出Y关于X1、X2、X3的线性回归方程：Y=-10.12799-0.05140X1+0.58695 X2+0.28685X3通过显著性检验得：?X1不显著，X2极为显著，X3显著，Y极为显著。(2)由于1X不显著，所以该回归方程不合理。输入：?lm.step-step(lm.sol)?运行得到：?Start:??AIC=-12.72?Y?~?X1?+?X2?+?X3?Df?Sum?of?Sq??????RSS??????AIC?-?X1????1????0.1278???1.8008?-13.9065?none????????????????1.6729?-12.7164?-?X3????1????1.8824???3.5554??-6.4238?-?X2????1????9.1967??10.8697???5.8689??Step:??AIC=-13.91?Y?~?X2?+?X3?????????Df?Sum?of?Sq?????RSS?????AIC?none????????????????1.801?-13.907-?X2????1?????9.651??11.452???4.443?-?X3????1???191.668?193.469??35.539输入：summary(lm.step)运行得到：Call:?lm(formula?=?Y?~?X2?+?X3,?data?=?import)??Residuals:??????Min???????1Q???Median???????3Q??????Max??-0.67365?-0.28064??0.05352??0.19187??0.73905?Coefficients:??????????????Estimate?Std.?Error?t?value?Pr(|t|)?????(Intercept)?-9.742740???1.059489??-9.196?1.58e-05?***?X2???????????0.596052???0.091028???6.548?0.000179?***?X3???????????0.212305???0.007276??29.180?2.06e-09?***?---?Signif.?codes:??0?‘***’?0.001?‘**’?0.01?‘*’?0.05?‘.’?0.1?‘?’?1?