简单线性回归

def solve():
    from sklearn.datasets import load_diabetes
    import statsmodels.api as sm
    data = load_diabetes()
    X = sm.add_constant(____)
    y = data.target
    model = sm.OLS(____).fit()
    print(____)