วิธีการพล็อตเส้นการถดถอยเชิงเส้นใน ggplot2 (พร้อมตัวอย่าง)

คุณสามารถใช้ไลบรารีการแสดงภาพ R ggplot2 เพื่อพล็อตโมเดลการถดถอยเชิงเส้นที่ติดตั้งโดยใช้ไวยากรณ์พื้นฐานต่อไปนี้:

 ggplot(data,aes(x, y)) +
  geom_point() +
  geom_smooth(method=' lm ')

ตัวอย่างต่อไปนี้แสดงวิธีใช้ไวยากรณ์นี้ในทางปฏิบัติ

ตัวอย่าง: การวางแผนเส้นการถดถอยเชิงเส้นใน ggplot2

สมมติว่าเราใส่แบบจำลอง การถดถอยเชิงเส้นอย่างง่าย เข้ากับชุดข้อมูลต่อไปนี้:

 #create dataset
data <- data.frame(y=c(6, 7, 7, 9, 12, 13, 13, 15, 16, 19, 22, 23, 23, 25, 26),
                   x=c(1, 2, 2, 3, 4, 4, 5, 6, 6, 8, 9, 9, 11, 12, 12))

#fit linear regression model to dataset and view model summary
model <- lm(y~x, data=data)
summary(model)

Call:
lm(formula = y ~ x, data = data)

Residuals:
    Min 1Q Median 3Q Max 
-1.4444 -0.8013 -0.2426 0.5978 2.2363 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 4.20041 0.56730 7.404 5.16e-06 ***
x 1.84036 0.07857 23.423 5.13e-12 ***
---
Significant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 1.091 on 13 degrees of freedom
Multiple R-squared: 0.9769, Adjusted R-squared: 0.9751 
F-statistic: 548.7 on 1 and 13 DF, p-value: 5.13e-12

รหัสต่อไปนี้แสดงวิธีแสดงภาพแบบจำลองการถดถอยเชิงเส้นที่ติดตั้ง:

 library (ggplot2)

#create plot to visualize fitted linear regression model
ggplot(data,aes(x, y)) +
  geom_point() +
  geom_smooth(method=' lm ') 

ตามค่าเริ่มต้น ggplot2 จะเพิ่มบรรทัดข้อผิดพลาดมาตรฐานลงในกราฟ คุณสามารถปิดการใช้งานได้โดยใช้อาร์กิวเมนต์ se=FALSE ดังนี้:

 library (ggplot2)

#create regression plot with no standard error lines
ggplot(data,aes(x, y)) +
  geom_point() +
  geom_smooth(method=' lm ', se= FALSE ) 

สุดท้ายนี้ เราสามารถปรับแต่งบางแง่มุมของแผนภูมิเพื่อให้ดูน่าดึงดูดยิ่งขึ้นได้:

 library (ggplot2)

#create regression plot with customized style
ggplot(data,aes(x, y)) +
  geom_point() +
  geom_smooth(method=' lm ', se= FALSE , color=' turquoise4 ') +
  theme_minimal() +
  labs(x=' X Values ', y=' Y Values ', title=' Linear Regression Plot ') +
  theme(plot.title = element_text(hjust=0.5, size=20, face=' bold ')) 

โปรดดู บทความนี้ เพื่อดูคำแนะนำฉบับสมบูรณ์เกี่ยวกับธีม ggplot2 ที่ดีที่สุด

แหล่งข้อมูลเพิ่มเติม

ความรู้เบื้องต้นเกี่ยวกับการถดถอยเชิงเส้นพหุคูณใน R
วิธีการพล็อตช่วงความเชื่อมั่นใน R