Проблемы с OLS для Python Pandas Statsmodel

Question

Я написал следующую функцию.Намерение состоит в том, чтобы получить модель, которая дает правильный набор функций с pvalues ​​<0,05. Функция не дает согласованных результатов.В крайних случаях это дает результаты, где значение больше 0,9. Может ли кто-нибудь помочь мне здесь? </p>

def run_reg(X,Y,df):
    #Now lets start the regression model
    lreg = sm.OLS(Y,X)
    feature_count=1
    while(feature_count !=0):
        model = lreg.fit()
        pvalues = model.pvalues
        print(pvalues)
        feature_count= len(pvalues[pvalues > 0.05] )
        print(pvalues)
        if(feature_count != 0):
            pvalues = pvalues[pvalues < 0.05]
            pvalues = pvalues.drop('const',axis=0)
            X = df[pvalues.index]
            X = sm.add_constant(X)
            lreg = run_reg(X,Y)
    return model