置信区间覆盖率可视化

def solve():
    import pandas as pd
    import numpy as np
    import matplotlib.pyplot as plt
    from scipy import stats
    np.random.seed(____)
    url = '____'
    df = pd.read_csv(____)
    mu = df['____'].mean()
    lows, highs, covers = [], [], []
    for _ in range(____):
        s = df['____'].sample(____)
        ci = stats.t.interval(____, df=____, loc=s.mean(), scale=stats.sem(____))
        lows.append(____); highs.append(____)
        covers.append(____)
    print(____)
    plt.figure(figsize=(____,____))
    for i,(____) in enumerate(zip(____)):
        plt.plot(____)
    plt.axvline(____)
    plt.title(____); plt.tight_layout(); plt.show()