分位数回归（含比较图）

def solve():
    import pandas as pd
    import statsmodels.api as sm
    import matplotlib.pyplot as plt
    
    df = pd.read_csv(____)
    X = sm.add_constant(____)
    
    quantiles = [____, ____, ____]
    coefs = []
    for q in quantiles:
        model_q = sm.QuantReg(____).fit(____)
        coefs.append(____)
    
    plt.figure(figsize=(____,____))
    plt.bar(____)
    plt.title(____)
    plt.xlabel(____)
    plt.ylabel(____)
    plt.grid(____)
    plt.tight_layout()
    plt.show()
    
    model = sm.QuantReg(____).fit(____)
    print(____)