So extrahieren sie p-werte aus der linearen regression in statistischen modellen

Sie können die folgenden Methoden verwenden, um p-Werte für Koeffizienten in einer linearen Regressionsmodellanpassung mithilfe des Statsmodels- Moduls in Python zu extrahieren:

 #extract p-values for all predictor variables
for x in range(0, 3):
    print ( model.pvalues [x])

#extract p-value for specific predictor variable name
model. pvalues . loc [' predictor1 ']

#extract p-value for specific predictor variable position
model. pvalues [0]

Die folgenden Beispiele zeigen, wie die einzelnen Methoden in der Praxis angewendet werden.

Beispiel: Extrahieren Sie P-Werte aus der linearen Regression in statistischen Modellen

Angenommen, wir haben den folgenden Pandas-DataFrame, der Informationen über die gelernten Stunden, die abgelegten Vorbereitungsprüfungen und die Abschlussnote der Schüler einer bestimmten Klasse enthält:

 import pandas as pd

#createDataFrame
df = pd. DataFrame ({' hours ': [1, 2, 2, 4, 2, 1, 5, 4, 2, 4, 4, 3, 6],
                   ' exams ': [1, 3, 3, 5, 2, 2, 1, 1, 0, 3, 4, 3, 2],
                   ' score ': [76, 78, 85, 88, 72, 69, 94, 94, 88, 92, 90, 75, 96]})

#view head of DataFrame
df. head ()

	hours exam score
0 1 1 76
1 2 3 78
2 2 3 85
3 4 5 88
4 2 2 72

Wir können die Funktion OLS() des Moduls statsmodels verwenden, um ein multiples lineares Regressionsmodell anzupassen, indem wir „Stunden“ und „Prüfungen“ als Prädiktorvariablen und „Punktzahl“ als Antwortvariable verwenden:

 import statsmodels. api as sm

#define predictor and response variables
y = df['score']
x = df[['hours', 'exams']]

#add constant to predictor variables
x = sm. add_constant (x)

#fit linear regression model
model = sm. OLS (y,x). fit ()

#view model summary
print ( model.summary ())

                            OLS Regression Results                            
==================================================== ============================
Dept. Variable: R-squared score: 0.718
Model: OLS Adj. R-squared: 0.661
Method: Least Squares F-statistic: 12.70
Date: Fri, 05 Aug 2022 Prob (F-statistic): 0.00180
Time: 09:24:38 Log-Likelihood: -38.618
No. Observations: 13 AIC: 83.24
Df Residuals: 10 BIC: 84.93
Df Model: 2                                         
Covariance Type: non-robust                                         
==================================================== ============================
                 coef std err t P>|t| [0.025 0.975]
-------------------------------------------------- ----------------------------
const 71.4048 4.001 17.847 0.000 62.490 80.319
hours 5.1275 1.018 5.038 0.001 2.860 7.395
exams -1.2121 1.147 -1.057 0.315 -3.768 1.344
==================================================== ============================
Omnibus: 1,103 Durbin-Watson: 1,248
Prob(Omnibus): 0.576 Jarque-Bera (JB): 0.803
Skew: -0.289 Prob(JB): 0.669
Kurtosis: 1.928 Cond. No. 11.7
==================================================== ============================

Standardmäßig zeigt die Funktion summary() die p-Werte jeder Prädiktorvariablen bis zu drei Dezimalstellen an:

• P-Wert für Schnittpunkt: 0,000
• P-Wert für Stunden: 0,001
• P-Wert für Prüfungen: 0,315

Wir können jedoch die vollständigen p-Werte für jede Prädiktorvariable mithilfe der folgenden Syntax aus dem Modell extrahieren:

 #extract p-values for all predictor variables
for x in range(0, 3):
    print ( model.pvalues [x])

6.514115622692573e-09
0.0005077783375870773
0.3154807854805659

Dadurch können wir p-Werte mit mehr Dezimalstellen sehen:

• P-Wert für Schnittpunkt: 0,00000000651411562269257
• P-Wert für Stunden: 0,0005077783375870773
• P-Wert für Prüfungen: 0,3154807854805659

Hinweis : Wir haben in unserer Funktion range() den Wert 3 verwendet, da unser Regressionsmodell drei Gesamtkoeffizienten enthielt.

Wir können auch die folgende Syntax verwenden, um den p-Wert speziell für die Variable „Stunden“ zu extrahieren:

 #extract p-value for 'hours' only
model. pvalues . loc [' hours ']

0.0005077783375870773

Oder wir könnten die folgende Syntax verwenden, um den p-Wert des Koeffizienten einer Variablen an einer bestimmten Position des Regressionsmodells zu extrahieren:

 #extract p-value for coefficient in index position 0
model. pvalues [0]

6.514115622692573e-09

Zusätzliche Ressourcen

Die folgenden Tutorials erklären, wie Sie andere häufige Aufgaben in Python ausführen:

So führen Sie eine logistische Regression in Python durch
So berechnen Sie den AIC von Regressionsmodellen in Python
So berechnen Sie das angepasste R-Quadrat in Python