كيفية إنشاء قطعة أرض متبقية في ggplot2 (مع مثال)

تُستخدم المخططات المتبقية لتقييم ما إذا كانت بقايا نموذج الانحدار موزعة بشكل طبيعي وما إذا كانت تظهر تغايرًا أم لا.

لإنشاء قطعة أرض متبقية في ggplot2، يمكنك استخدام بناء الجملة الأساسي التالي:

 library (ggplot2)

ggplot(model, aes(x = .fitted, y = .resid)) +
  geom_point() +
  geom_hline(yintercept = 0 )

يوضح المثال التالي كيفية استخدام بناء الجملة هذا عمليًا.

مثال: إنشاء قطعة أرض متبقية في ggplot2

في هذا المثال، سوف نستخدم مجموعة بيانات mtcars المضمنة في R:

 #view first six rows of mtcars dataset
head(mtcars)
                   mpg cyl disp hp drat wt qsec vs am gear carb
Mazda RX4 21.0 6 160 110 3.90 2.620 16.46 0 1 4 4
Mazda RX4 Wag 21.0 6 160 110 3.90 2.875 17.02 0 1 4 4
Datsun 710 22.8 4 108 93 3.85 2.320 18.61 1 1 4 1
Hornet 4 Drive 21.4 6 258 110 3.08 3.215 19.44 1 0 3 1
Hornet Sportabout 18.7 8 360 175 3.15 3.440 17.02 0 0 3 2
Valiant 18.1 6 225 105 2.76 3,460 20.22 1 0 3 1

أولاً، سنلائم نموذج الانحدار باستخدام mpg كمتغير الاستجابة و qsec كمتغير متوقع:

 #fit regression model
model <- lm(mpg ~ qsec, data=mtcars)

بعد ذلك، سوف نستخدم بناء الجملة التالي لإنشاء قطعة أرض متبقية في ggplot2:

 library (ggplot2)

#create residual plot
ggplot(model, aes(x = .fitted, y = .resid)) +
  geom_point() +
  geom_hline(yintercept = 0 ) 

يعرض المحور السيني القيم المجهزة ويعرض المحور الصادي القيم المتبقية.

ويبدو أن البقايا متناثرة بشكل عشوائي حول الصفر مع عدم وجود نمط واضح، مما يشير إلى أن افتراض التماثل محقق.

بمعنى آخر، يجب أن تكون معاملات نموذج الانحدار موثوقة ولا نحتاج إلى إجراء أي تحويل على البيانات.

لاحظ أيضًا أنه يمكننا استخدام الدالة labs() لإضافة عنوان وتسميات محاور إلى المخطط المتبقي:

 library (ggplot2)

#create residual plot with title and axis labels
ggplot(model, aes(x = .fitted, y = .resid)) +
  geom_point() +
  geom_hline(yintercept = 0 ) +
  labs(title=' Residual vs. Fitted Values Plot ', x=' Fitted Values ', y=' Residuals ') 

مصادر إضافية

تشرح البرامج التعليمية التالية كيفية تنفيذ المهام الشائعة الأخرى في R:

كيفية حساب المخلفات الموحدة في R
كيفية حساب المخلفات الطالبية في R
كيفية إنشاء رسم بياني للبقايا في R