名师讲堂｜使用 R 语言计算专利文本相似度及城市间知识溢出指标

⚠️注意：如果不是走投无路了，非常不建议使用这种方法，运算量过大。

引言

知识溢出（Knowledge Spillover）是区域创新研究的核心概念，指一个地区的知识活动对其他地区创新产出产生的外部性影响。准确测度城市间的知识溢出，对于理解区域创新网络、制定创新政策具有重要意义。

传统的知识溢出测度方法主要依赖地理距离、专利申请引用等间接指标，但这些方法存在明显局限：

• 地理距离无法捕捉非空间的知识传播。
• 专利引用存在严重的漏报和时滞问题
• 无法反映知识内容的相似性

随着文本挖掘技术的发展，专利摘要文本相似度为直接测度知识溢出提供了新的可能。本文档介绍如何使用 TF-IDF（Term Frequency-Inverse Document Frequency） 方法计算专利文本相似度，并构建城市间知识溢出指标。

1. TF-IDF 的原理

TF-IDF 是一种用于信息检索和文本挖掘的常用加权技术，由两个部分组成：

为什么 TF-IDF 适合专利文本相似度计算？

2.1 自动降噪：过滤无意义高频词

专利摘要中包含大量通用表述（如”本发明涉及…”、”其特征在于…”），这些词在几乎所有专利中都会出现。TF-IDF 通过 IDF 权重自动降低这些”停用词”的贡献，突出真正有区分度的技术术语。

对比：单纯的词频（TF）或二进制向量（词是否出现）会将”的”和”深度学习”同等对待，显然不合理。

2.2 区分通用技术与专有技术

• 通用技术词（如”控制器”、”系统”）：IDF 值低，对相似度贡献小
• 专有技术词（如”卷积神经网络”、”石墨烯”）：IDF 值高，对相似度贡献大

这种加权方式恰好符合专利相似度判断的需求：关注技术内容的相似性，而非表述格式的相似性。

2.3 稀疏性友好：适合大规模计算

TF-IDF 向量通常是高维稀疏向量（数万维特征，但每个专利只有数百个非零值）。这种稀疏性使得：

• 存储效率高（可用 Matrix::sparseMatrix 存储）
• 计算速度快（稀疏矩阵乘法优化）
• 内存占用小（适合百万级专利数据）

2.4 余弦相似度的自然适配

TF-IDF 向量经过 L2 归一化（每个向量的欧氏长度变为1）后，两个向量的点积就等于它们的余弦相似度：

这一性质使得相似度计算可以转化为高效的矩阵乘法操作，避免了逐一计算余弦相似度的高昂成本。

3. 与其他方法的对比

方法	优点	缺点	适用场景
TF-IDF + 余弦相似度	计算效率高，可解释性强，适合大规模数据	无法捕捉语义相似（同义词、上下位关系）	初步筛选、大规模专利相似度计算
Word2Vec/Doc2Vec	捕捉语义相似，考虑词序信息	需要大量训练数据，计算成本高	小规模精细分析、语义相似度研究
BERT/Transformer	最先进的语义理解能力	计算成本极高，需要GPU	高精度需求、小规模案例研究

结论：对于数百万级专利对的知识溢出计算，TF-IDF 是在计算效率和准确度之间的最佳平衡点。

4. 论文中的谷歌专利向量方法

4.1 什么是谷歌专利文本向量？

谷歌专利文本向量（Google Patent Vectors） 是由 Google Patents 团队基于海量专利全文文本，利用机器学习模型生成的高维语义向量。

技术特点：

• 基于专利全文（而非仅摘要）训练，技术信息更全面
• 使用深度学习模型捕捉语义相似度（如同义词、上下位关系）
• 向量维度高（通常 512 或 1024 维），语义表示能力强
• 已公开供科学研究使用

4.2 为什么原论文使用谷歌专利向量？

王慧扬等（2026）的论文选择谷歌专利向量，主要基于以下考虑：

1. 专利摘要信息有限：与专利代理人访谈发现，专利摘要传达的技术信息十分有限，无法全面反映技术方案
2. 全文文本获取困难：逐一获取数百万件专利的全文文本不具有可行性
3. 模型训练成本高：生成高质量专利文本向量需要基于数以百万计专利语料训练的机器学习模型，对资金和技术要求高
4. 谷歌已提供公开数据：Google Patents 团队已经完成了上述工作，并将专利文本向量公开供科学研究使用
5. 已有研究验证有效性：Higham 等（2021）、De Rassenfosse 和 Palangkaraya（2023）等研究已使用谷歌专利向量并验证了其有效性

4.3 为什么本实现没有使用谷歌专利向量？

尽管谷歌专利向量在语义理解上具有优势，但本实现选择 TF-IDF 方法，主要基于以下原因：

考量维度	谷歌专利向量	TF-IDF	本实现的选择理由
数据可得性	需要匹配中国专利与谷歌数据库	直接使用专利摘要文本	中国专利与谷歌数据库的匹配存在困难
数据覆盖范围	主要覆盖美国等发达国家专利	适用于任何文本的专利	中国专利的覆盖率可能不足
方法透明度	黑盒模型（深度学习）	白盒模型（统计方法）	学术研究需要可解释性
计算资源需求	需要下载海量向量数据（GB 级）	可实时计算	降低计算门槛
中文处理适应性	主要针对英文专利训练	直接适用于中文分词	避免跨语言语义表示的偏差
方法简洁性	依赖外部数据和 API	完全自包含	便于方法复现和推广

核心原因总结：

1. 数据匹配挑战：将中国专利号（如 CN107519921B）与 Google Patents 的专利 ID 进行匹配需要额外的工作，且覆盖率不确定
2. 数据获取成本：谷歌专利向量数据量庞大（数 GB），下载和存储成本高
3. 中国专利适应性：谷歌专利向量主要基于英文专利训练，对中文专利的语义表示效果有待验证（原论文也指出需要在发展中国家情境中检验其有效性）
4. 方法可复现性：TF-IDF 是完全自包含的统计方法，不依赖任何外部数据或专有模型，更利于学术研究的开放性和复现性
5. 计算效率满足需求：对于城市间知识溢出这种大规模计算任务，TF-IDF 的效率优势更为明显

4.4 两种方法的结果对比

如果使用谷歌专利向量，相似度计算代码将类似：

# 谷歌专利向量的使用方法（示例）
# 假设已从 Google Patents API 获取向量
patent_vectors <- read_csv("google_patent_vectors.csv")

# 计算余弦相似度
similarity <- cosine(patent_vectors[i, ], patent_vectors[j, ])

# 对比：TF-IDF 方法
similarity <- as.numeric(dtm_norm[i, ] %*% dtm_norm[j, ])

实证研究建议：

• 初步筛查/大规模计算：使用 TF-IDF（本实现方法）
• 精细分析/小样本研究：可尝试获取谷歌专利向量进行对比验证
• 方法创新方向：训练基于中文专利全文的专用词向量模型

注意：原论文使用谷歌专利向量是为了验证方法的有效性。如果研究者有条件获取中国专利的谷歌专利向量，可以进行对比分析，进一步增强研究的可信度。

指标来源与定义

本讲义的指标构建方法参考论文：

《创新知识溢出的测度与检验——基于机器学习生成专利文本相似度的证据》

城市间知识溢出指标定义

时间条件的经济学含义

论文特别强调：来源专利的公开时间必须早于接收专利的申请时间。

这一条件保证了知识溢出的时间逻辑：

• 创新知识先在特定城市产生并公开（publish_year）
• 其他城市的创新主体学习这些知识后，申请自己的专利（apply_year）
• 如果公开时间晚于申请时间，则不可能存在知识溢出

代码实现：

sp <- df_text %>% filter(publish_year >= t - 3, publish_year <= t - 1)
tp <- df_text %>% filter(apply_year == t)

---

完整代码实现

本节将完整展示计算城市间知识溢出指标的 R 代码，并按照功能模块逐步讲解。

第0步：加载包与创建输出文件夹

在开始计算之前，我们需要加载必要的 R 包，并创建一个用于保存中间结果的文件夹。这样做有两个好处：一是避免中间结果丢失，二是支持断点续算（如果计算中断，下次可以从断点继续）。

另外，我们还会清理可能存在的旧缓存文件，避免读取到错误的数据。

完整代码

# ============================================================================
# 城市间知识溢出指标计算 - 支持断点续算（使用 jiebaR 分词）
# 双维度保存：按年份 + 按 source_city
# 中间文件保存到 res/ 文件夹（.rds 格式），使用 fs 包索引文件
#
# 论文公式：knowledgespillover_ab = Σᵢ Σⱼ sim(pᵢ,qⱼ) / (n_a × n_b)
# ============================================================================

library(tidyverse)
library(tidytext)
library(Matrix)
library(readr)
library(fs)
library(jiebaR)  # 添加 jiebaR 包用于中文分词

# ============================================================================
# 第0步：创建 res 文件夹，清理旧缓存
# ============================================================================
if (!dir.exists("res")) {
  dir.create("res")
  cat("===== 已创建 res 文件夹 =====\n\n")
}

# 清理工作目录中的旧缓存文件（如果存在）
old_cache_files <- c("dtm_norm.rds", "df_text_idx.rds")
for (f in old_cache_files) {
  if (file.exists(f)) {
    file.remove(f)
    cat(sprintf("  已删除旧缓存文件：%s\n", f))
  }
}

首先加载需要的 R 包：tidyverse 用于数据清洗，tidytext 用于文本挖掘，Matrix 用于稀疏矩阵计算，readr 用于高效读写数据，fs 用于文件系统操作。特别要提到的是 jiebaR 包，这是我们用来进行中文分词的工具。

接着创建 res/ 文件夹，所有中间结果都会保存在这里（.rds 格式），这样既便于断点续算，也能避免中间文件混乱。

最后，我们清理工作目录中可能存在的旧缓存文件（dtm_norm.rds 和 df_text_idx.rds），这些文件如果是之前计算遗留的，可能会导致错误。

---

第一步：读取数据

这一步很简单，就是读取预处理好的专利数据文件 df_text_sim.rds。这个文件应该包含以下字段：专利唯一标识（newipzlid）、申请年份（apply_year）、公开年份（publish_year）、申请人所在城市（city），以及专利摘要（摘要）。

读取后我们简单查看一下数据的维度、列名和前几行，确认数据格式符合预期。

完整代码

# ============================================================================
# 第一步：读取数据
# ============================================================================
cat("\n===== 第一步：读取 df_text_sim.rds =====\n")

df_text <- readr::read_rds("df_text_sim.rds")
cat(sprintf("  数据：%d 行 × %d 列\n", nrow(df_text), ncol(df_text)))
cat("  列名：", paste(names(df_text), collapse = ", "), "\n")

# 查看数据前几行
df_text

使用 readr::read_rds() 读取 RDS 格式的二进制文件，这种方法比 base R 的 readRDS() 更快。读取后打印数据维度和列名，确认数据加载成功，然后显示数据前几行，检查数据内容是否正确。

---

第二步：TF-IDF 计算与稀疏矩阵构建（含缓存机制）

这一步是整个流程的核心，我们需要完成以下任务：对专利摘要进行分词，计算每个词的 TF-IDF 值，然后构建稀疏矩阵并归一化（这样后续计算余弦相似度就可以通过简单的矩阵乘法实现）。

考虑到 TF-IDF 计算非常耗时，我们设计了缓存机制：计算结果会保存到 res/ 文件夹中（dtm_norm.rds 和 df_text_idx.rds）。下次运行时，如果缓存文件存在且有效，就直接加载，跳过重复计算。

完整代码

此处代码需下载讲义材料查看～

---

第三步：断点续算检查

在开始大规模计算之前，我们先检查一下 res/ 文件夹中是否已经存在部分结果。这样可以跳过已完成的部分，从断点继续计算。

这里使用了双维度保存策略：按年份保存（spillover_year_YYYY.rds）和按来源城市+年份保存（spillover_source_城市_year_YYYY.rds）。

完整代码

# ============================================================================
# 第三步：断点续算 - 检测已完成的任务（双维度）
# ============================================================================
cat("\n===== 第三步：断点续算检查 =====\n")

# 维度1：按年份（使用 .rds 格式）
year_files <- fs::dir_ls("res", regexp = "spillover_year_\\d+\\.rds$", type = "file")
years_done <- gsub(".*spillover_year_|\\.rds$", "", year_files) %>% as.integer()

if (length(years_done) > 0) {
  cat(sprintf("  已完成的年份：%s\n",
              paste(sort(years_done), collapse = ", ")))
} else {
  cat("  无已完成年份\n")
}

# 维度2：按 source_city（使用 .rds 格式）
source_files <- fs::dir_ls("res", regexp = "spillover_source_.*\\.rds$", type = "file")
sources_done <- gsub(".*spillover_source_|\\.rds$", "", source_files)

if (length(sources_done) > 0) {
  cat(sprintf("  已完成的来源城市文件：%d 个\n", length(sources_done)))
} else {
  cat("  无已完成的来源城市文件\n")
}

使用 fs::dir_ls() 列出 res/ 文件夹中匹配特定模式的文件，通过正则表达式匹配按年份保存的文件和按来源城市保存的文件，然后从文件名中提取年份或城市名称，最后打印已完成的任务，让用户了解计算进度。

---

第四步：主循环 - 计算城市间知识溢出指标

这是整个流程最核心的部分，我们需要计算每一年、每一个城市对的知识溢出指标。

具体逻辑是：遍历每一年，对每一年确定来源专利（t-3 到 t-1 年公开的专利）和目标专利（t 年申请的专利），然后遍历每个来源城市和目标城市组合，计算所有专利对的文本相似度，最后汇总得到城市间的知识溢出指标。

考虑到计算量巨大，我们设计了断点续算机制：如果某一年已经计算过（存在 res/spillover_year_YYYY.rds），就跳过；每计算完一年，立即保存结果。这样即使中断运行，下次也可以从断点继续。

完整代码

此处代码需下载讲义材料查看～

---

第五步：合并所有结果

所有年份计算完成后，我们需要将分散在 res/ 文件夹中的中间文件合并为一个完整的结果文件。这里使用 fs::dir_ls() 索引所有年份文件，然后逐个读取并合并。

完整代码

# ============================================================================
# 第五步：合并所有结果（使用 fs 包索引文件）
# ============================================================================
cat("\n===== 第五步：合并所有结果 =====\n")

# 按年份文件合并（使用 .rds 格式）
year_files <- fs::dir_ls("res", regexp = "spillover_year_\\d+\\.rds$", type = "file")
cat(sprintf("  找到 %d 个年份文件\n", length(year_files)))

all_results <- vector("list")
for (f in year_files) {
  df <- readr::read_rds(f)
  if (nrow(df) > 0) {
    all_results[[length(all_results) + 1]] <- df
  }
}

if (length(all_results) > 0) {
  spillover_df <- bind_rows(all_results)
  # 最终合并结果保存到工作目录（根目录），使用 CSV 格式便于查看
  readr::write_csv(spillover_df, "城市间知识溢出指标_tidytext.csv")
  cat(sprintf("  合并完成：%d 条记录已保存到 城市间知识溢出指标_tidytext.csv\n",
              nrow(spillover_df)))
} else {
  spillover_df <- tibble()
  cat("  无有效结果可合并\n")
}

使用 fs::dir_ls() 列出所有匹配的年份文件，逐个读取文件并检查是否有有效数据，然后使用 bind_rows() 合并所有数据框，最后保存为 CSV 格式（城市间知识溢出指标_tidytext.csv），便于查看和后续分析。

---

第六步：汇总与统计

所有年份计算完成后，我们需要对合并后的结果进行汇总，计算每年、每个城市对的知识溢出指标，并生成统计报告。

具体来说，我们按 year、source_city、target_city 分组，每个城市对只保留一行（因为第四步已经计算了平均相似度），然后生成年度统计：每年的城市对数量、平均知识溢出、标准差。

完整代码

# ============================================================================
# 第六步：汇总
# ============================================================================
cat("\n===== 第六步：汇总 =====\n")

if (nrow(spillover_df) > 0) {
  spillover_summary <- spillover_df %>%
    group_by(year, source_city, target_city) %>%
    summarise(
      n_source             = dplyr::first(n_source),
      n_target             = dplyr::first(n_target),
      knowledgespillover_ab = dplyr::first(knowledgespillover_ab),
      .groups = "drop"
    ) %>%
    arrange(year, source_city, target_city)

  # 汇总文件保存到工作目录
  readr::write_csv(spillover_summary, "城市间知识溢出指标_tidytext_汇总.csv")

  summary_stats <- spillover_summary %>%
    group_by(year) %>%
    summarise(
      n_pairs       = n(),
      mean_spillover = mean(knowledgespillover_ab, na.rm = TRUE),
      sd_spillover   = sd(knowledgespillover_ab, na.rm = TRUE),
      .groups = "drop"
    )

  # 统计文件保存到工作目录
  readr::write_csv(summary_stats, "知识溢出年度统计_tidytext.csv")

  cat("  汇总文件已保存到工作目录：\n")
  cat("    - 城市间知识溢出指标_tidytext.csv\n")
  cat("    - 城市间知识溢出指标_tidytext_汇总.csv\n")
  cat("    - 知识溢出年度统计_tidytext.csv\n")

  # 打印结果
  cat("\n===== 年度统计 =====\n")
  print(as.data.frame(summary_stats))

  cat("\n===== Top 10 知识溢出城市对 =====\n")
  top10 <- spillover_summary %>%
    arrange(desc(knowledgespillover_ab)) %>%
    head(10)
  print(as.data.frame(top10))
} else {
  cat("  无数据可汇总\n")
}

这段代码主要做三件事：首先按年份、来源城市、目标城市分组汇总数据；然后生成年度统计（城市对数量、平均知识溢出、标准差）；最后打印结果，包括年度统计表格和知识溢出最高的前10个城市对。

---

断点续算说明

代码末尾打印了断点续算的使用说明和指标定义，方便用户理解。

完整代码

此处代码需下载讲义材料查看～

点击这里跳转到 RStata 短书平台获取附件：名师讲堂｜使用 R 语言计算专利文本相似度及城市间知识溢出指标