多元线性回归

def solve():
    import numpy as np
    from sklearn.datasets import load_diabetes
    from sklearn.preprocessing import PolynomialFeatures
    from sklearn.linear_model import LinearRegression
    from sklearn.metrics import r2_score
    data = load_diabetes()
    X = data.data[:, ____].reshape(____)
    y = data.target
    poly = PolynomialFeatures(____)
    X_poly = poly.fit_transform(____)
    model = LinearRegression().fit(____)
    r2 = r2_score(y, model.predict(____))
    print(____)