ANOVA方差分解表

def solve():
    import pandas as pd
    import numpy as np
    import matplotlib.pyplot as plt
    from scipy.stats import f_oneway
    url = '____'
    df = pd.read_csv(____)
    gs = [df[df['____']==g]['____'].dropna() for g in [____,____]]
    grand = df['____'].mean()
    N = sum(len(____) for g in gs)
    SSB = sum(len(____)*(g.mean()-grand)**____ for g in gs)
    SSW = sum(((g-g.mean())**____).sum() for g in gs)
    MSB, MSW = SSB/____, SSW/(____)
    F = MSB/MSW
    _, p = f_oneway(____)
    print(____)
    plt.boxplot(gs, labels=['____','____'], patch_artist=True, boxprops=dict(____))
    plt.ylabel(____); plt.title(____); plt.tight_layout(); plt.show()