R语言建模lesson3

载入数据

library(MASS)
library(ISLR2)
#install.packages("ISLR2")
head(Boston)

简单线性回归

head(Boston)# 显示前五行数据
lm.fit <- lm(medv ~ lstat, data = Boston)
attach(Boston)
lm.fit <- lm(medv ~ lstat)
lm.fit
summary(lm.fit)

names(lm.fit)# 查看变量名称
coef(lm.fit)# 获得这个模型每个自变量所对应的回归系数
confint(lm.fit)# 为了获得系数估计的置信区间

#该函数可用于生成置信区间和预测区间，用于预测给定值 的
predict(lm.fit, data.frame(lstat = (c(5, 10, 15))),
        interval = "confidence")
predict(lm.fit, data.frame(lstat = (c(5, 10, 15))),
    interval = "prediction")
    
plot(lstat, medv)
abline(lm.fit)

绘制带斜率和截距的线

plot(lstat, medv)
abline(lm.fit, lwd = 3)
abline(lm.fit, lwd = 3, col = "red")

点图参数

plot(lstat, medv, col = "red")

在 R 语言中，pch 是绘制散点图时用来指定点的形状的参数。默认情况下，pch 的取值为 1，表示绘制普通的圆点。您可以使用不同的整数值来指定不同的点形状，以下是一些常用的 pch 参数取值及对应的点形状：

pch = 0：方块
pch = 1：圆点（默认值）
pch = 2：三角形点
pch = 3：十字叉
pch = 4：X 叉
pch = 5：菱形点
pch = 6：六边形点
pch = 7：十字十
pch = 8：星号
pch = 9：五角星
pch = 10：十字加号

plot(lstat, medv, pch = 20)
plot(lstat, medv, pch = "+")
plot(1:20, 1:20, pch = 1:20)

分屏

par(mfrow = c(2, 2))
plot(lm.fit)

作业：使用Auto数据集进行简单线性回归

attach(Auto)
View(Auto)
head(auto)# 显示前五行数据
lm.fit <- lm(mpg ~ horsepower, data = Auto)
lm.fit <- lm(mpg ~ horsepower)
lm.fit
summary(lm.fit)