regnam r语言_R语言一元线性回归

2020-06-25

例子来源：数学建模的三十二种常规方法

exam1：合金的强度 y 与其中的碳含量 x 有比较密切的关系，从生产中收集了一批数据如下表 。

x 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18

y 42.0 41.5 45.0 45.5 45.0 47.5 49.0 55.0 50.0

试先拟合一个函数 y(x) ，再用回归分析对它进行检验。

求解：

(1)先画出散点图

x

y

plot(x, y)

(2)设回归模型为y = β0 + β1*x

lm.reg

(3)结果分析

summary(lm.reg)

confint(lm.reg, level=0.95)

plot(lm.reg)

结果如下：

Call:

lm(formula= y ~ 1 + x)

Residuals:

Min     1Q  Median      3Q    Max

-2.2222-1.0972 -0.4722  0.7778  4.1528

Coefficients:

Estimate Std. Error t value Pr(>|t|)

(Intercept)   27.472     3.716   7.393  0.00015 ***

x               137.500    26.103  5.268  0.00116 **

---

Signif.codes:

0 ‘***’0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 2.022 on 7 degrees of freedom

Multiple R-squared:  0.7985,   Adjusted R-squared:  0.7698

F-statistic:27.75 on 1 and 7 DF,  p-value: 0.001164

> confint(lm.reg, level=0.95)

2.5%        97.5 %

(Intercept) 18.68506    36.25939

x              75.77550    199.22450

说明：

β0 = 27.4722 ， β1 = 137.5000

β0 的置信区间是[18.6851,36.2594]， β1 的置信区间是[75.7755,199.2245]；

概率p-value: 0.001164，模型成立