遗漏变量偏误

def solve():
    import pandas as pd
    import statsmodels.api as sm
    import matplotlib.pyplot as plt
    import numpy as np
    
    df = pd.read_csv(____)
    
    X1 = sm.add_constant(____)
    model1 = sm.OLS(____).fit()
    X2 = sm.add_constant(____)
    model2 = sm.OLS(____).fit()
    
    plt.figure(figsize=(____,____))
    models = ['____', '____']
    coefs = [model1.params['____'], model2.params['____']]
    
    x_pos = np.arange(len(____))
    plt.bar(____)
    plt.xticks(____)
    plt.title(____)
    plt.ylabel(____)
    for i, v in enumerate(____):
        plt.text(____)
    plt.grid(____)
    plt.tight_layout()
    plt.show()
    
    diff = model1.params['____'] - model2.params['____']
    print(round(____))