Auto数据集多元线性回归

import pandas as pd
import statsmodels.formula.api as smf
auto = pd.read_csv('https://liangdaima.com/static/data/statistics/Auto.csv')
auto['horsepower'] = pd.to_numeric(auto['horsepower'], errors='coerce')
auto = auto.dropna()
model = smf.ols('mpg ~ horsepower + weight + acceleration', data=auto).fit()
print('R-squared:', model.rsquared)
print('回归系数:')
print(model.params)