1995～2023 年上市公司参控股公司专利创新突破度 CD 指数计算结果

1277 字 约 4 分钟 数据资料 专利上市公司创新与知识

《Nature》论文：「Papers and patents are becoming less disruptive over time」中使用了专利创新突破度（CD 指数）数据：

CD 指数的计算思路如下：对于专利 i，首先找出其后向引用（即专利 i 引用了哪些更早期的专利），然后考察后续专利 j（申请日期晚于专利 i）的引用行为：

  1. 若专利 j 既引用了专利 i，又引用了专利 i 的后向引用专利，说明专利 j 沿用了专利 i 的既有知识脉络，此时专利 i 的该项引用 CD 值为 -1；
  2. 若专利 j 只引用专利 i 但未引用其后向引用专利，说明专利 i 具有知识突破性，CD 值为 1；
  3. 若专利 j 没有引用专利 i 但引用其后向引用专利，CD 值为 0。

将所有后续专利对专利 i 的引用 CD 值取平均，即得到专利 i 的 CD 指数。

关于 CD 指数的计算方法，之前也讲解过：

名师讲堂｜使用 R 语言测算专利的创新突破度：CD 指数：https://rstata.duanshu.com/#/course/3b00a43418ee43f194df7fa664ef0c7d

名师讲堂｜使用 Stata 测算专利的创新突破度：CD指数：https://rstata.duanshu.com/#/course/ca92eb2b71854c1b8a6cd74c66bd6c5e

今天再给大家分享一份 1995～2023 年上市公司参控股公司专利创新突破度（CD 指数）的计算结果。

数据概览

该数据由上市公司参控股公司专利申请与授权信息匹配全部专利的引用与被引用关系计算得到。

主要变量说明：

变量名 说明
ckggsid 参控股公司数据观测值编号
年份 会计年度（1995-2023）
newipzlid 专利数据观测值编号
cd_index CD指数（范围 -1 到 1）

图表展示

下图展示了历年上市公司参控股公司平均专利创新突破度（CD 指数）变化趋势：

下图展示了 CD 指数的分布特征：

按时间段分组的箱线图显示：

观测数量变化趋势：

各年度平均CD指数对比：

CD指数累积分布函数：

创新类型占比变化（破坏性/渐进性/突破性）：

各年度CD指数分布详情：

带线性趋势的散点图：

处理代码

为方便大家学习，附件中也提供了数据处理代码以供参考：

#*- 读取上市公司参控股专利数据
readr::read_csv("/Volumes/ADT/常用大数据/1995~2023年上市公司参控股公司专利申请与授权信息.csv") -> df1

df1 %>%
  left_join(unique_newipzlid %>% mutate(value = 1)) -> df1b

#*- 根据公开公告号去重
df1b %>%
  mutate(年份 = 会计年度) %>%
  mutate(公开公告号 = str_remove(公开公告号, "[A-Z]$")) %>%
  mutate(value = if_else(is.na(value), 0, value)) %>%
  group_by(ckggsid, 年份, 公开公告号) %>%
  filter(value == max(value, na.rm = T)) %>%
  ungroup() %>%
  distinct(ckggsid, 年份, 公开公告号, .keep_all = T) %>%
  select(-公开公告号) %>%
  #*- 再根据申请号去重
  group_by(ckggsid, 年份, 申请号) %>%
  filter(value == max(value, na.rm = T)) %>%
  ungroup() %>%
  distinct(ckggsid, 年份, 申请号, .keep_all = T) %>%
  select(-申请号, -value) %>%
  select(ckggsid, newipzlid, 年份) -> patent_info

#*- 和CD指数计算结果链接
patent_info %>%
  left_join(bind_rows(df0, nobackrefs)) -> dfres

dfres %>%
  mutate(cd_index = if_else(is.na(cd_index), 0, cd_index)) -> dfres

#*- 保存主数据文件
dfres %>%
  haven::write_dta("1995～2023年上市公司参控股公司专利创新突破度 CD 指数计算结果.dta",
                   label = "数据处理：微信公众号 RStata")

#*- 保存分公司年份平均数据
dfres %>%
  group_by(ckggsid, 年份) %>%
  summarise(cd_index = mean(cd_index, na.rm = T), .groups = "drop") %>%
  haven::write_dta("1995～2023年上市公司参控股公司专利创新突破度 CD 指数计算结果（分公司年份平均）.dta",
                   label = "数据处理：微信公众号 RStata")

附件中也提供了该数据的处理代码供参考：

数据引用格式

由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：

RStata 数据中心: 1995～2023年上市公司参控股公司专利创新突破度 CD 指数计算结果. 2026. https://tidyfriday.cn/rsdb2/

英文文献可以使用下面的格式引用：

RStata Data Center: CD Index of Patent Innovation Breakthrough for Listed Companies’ Subsidiaries and Controlled Companies, 1995–2023. 2026. https://tidyfriday.cn/rsdb2/

关联课程/数据推荐

名师讲堂｜使用 R 语言测算专利的创新突破度：CD 指数：https://rstata.duanshu.com/#/course/3b00a43418ee43f194df7fa664ef0c7d

名师讲堂｜使用 Stata 测算专利的创新突破度：CD指数：https://rstata.duanshu.com/#/course/ca92eb2b71854c1b8a6cd74c66bd6c5e

1985～2024年上市公司渐进性创新与突破性创新相关指标计算结果：https://rstata.duanshu.com/#/course/786e14a581564003b55f8bf340cfd205

1985～2024 年各城市渐进性创新与突破性创新相关指标计算结果：https://rstata.duanshu.com/#/course/601a1c3da3954ca39283e5e85d9a59d3

1995～2023 年上市公司参控股公司专利申请与授权信息：https://rstata.duanshu.com/#/course/e5e68a65a0234f18b1338ad1caadf33d

如果有需要 RStata 定制相关数据的，可以联系李老师付费定制。

点击这里跳转到 RStata 短书平台获取附件：1995～2023 年上市公司参控股公司专利创新突破度 CD 指数计算结果

相关文章

评论