R'de eklenen değişken grafikleri nasıl oluşturulur

İstatistiklerde, eklenen değişken grafikleri, bir çoklu doğrusal regresyon modelinde bir yanıt değişkeni ile bir yordayıcı değişken arasındaki ilişkiyi görüntülerken, modeldeki diğer yordayıcı değişkenlerin varlığını kontrol eden ayrı grafiklerdir.

Not: Bazen bu grafiklere “kısmi regresyon grafikleri” de denir.

Bu tür çizim, diğer öngörücü değişkenleri sabit tutarken, bir modeldeki her bir yordayıcı değişken ile yanıt değişkeni arasındaki ilişkiyi gözlemlememize olanak tanır.

R’de eklenen değişkenlerin grafiklerini oluşturmak için araba paketindeki avPlots() fonksiyonunu kullanabiliriz:

 #load car package
library (car) 

#fit multiple linear regression model
model <- lm(y ~ x1 + x2 + ..., data = df)

#create added variable plots
avPlots(model)

Aşağıdaki örnek, bu sözdiziminin pratikte nasıl kullanılacağını gösterir.

Örnek: R’ye Değişken Grafikler Ekleme

Mtcars veri kümesindeki verileri kullanarak aşağıdaki çoklu doğrusal regresyon modelini R’ye uydurduğumuzu varsayalım:

 #fit multiple linear regression model
model <- lm(mpg ~ disp + hp + drat, data = mtcars)

#view summary of model
summary(model)

Call:
lm(formula = mpg ~ disp + hp + drat, data = mtcars)

Residuals:
    Min 1Q Median 3Q Max 
-5.1225 -1.8454 -0.4456 1.1342 6.4958 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)   
(Intercept) 19.344293 6.370882 3.036 0.00513 **
available -0.019232 0.009371 -2.052 0.04960 * 
hp -0.031229 0.013345 -2.340 0.02663 * 
drat 2.714975 1.487366 1.825 0.07863 . 
---
Significant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 3.008 on 28 degrees of freedom
Multiple R-squared: 0.775, Adjusted R-squared: 0.7509 
F-statistic: 32.15 on 3 and 28 DF, p-value: 3.28e-09

Yanıt değişkeni “mpg” ile modeldeki her bir tahmin değişkeni arasındaki ilişkiyi görselleştirmek için avPlots() işlevini kullanarak eklenen değişkenlerin grafiklerini üretebiliriz:

 #load car package
library (car)

#produce added variable plots
avPlots(model)

Her grafiği nasıl yorumlayacağınız aşağıda açıklanmıştır:

• X ekseni tek bir öngörücü değişkeni görüntüler ve y ekseni yanıt değişkenini görüntüler.
• Mavi çizgi , diğer tüm öngörücü değişkenlerin değerini sabit tutarken , yordayıcı değişken ile yanıt değişkeni arasındaki ilişkiyi gösterir.
• Her grafikteki etiketli noktalar, en büyük artıklara sahip iki gözlemi ve en büyük kısmi kaldıraca sahip iki gözlemi temsil eder.

Her grafikteki çizginin açısının, tahmin edilen regresyon denkleminin katsayısının işaretine karşılık geldiğine dikkat edin.

Örneğin, modeldeki her bir yordayıcı değişken için tahmin edilen katsayılar şunlardır:

• ekran: -0.019232
• kanal: -0.031229
• tarih: 2.714975

Drat için eklenen değişken grafiğinde çizginin açısının pozitif olduğunu, disp ve hp için ise negatif olduğunu ve bunların tahmini katsayılarının işaretlerine karşılık geldiğini unutmayın:

Bu grafikler, her bir öngörücü değişken ile yanıt değişkeni arasındaki ilişkiyi kolayca görselleştirmemize olanak tanır.

Ek kaynaklar

R’de basit doğrusal regresyon nasıl gerçekleştirilir
R’de çoklu doğrusal regresyon nasıl gerçekleştirilir
R’de lojistik regresyon nasıl gerçekleştirilir