Scatterplot with regression line

こういう図を描く

# 練習データ
data <- matrix(c(1001,34,31,1002,13,4,1003,41,39,1004,37,36,1005,40,38,
                 1006,12,14,1007,48,41,1008,46,50,1009,38,44,1010,38,29,
                 1011,17,13,1012,20,25,1013,12,6,1014,17,15,1015,22,15,
                 1016,19,10,1017,14,5,1018,32,26,1019,41,37,1020,16,12,
                 1021,34,28,1022,25,20,1023,29,33,1024,39,34,1025,33,24,
                 1026,32,23,1027,38,41,1028,43,43,1029,28,21,1030,38,35,
                 1031,41,31,1032,48,50,1033,22,25,1034,37,34,1035,16,20,
                 1036,50,47,1037,12,17,1038,46,50,1039,49,45,1040,20,10,
                 2001,43,44,2002,19,27,2003,35,36,2004,43,41,2005,17,20,
                 2006,16,24,2007,49,45,2008,12,20,2009,30,32,2010,49,45,
                 2011,17,26,2012,14,23,2013,28,40,2014,20,31,2015,10,28,
                 2016,45,49,2017,25,26,2018,43,40,2019,31,32,2020,25,37,
                 2021,44,50,2022,20,23,2023,24,25,2024,28,33,2025,17,24,
                 2026,31,32,2027,22,20,2028,34,40,2029,19,28,2030,25,22,
                 2031,33,36,2032,48,44,2033,42,42,2034,46,45,2035,32,34,
                 2036,10,29,2037,45,45,2038,35,41,2039,12,23,2040,28,24),
               nrow = 80, ncol = 3, byrow = TRUE)

data.df <- data.frame(data)
colnames(data.df) <- c("id", "x", "y")
DT::datatable(data.df)

# 散布図
plot(0, 0, type = "n", xlim = c(0, 60), ylim = c(0, 60),
     xlab = "X", ylab = "Y")
points(data.df$x, data.df$y, col = 1, pch = 1)

# 回帰直線
lm01 <- lm(data.df$y ~ data.df$x)
abline(lm01, col=1, lty=1)

# 画像として保存する
quartz(type = "png", file = "./Plot/Scatter_Abline.png", 
       width = 5, height = 5) 
par(family = "HiraKakuProN-W3")
plot(0, 0, type = "n", xlim = c(0, 60), ylim = c(0, 60),
     xlab = "X", ylab = "Y")
points(data.df$x, data.df$y, col = 1, pch = 1)

# 回帰直線
lm01 <- lm(data.df$y ~ data.df$x)
abline(lm01, col=1, lty=1)
dev.off()

## quartz_off_screen 
##                 2

出力したpngファイル