可变膨胀因子不适用于dataframes python

import pandas as pd
import statsmodels.api as sm

# 假设你有一个名为df的数据框，包含自变量和因变量
# ...

# 添加常数列作为回归模型的截距
df['intercept'] = 1

# 使用statsmodels进行回归分析
model = sm.OLS(df['因变量'], df[['自变量1', '自变量2', 'intercept']])
results = model.fit()

# 计算可变膨胀因子
vif = pd.DataFrame()
vif["特征"] = df[['自变量1', '自变量2', 'intercept']].columns
vif["VIF因子"] = [variance_inflation_factor(df[['自变量1', '自变量2', 'intercept']].values, i) for i in range(df[['自变量1', '自变量2', 'intercept']].shape[1])]

# 打印可变膨胀因子
print(vif)