1998～2013 年每年各种专利类型的申请数量和专利总数量的绘图展示：R 语言 & Stata
在之前的推文 1998～2013年工企与专利匹配、各工企专利申请数量等数据（大数据！89 GB + ） 中我绘制了这样的一幅图：
很多小伙伴感觉这幅图很漂亮，好奇这个图是怎么绘制的。今天就来解答一下！
首先大家可以下载1998～2013年工企与专利匹配、各工企专利申请数量等数据（大数据！88 GB + ） 里面介绍的数据（购买会员后即可从 RStata 平台上下载：https://rstata.duanshu.com/#/course/83e6377f39124f5ca1b513e212b18b32 ），推文中提供了多种格式，下面使用 1998~2013年各工企业专利数量.dta 数据。
首先读取这份数据：
# 加载 tidyverse
library(tidyverse)
# 读取 1998~2013年各工企业专利数量数据
haven::read_dta('1998~2013年各工企业专利数量.dta') -> cdf
下面我们统计每年每种专利的数量：
# 汇总每年每种专利的数量cdf %>% group_by(年份) %>% summarise(设计型专利 = sum(设计型专利), 实用型专利 = sum(实用型专利), 发明型专利 = sum(发明型专利), 总专利数量 = sum(总专利数量)) -> cdfcdf#> # A tibble: 16 x 5#> 年份 设计型专利 实用型专利 发明型专利 总专利数量#> <chr> <dbl> <dbl> <dbl> <dbl>#> 1 1998 5489 3551 770 9810#> 2 1999 8187 5141 1290 14618#> 3 2000 9249 6107 2097 17453#> 4 2001 10275 7611 3411 21297#> 5 2002 14077 10593 6793 31463#> 6 2003 14839 13646 10214 38699#> 7 2004 19771 16984 14883 51638#> 8 2005 22633 20593 20872 64098#> 9 2006 28729 28318 30395 87442#> 10 2007 34249 37106 39442 110797#> 11 2008 45657 57037 50639 153333#> 12 2009 52819 80686 60881 194386#> 13 2010 51062 77831 75249 204142#> 14 2011 72751 128677 109636 311064#> 15 2012 110226 176018 153316 439560#> 16 2013 101224 195871 189561 486656然后就可以使用这个数据绘图了，我使用的是 highcharter 包绘制：
* 安装 plotplain 主题：ssc install blindschemes
use cdf.dta, clear
destring _all, replace
tw bar 总专利数量 年份, barw(0.4) ///
xla(1998(1)2013) color('158 203 193') ///
yaxis(2) || ///
conn 实用型专利 年份, color('55 92 178') m(O) lp(solid) yaxis(1) || ///
conn 设计型专利 年份, color('216 69 53') m(O) lp(solid) yaxis(1) || ///
conn 发明型专利 年份, color('226 180 94') m(O) lp(solid) yaxis(1) ///
leg(order(2 '实用型专利' 3 '设计型专利' 4 '发明型专利') pos(6) row(1)) ///
xti('') yti('各种类型的数量', axis(1)) ///
yti('总数量', axis(2) color('141 211 199')) ///
yla(0 '0' 60000 '60k' 120000 '120k' ///
180000 '180k' 240000 '240k', axis(1)) ///
yla(0 '0' 150000 '150k' 300000 '300k' 450000 '450k' 600000 '600k', axis(2) labcolor('141 211 199') tlcolor('141 211 199')) ///
ysc(lc('141 211 199') axis(2)) ///
ti('工企与专利数据库匹配结果中的专利数量', size(*1.3)) ///
subti('时间范围：1998~2013 年') ///
caption('数据来源：工企数据库，Chinese Patent Data Project', size(*0.8)) ///
note('绘制：微信公众号 RStata', size(*0.8)) ///
scheme(plotplain)