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) # 展示模拟模型的详细结果

猜你喜欢

转载自www.cnblogs.com/shanger/p/12090996.html