R 语言中的多元线性回归

示例

sessionInfo()  #  查询版本及系统和库等信息

# 工作目录设置

getwd()

path <- "E:/RSpace"

setwd(path)

rm(list=ls()) # 清空内存中的变量

# state.x77 # 展示基础安装中的 women 数据集

str(state.x77) # 查看 state.x77 的数据结构

# 查看数据集 state.x77 的相关属性

attributes(state.x77)

typeof(state.x77) # 数据类型

class(state.x77) # 数据结构

dim(state.x77) # 维度

# 从矩阵中抽取数，创建数据框

states<- as.data.frame(state.x77[,c("Murder", "Population","Illiteracy", "Income", "Frost")])

summary(states) # 查看 states 的摘要统计量

# 检测二变量关系

cor(states) # 相关系数矩阵

# 添加包

install.packages("car")

library(car)

# 创建散点图矩阵

scatterplotMatrix(states, spread=FALSE, smoother.arg=list(lty=2), main="散点图矩阵")

# 多元线性回归

# 研究州犯罪率跟人口、文盲率、平均收入和结霜期的关系。

fit <- lm(Murder ~ Population + Illiteracy + Income + Frost, data=states)

fit # 查看模型的结果对象列表

summary(fit) # 展示模拟模型的详细结果