python线性回归

import numpy as np
import pandas as pd
from sklearn.datasets import load_boston
import statsmodels.formula.api as smf

filename = r"D:\RUI\MathModeling\2021B乙醇\装料方式1数据.xlsx"
data = pd.read_excel(filename)
print(data)
print(data.columns)
mod = smf.ols(formula='t~x0+x1+x2+x3+x4+x0x0+x0x1+x0x2+x0x3+x0x4+x1x1+x1x2+x1x3+x1x4+x2x2+x2x3+x2x4+x3x3+x3x4+x4x4', data=data)
res = mod.fit()
print(res.summary())

mod2 = smf.ols(formula='t~x0+x1+x0x0+x0x1+x0x4+x1x4+x3x3+x3x4', data=data)
res2 = mod2.fit()
print(res2.summary())