R语言ggplot2折线图(line plot)添加置信区间（CI）展示学术论文作者数量的变化趋势

用户7010445

发布于 2021-05-07 10:53:15

7.5K1

发布于 2021-05-07 10:53:15

文章被收录于专栏：小明的数据分析笔记本小明的数据分析笔记本

非常有意思的数据可视化案例，原文提出的问题是 学术论文中的作者数量有逐年增加的趋势 ；于是利用R语言里的 rplos 包抓取了 Plos 系列的6本期刊的2006年至2013年的每篇论文里的作者数量进行可视化展示

原文链接是 https://benjaminlmoore.wordpress.com/2014/04/06/author-inflation-in-academic-literature/

完整的代码链接

https://github.com/blmoore/blogR

原始代码中抓取数据的部分好像不能用了，我稍微改动了一下，选取了2006年到2020年的数据，获取数据的代码这里就不放了，如果需要本文的示例数据可以知己在文末留言

数据已经存储到了文件里首先是读取数据

library(readr)
df<-read_tsv("author_number.tsv")
head(df)

image.png

期刊分别是

table(df$journal)

image.png

作图代码

ibrary(ggplot2)
ggplot(df, aes(x=year, y=auth_num, col=journal, fill=journal)) +
  stat_summary(fun.data="mean_cl_boot", geom="ribbon",
               #width=.2, 
               alpha=I(.5)) +
  stat_summary(fun="mean", geom="line") +
  labs(x="Year", y="Mean number of authors per paper") +
  theme_bw() +
  theme(legend.position="top") +
  scale_fill_brewer(type="qual", palette=2, 
                    guide=guide_legend(direction="vertical",
                                       label.position="bottom",
                                       title=NULL, 
                                       ncol=6,
                                       label.hjust=0.5)) +
  scale_color_brewer(type="qual", palette=2, guide="none")+
  facet_wrap(~journal,ncol=23)