R version (2017-11-30) -- \
Copyright (C) 2017 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit)
R是自由软件,不带任何担保。 在某些条件下你可以将其自由散布。
用'license()'或'licence()'来看散布的详细条件。
R是个合作计划,有许多人为之做出了贡献. 用'contributors()'来看合作者的详细情况
用'citation()'会告诉你如何在出版物中正确地引用R或R程序包。
用'demo()'来看一些示范程序,用'help()'来阅读在线帮助文件,或 用'()'通过HTML浏览器来看帮助文件。 用'q()'退出R.
[原来保存的工作空间已还原]
> h=(\
Error in (file = file, header = header, sep = sep, quote = quote, :
找不到对象'true' > h=(\ > h
地区 x1 x2 x3 x4 x5 x6 x7 x8 x9 y 1 北京 7535 2639 1971 1658 3696 84742 87475 240462 天津 7344 1881 1854 1556 2254 61514 93173 200243 河北 4211 1542 1502 1047 1204 38658 36584 125314 山西 3856 1529 1439 906 1506 44236 33628 122125 内蒙古 5463 2730 1584 1354 1972 46557 63886 177176 辽宁 5809 2042 1433 1310 1844 41858 56649 165947 吉林 4635 2045 1594 1448 1643 38407 43415 146148 黑龙江 4687 1807 1337 1181 1217 36406 35711 129849 上海 9656 2111 1790 1017 3724 78673 85373 2625310 江苏 6658 1916 1437 1058 3078 50639 68347 1882511 浙江 7552 2110 1552 1228 2997 50197 63374 2154512 安徽 5815 1541 1397 1143 1933 44601 28792 1501213 福建 7317 1634 1754 773 2105 44525 52763 1859314 江西 5072 1477 1174 671 1487 38512 28800 1277615 山东 5201 2197 1572 1005 1656 41904 51768 15778
16 河南 4607 1886 1191 1085 1525 37338 31499 13733 17 湖北 5838 1783 1371 1030 1652 39846 38572 14496 18 湖南 5442 1625 1302 918 1738 38971 33480 14609 19 广东 8258 1521 2100 1048 2954 50278 54095 22396 20 广西 5553 1146 1377 884 1626 36386 27952 1424421 海南 6556 865 1521 993 1320 39485 32377 1445722 重庆 6870 2229 1177 1102 1471 44498 38914 1657323 四川 6074 1651 1284 773 1587 42339 29608 1505024 贵州 4993 1399 1014 655 1396 41156 19710 1258625 云南 5468 1760 974 939 1434 37629 22195 1388426 西藏 5518 1362 845 467 550 51705 22936 1118427 陕西 5551 1789 1322 1212 2079 43073 38564 1533328 甘肃 4602 1631 1288 1050 1388 37679 21978 1284729 青海 4667 1512 1232 906 1097 46483 33181 1234630 宁夏 4769 1876 1193 1063 1516 47436 36394 1406731 新疆 5239 2031 1167 1028 1281 44576 33796 13892> lm=lm(y~x1+x2+x3+x4+x5+x6+x7+x8+x9,data=h) > lm Call:
lm(formula = y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9, data = h)
Coefficients:
(Intercept) x1 x2 x3 x4 x5 x6 x7 x8 x9
> summary(lm) Call:
lm(formula = y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9, data = h)
Residuals:
Min 1Q Median 3Q Max
Coefficients:
Estimate Std. Error t value Pr(>|t|) (Intercept) +02 +03
x1 +00 *** x2 +00 *** x3 +00 *** x4 x5 +00 *** x6 x7 x8 +01 +01 x9 +01 +02 ---
Signif. codes: 0 ‘***’ ‘**’ ‘*’ ‘.’ ‘ ’
Residual standard error: on 21 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: on 9 and 21 DF, p-value: <
> pre=(lm)
> res=residuals(lm) > sd(res) [1]
1