كيفية رسم منحنى الانحدار متعدد الحدود في r

الانحدار متعدد الحدود هو أسلوب انحدار نستخدمه عندما تكون العلاقة بين متغير التوقع ومتغير الاستجابة غير خطية.

يشرح هذا البرنامج التعليمي كيفية رسم منحنى الانحدار متعدد الحدود في R.

ذات صلة: الأنواع السبعة الأكثر شيوعًا للانحدار

مثال: رسم منحنى انحدار متعدد الحدود في R

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

 #define data
x <- runif(50, 5, 15)
y <- 0.1*x^3 - 0.5 * x^2 - x + 5 + rnorm(length(x),0,10) 
 
#plot x vs. y
plot(x, y, pch= 16 , cex= 1.5 ) 
 
#fit polynomial regression model
fit <- lm(y ~ x + I(x^2) + I(x^3))
 
#use model to get predicted values
pred <- predict(fit)
ix <- sort(x, index. return = T )$ix

#add polynomial curve to plot
lines(x[ix], pred[ix], col=' red ', lwd= 2 )

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

 #define data
x <- runif(50, 5, 15)
y <- 0.1*x^3 - 0.5 * x^2 - x + 5 + rnorm(length(x),0,10) 
 
#plot x vs. y
plot(x, y, pch=16, cex=1.5) 
 
#fit polynomial regression model
fit <- lm(y ~ x + I(x^2) + I(x^3))
 
#use model to get predicted values
pred <- predict(fit)
ix <- sort(x, index. return = T )$ix

#add polynomial curve to plot
lines(x[ix], pred[ix], col=' red ', lwd= 2 )

#get model coefficients
coeff <- round(fit$coefficients, 2)

#add fitted model equation to plot
text(9, 200 , paste("Model: ", coeff[1], " + ", coeff[2],
                    "*x", "+", coeff[3], "*x^2", "+", coeff[4], "*x^3"), cex= 1.3 )

لاحظ أن الوسيطة cex تتحكم في حجم خط النص. القيمة الافتراضية هي 1، لذلك اخترنا استخدام القيمة 1.3 لتسهيل قراءة النص.

مصادر إضافية

مقدمة للانحدار متعدد الحدود
كيفية احتواء منحنى متعدد الحدود في إكسيل
كيفية إجراء الانحدار متعدد الحدود في بايثون