多重共线性检验VIF

def solve():
    import pandas as pd
    import matplotlib.pyplot as plt
    from statsmodels.stats.outliers_influence import variance_inflation_factor
    import numpy as np
    url = '____'
    df = pd.read_csv(____)
    X = df[['____','____','____']].dropna().values
    Xa = np.column_stack([np.ones(len(____)), X])
    vifs = []
    for i,c in enumerate(____):
        vif = variance_inflation_factor(____)
        vifs.append(____)
        print(____)
    plt.bar(____)
    plt.axhline(____, color='____', ls='____', label='____')
    plt.ylabel(____); plt.title(____); plt.legend()
    plt.tight_layout(); plt.show()