残差分析

def solve():
    from sklearn.datasets import load_iris
    import matplotlib.pyplot as plt
    import numpy as np
    iris = load_iris()
    x = iris.data[:, ____]
    y = iris.data[:, ____]
    x_mean, y_mean = np.mean(____), np.mean(____)
    slope = np.sum((____) * (____)) / np.sum((____) ** ____)
    intercept = y_mean - slope * x_mean
    y_pred = slope * x + intercept
    residuals = y - y_pred
    plt.figure(figsize=(____, ____))
    plt.scatter(____)
    plt.axhline(____, color='____', linestyle='____', linewidth=____)
    plt.xlabel(____)
    plt.ylabel(____)
    plt.title(____)
    plt.grid(____)
    plt.show()