Auto数据集回归诊断图

import pandas as pd
import seaborn as sns
import matplotlib.pyplot as plt
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()
fig, axes = plt.subplots(1, 2, figsize=(14, 5))
sns.regplot(x='horsepower', y='mpg', data=auto, ax=axes[0])
axes[0].set_title('mpg vs horsepower 回归图')
sns.residplot(x='horsepower', y='mpg', data=auto, ax=axes[1])
axes[1].set_title('残差图')
plt.show()