如何用statnet/igraph创建包含边与孤立节点的网络？

嗨，我来帮你搞定这个网络构建的问题！你遇到的核心问题是没有把全部228个顶点纳入初始网络，以及后续加边时的顶点匹配错误。下面分别用igraph和statnet两个包给出具体解决方案：

一、使用igraph包的解决方案 #

方法1：直接用全量顶点表创建网络（最简便） #

你之前的问题是vertices参数只传了164个顶点的属性表，导致孤立节点被遗漏。只要把包含全部228个顶点属性的表传给这个参数，igraph会自动保留未出现在边表中的孤立节点：

# 假设fdi.attr.all是包含228个顶点所有属性的data.frame，其中必须有一列（比如"name"）作为顶点唯一ID，且与边表的from/to列匹配
fdi.graph.2003 <- graph_from_data_frame(
  d = fdi.edge.2003, 
  directed = TRUE, 
  vertices = fdi.attr.all  # 这里传全量顶点属性表，不是仅164个的子集
)

# 验证顶点数量，应该返回228
vcount(fdi.graph.2003)

方法2：先初始化全顶点空网络，再添加边（适合分步操作） #

如果之前分步操作报错，大概率是顶点ID的格式不匹配（比如边表是字符ID，空网络用的是数值索引）。按以下步骤操作可以避免：

• 从全量顶点属性表创建无任何边的空网络：

fdi.graph.empty <- graph_from_data_frame(
  d = NULL, 
  directed = TRUE, 
  vertices = fdi.attr.all
)

• 将边表的from/to转成igraph可识别的顶点ID向量（注意要保证from/to与顶点表的name列完全匹配）：

edge_pairs <- as.vector(t(fdi.edge.2003[, c("from", "to")]))

• 添加边及边属性：

fdi.graph.2003 <- add_edges(
  fdi.graph.empty, 
  edges = edge_pairs,
  # 传递边表中除from/to外的所有属性
  attr = fdi.edge.2003[, setdiff(colnames(fdi.edge.2003), c("from", "to"))]
)

关键注意点：如果边表的from/to是数值类型，顶点表的name列也要是数值类型（或者统一转成字符），否则会出现“Illegal vertex reference”这类匹配错误。

二、使用statnet包（network包）的解决方案 #

步骤1：创建包含全量顶点的空网络并添加属性 #

library(network)

# 创建有向空网络，指定顶点数量为228
fdi.net.2003 <- network.initialize(n = 228, directed = TRUE)

• 添加顶点ID属性：确保fdi.attr.all包含所有228个顶点的属性，且有唯一ID列（比如"id"）

set.vertex.attribute(fdi.net.2003, "id", fdi.attr.all$id)

• 批量添加其他顶点属性（比如国家名称、所属区域等）：

for (attr_col in setdiff(colnames(fdi.attr.all), "id")) {
  set.vertex.attribute(fdi.net.2003, attr_col, fdi.attr.all[[attr_col]])
}

步骤2：添加边及边属性 #

statnet用顶点索引（从1开始）来识别顶点，所以需要先把边表的ID转换成对应的索引：

• 创建顶点ID到网络索引的映射：

from_idx <- match(fdi.edge.2003$from, fdi.attr.all$id)
to_idx <- match(fdi.edge.2003$to, fdi.attr.all$id)

• 检查映射是否有NA（如果有，说明边表中有不存在的顶点ID，需要先清理数据）：

if (any(is.na(from_idx)) || any(is.na(to_idx))) {
  stop("边表中存在未在全量顶点表中出现的ID，请检查数据一致性！")
}

• 添加边：

add.edges(fdi.net.2003, tail = from_idx, head = to_idx)

• 添加边属性（比如权重）：

if ("weight" %in% colnames(fdi.edge.2003)) {
  set.edge.attribute(fdi.net.2003, "weight", fdi.edge.2003$weight)
}

核心总结 #

• 不管用哪个包，核心原则都是：先确保网络包含所有228个顶点，再添加边。
• 顶点ID的格式一致性是关键，这是导致“Illegal vertex reference”报错的最常见原因，一定要保证边表中的起点/终点ID在全量顶点表中都存在，且类型统一。

内容的提问来源于stack exchange，提问作者Filippo Santi