绘制F分布

def solve():
    import numpy as np
    from sklearn.datasets import load_iris
    iris = load_iris()
    y = iris.target
    p_hat = np.mean(y == 0)
    n = len(y)
    import scipy.stats as stats
    se = np.sqrt(p_hat*(1-p_hat)/n)
    ci = stats.norm.interval(0.95, loc=p_hat, scale=se)
    print(f'{ci[0]:.4f},{ci[1]:.4f}')