ANOVA方差分析

def solve():
    from sklearn.datasets import load_iris
    import matplotlib.pyplot as plt
    import numpy as np
    from scipy import stats
    iris = load_iris()
    setosa = iris.data[iris.target == ____, ____]
    versicolor = iris.data[iris.target == ____, ____]
    virginica = iris.data[iris.target == ____, ____]
    print(____)
    print(____)
    print(____)
    print(____)
    f_stat, p_value = stats.f_oneway(____)
    print(____)
    print(____)
    print(____)
    plt.figure(figsize=(____, ____))
    plt.boxplot(____)
    plt.ylabel(____)
    plt.title(____)
    plt.grid(____)
    plt.show()