预测新观测值

def solve():
    import pandas as pd
    import statsmodels.api as sm
    import matplotlib.pyplot as plt
    import numpy as np
    
    df = pd.read_csv(____)
    X = sm.add_constant(____)
    model = sm.OLS(____).fit()
    
    pred = model.predict(____)
    
    plt.figure(figsize=(____,____))
    plt.scatter(____)
    plt.scatter(____)
    
    plt.title(____)
    plt.xlabel(____)
    plt.ylabel(____)
    plt.legend()
    plt.grid(____)
    plt.tight_layout()
    plt.show()
    
    print(round(____))