R 《回归分析与线性统计模型》page164 单变量、多变量多项式模型

rm(list = ls())
library(openxlsx)
library(leaps)

#单变量多项式模型#
data = read.xlsx("table7-1.xlsx")
head(data)
plot(data$f,data$mw)

attach(data)
new_data = as.data.frame(cbind(mw,f,f**2,f**3))
names(new_data) = c("mw","f","f2","f3")
fm_3 = lm(mw~.,new_data) #拟合一个三次多项式模型
summary(fm_3) #用电负荷峰值三次多项式模型的OLS估计结果
anova(fm_3) #用电负荷峰值三次多项式模型的方差分析表
detach(data)

new_data = new_data[,-4] 
fm_2 = lm(mw~.,new_data) #拟合一个二次多项式模型
summary(fm_2) #用电负荷峰值二次多项式模型的OLS估计结果
anova(fm_2) #用电负荷峰值二次多项式模型的方差分析表


#多变量多项式模型#
rm(list = ls())
data = read.xlsx("table7-2.xlsx")
head(data)
attach(data)
new_data =as.data.frame(cbind(y,x1,x2,x1*x2,x1**2,x2**2))
fm = lm(y~.,new_data) #二次多项式模型
summary(fm) #温度、压强及成品质量数据OLS估计结果
anova(fm) #方差分析