主题：【原创】人类不感谢罗辑--谁是下一位老中医 -- 泉畔人家

library(ggpubr)

GDP <- read.table("https://www.dropbox.com/s/a5ljh29ebo82fld/GDP.txt?raw=1",header = T)

GDP$Log10.GDP = log10(GDP$GDP)

g <- ggscatter(GDP, x = "Year", y = "Log10.GDP",

title = "China GDP since year 1949",

conf.int = TRUE, # Add confidence interval

color = "before.after.78", palette = "jco")+

geom_smooth(aes(color = "before.after.78"), method = lm)+

stat_cor(aes(label = paste(..rr.label.., ..p.label.., sep = "~`,`~"),

color = "before.after.78"))+

theme(plot.title = element_text(hjust = 0.5),

legend.position = "none")

g

原始数据直接复制泉畔人家的帖子，稍作整理。