大样本模拟-一致性

def solve():
    import numpy as np
    import statsmodels.api as sm
    import matplotlib.pyplot as plt
    
    np.random.seed(____)
    n = ____
    x = np.random.randn(____)
    y = ____ + ____*x + np.random.randn(____)
    X = sm.add_constant(____)
    model = sm.OLS(____).fit()
    
    sample_sizes = [____, ____, ____, ____, ____]
    betas = []
    for ns in sample_sizes:
        x_temp = np.random.randn(____)
        y_temp = ____ + ____*x_temp + np.random.randn(____)
        X_temp = sm.add_constant(____)
        betas.append(sm.OLS(____).fit().params[____])
    
    plt.figure(figsize=(____,____))
    plt.plot(____)
    plt.axhline(____)
    plt.title(____)
    plt.xlabel(____)
    plt.ylabel(____)
    plt.legend()
    plt.grid(____)
    plt.tight_layout()
    plt.show()
    
    print(round(____))