R 语言中的多项式回归

代码

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

# 工作目录设置
getwd()
path <- "E:/RSpace/R_in_Action"
setwd(path)

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

women # 展示基础安装中的 women 数据集
str(women) # 查看 women 的数据结构
summary(women) # 查看 women 的摘要统计量

# 通过添加高次项来提高回归模型的精度
fit2 <- lm(weight ~ height + I(height^2), data=women)
fit2 # 查看模型的结果对象列表

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

women$height
fitted(fit2) # 列出拟合模型的预测值
residuals(fit2) # # 列出拟合模型的残差值

# 绘制身高体重的散点图
plot(women$height, women$weight,
     xlab="Height (in inches)",
     ylab="Weight (in pounds)")

lines(women$height, fitted(fit2)) # 添加回归线