简介这份针对2018年国际贸易网络分布分析的R语言资源面向经济学、社会学及复杂网络研究者解决从原始贸易数据到网络指标计算与模型解读的完整流程落地问题。资源包内仅有1个R脚本却集中覆盖网络密度、平均路径长度、传递性、互易性、分类性、自旋玻璃社区检测、结构等效性以及指数随机图模型ERGM等多类方法实现压缩包大小仅7KB轻量易用。脚本按分析模块清晰组织可对照理论定义分段运行并拆解输出既能直观理解密度、路径长度等基础指标在全球贸易格局中的含义也可借助社区检测识别区域贸易板块再通过ERGM探究贸易关系形成的影响机制。当前已有233人学习适合具备一定R基础并希望结合真实贸易数据进阶网络分析的中高级用户。 2018年全球国家间贸易数据用R跑一遍网络分析。这篇博文围绕TradeNetworkDistributionsAnalysis项目把进出口网络的密度、平均路径长度、传递性、互易性、分类性、自旋玻璃社区检测、结构等效性和ERGM模型的完整流程写出来。适合正在学网络分析的R用户以及想从网络视角看国际贸易的研究者。我会按实际操作的顺序把每一步的思路、R代码和踩过的坑都摊开讲。1. 项目整体设计与网络构建1.1 为什么把贸易关系看成网络国际贸易数据天然是一张有向加权网络每个国家或地区是节点国家之间的出口和进口流是边。传统统计只能回答“谁和谁贸易多”而网络分析能回答更深一层的问题这个贸易体系是紧密还是松散是否存在核心—边缘结构贸易关系是如何形成的对于2018年全球贸易可以用公开贸易流数据比如UN Comtrade或世界银行WITS构建邻接矩阵行是出口国列是进口国单元格是贸易额。这里有一个关键选择网络必须是有向的因为A国对B国出口与B国对A国出口是两条不同方向的边只有同时记录方向才能正确计算互易性、入度/出度等指标。同时贸易额是权重边权重取出口额还是进口额取决于你关注哪个视角。我做项目时直接采用出口矩阵因为各国上报的出口数据通常比进口数据更准确。2018年是个不错的分析窗口全球贸易总量处于周期波动中区域贸易协定密集网络结构相对稳定又有足够的公开数据支撑。分析结果既能反映当年的贸易格局也能作为后续年份对比的基线。1.2 R语言中构建网络对象在R里最常用的网络分析包是igraph和statnet。我的做法是先用igraph做指标计算和社区检测再用statnet里的network对象做ERGM。igraph建网很简单如果手头有edge list三行代码就能得到网络对象library(igraph) edges - read.csv(trade_edges_2018.csv, header TRUE) g - graph_from_data_frame(edges, directed TRUE)这里的edges至少包含三列from、to、weight分别代表出口国、进口国和贸易额。from和to最好先转成字符型否则graph_from_data_frame会做隐式转换可能导致后续匹配节点属性时出错。建完网络后我习惯立刻检查一下基本规模summary(g)这个输出会列出节点数、边数、是否有向、是否有权重。节点数一般在150到200之间取决于你纳入多少国家。如果某些小经济体数据缺失严重建议直接过滤掉避免网络被散点拖累。网络构建是后面所有分析的地基最容易掉坑的是重复边和缺失值。如果同一个出口国-进口国组合出现多次必须提前聚合否则网络会包含平行边。igraph中有simplify()可以合并但更稳妥的是在数据处理阶段用dplyr聚合library(dplyr) edges_agg - edges %% group_by(from, to) %% summarise(weight sum(weight, na.rm TRUE), .groups drop)缺失值要么删掉要么用0填充。0表示没有贸易关系而NA会导致后续计算直接报错。另外有些数据源会把“地区”或“未分类”行也包含进来这些行不属于真正的国家节点必须提前剔除。还有一个特别隐蔽的坑自环。有些出口数据会包含国家对自身的经济体报告如果不加过滤graph_from_data_frame会创建自环自环会影响密度、路径长度、传递性等多个指标。处理方法是直接过滤edges_clean - edges_agg[edges_agg$from ! edges_agg$to, ]2. 基础网络指标密度、路径、传递与互易2.1 网络密度与平均路径长度网络密度是所有实际存在的边数除以最大可能边数。在有向网络中最大可能边数是n*(n-1)因为贸易网络一般不考虑国家内部的自我贸易。R中直接用edge_density(g)计算。密度越接近1说明国家之间的贸易联系越紧密。2018年全球贸易网络的密度大概在0.4~0.6之间具体取决于纳入的国家数量和数据过滤阈值。如果网络太稀疏后面的社区检测和ERGM都会不稳定。平均路径长度average path length是衡量网络“小世界”性质的核心指标它表示任意两个可达节点之间的最短路径的平均步数。在R里mean_distance(g)对于有向图如果存在不可达节点对mean_distance默认会忽略这些对并且会在结果中带一个告警。你可以显式设置unconnectedTRUE来确认行为。很多人在这一步会忽略一个重要问题贸易网络中存在少数孤立或只出不进的国家这些国家会拉低网络连通性。建议先找出最大连通分量在核心网络上计算路径指标否则平均路径长度会被严重低估或计算失败。贸易网络的路径长度通常很短2到4步之间。这说明国际贸易高度全球化A国通过B国间接与C国贸易是很常见的。路径长度还可以结合最短路径的中介性betweenness来识别“中间人国家”这类分析对供应链研究特别有用。2.2 传递性与互易性的业务含义传递性transitivity也叫聚类系数衡量的是“朋友的朋友也是朋友”的程度。在贸易网络中如果A国大量出口到B国B国又大量出口到C国那么A国和C国之间也倾向于有贸易往来这样会形成三角形关系。R中transitivity(g)返回全局传递系数transitivity(g, typelocal)返回每个节点的局部聚类系数。贸易网络的传递性通常较高说明贸易关系存在区域集聚效应例如地理邻近或共同协定国家之间更容易形成紧密的三角贸易圈。互易性reciprocity是有向网络特有的指标它表示双向边的比例网络中既存在A→B也存在B→A的边占所有有向边的比例。R中reciprocity(g)默认计算“互易边数/总边数”的比率。贸易网络的互易性往往很高因为大多数国家之间都是既出口又进口只是贸易额不同。如果使用二值化网络只看边是否存在互易性可能超过0.7如果使用加权网络定义互易性会更复杂需要比较两个方向的权重。igraph里的reciprocity()只支持无权网络如果你需要加权互易性建议自己写一个简单函数计算两个方向权重的镜像程度。我常用的做法是# 加权互易性双向一致边权重占总权重比例 w - as_adjacency_matrix(g, attr weight, sparse FALSE) recip_weight - sum(pmin(w, t(w))) / sum(w)这种加权指标可以更真实地反映贸易的对称性而非仅仅看“有没有”双向关系。2.3 分类性核心—边缘结构的直接证据分类性assortativity也叫同配性衡量的是“喜欢和相似的人在一起”的程度。在网络分析里通常用度—度同配系数判断高贸易额国家是否倾向于连接高贸易额国家。R中assortativity_degree(g)取值范围是[-1,1]正值表示同配负值表示异配。做2018年贸易网络时我发现二值网络的分类性通常是正的说明少数贸易大国之间形成紧密的核心—核心连接而小国主要连接大国形成一个典型的“核心—边缘”结构。不过分类性只依赖度即连接数量而贸易网络更重要的是边的权重。如果你用强度strength替代度即把每个国家的总出口额和总进口额作为“重要性”指标也许会得到不一样的结果。我建议同时计算两个版本一个是基于度的分类性另一个是基于强度的分类性。实际项目中强度分类性可能比度分类性更能揭示资本和商品流向的集中趋势。3. 社区检测与结构等效性3.1 自旋玻璃社区检测算法原理社区检测是网络分析的重头戏目标是找出网络中联系紧密的国家组团。我选了自旋玻璃Spin Glass算法理由是它基于Potts自旋模型能够处理带权重的网络并且通过模拟退火找到全局最优的社区划分。在igraph里sg - cluster_spinglass(g, weights E(g)$weight, spins 8)运行后用membership(sg)取出每个国家的社区标签。这个算法的一个明显优势是允许节点带权重可以充分利用贸易额信息。但它有几个非常实际的前提条件第一算法要求网络是连通的所以要先删掉孤立节点或者只在最大的连通分量上运行。第二spins参数指定允许的“自旋”数也就是你预期组团的个数。设置得过多会导致社区过度分裂过少又会把不同社区强行合并。我的经验是先设定spins8再根据结果调整。如果网络规模较大还可以调大start.temp让模拟退火有更充分的搜索过程。跑完以后用V(g)$community - membership(sg)把社区信息保存在网络对象里。2018年的贸易网络通常能分出一批高度抱团的区域经济集团但也可能混着一些全球性贸易大国因为它们和很多区域都有联系社区归属会比较模糊。这时候不要怀疑算法反而应该把这看作全球贸易“多重区域化”的证据。另外自旋玻璃算法结果受随机种子影响比较大最好设置set.seed()固定种子并多跑几次看社区划分的稳定性。如果两次运行结果差异很大说明网络本身没有清晰的社区结构这时候即使有输出也不能过度解读。3.2 结构等效性分析结构等效性Structural Equivalence度量的是两个节点在网络中扮演的角色是否相同。简单来说如果两个国家的贸易伙伴集合高度重叠即使它们之间没有直接贸易也可以认为它们在网络中的“位置”等效。计算结构等效性通常有两种方法一是计算节点邻接向量的欧氏距离或余弦相似度二是用sna包中的sedegree或correlation距离。igraph中自带similarity()函数搭配methodjaccard可以计算基于共同邻居的相似度。但要注意这个函数默认把网络当成无权网络且只计算有边连接的相似性对于没有共同邻居的节点会直接给0。更严谨的做法是自建一个矩阵把每个节点的出边和入边分别提取出来组成一个2k维向量然后计算两两之间的欧氏距离adj - as_adjacency_matrix(g, sparse FALSE) node_vec - cbind(adj, t(adj)) # 每个节点对应一个2n维向量 dist_mat - as.matrix(dist(node_vec, method euclidean))距离越小说明两个国家在整个网络中的结构位置越接近。比如一个内陆小国和一个同体量的岛国可能彼此没有直接大量贸易但它们都主要依赖某个大国进行进出口那么它们的结构等效性就很高在供应链中扮演的角色类似。实际操作中我还会把结构等效性和社区检测结果放在一起画热图能直观看到哪些国家在贸易网络中扮演相似的“桥接型”或“边缘型”角色。这种描述性分析能为后续ERGM的节点属性选择提供灵感比如如果你发现某些国家结构等效可能意味着存在相似的经济体量或区域特色。4. 指数随机图模型ERGM4.1 ERGM到底在拟合什么前面算的密度、互易、传递、分类性都是描述统计只能告诉我们网络“是什么样”而指数随机图模型ERGM解决的则是“为什么网络会呈现这个样子”。ERGM把网络视为随机变量通过一组网络统计量边数、互易边数、三角形数量、度数分布等来估计网络形成的概率分布。你可以把它类比成逻辑回归在网络中每一条可能的边都有一个出现概率而这个概率由各种局部结构项共同决定。对于贸易网络最值得考察的生成机制有三类互易机制国家之间倾向于相互开放市场而不是单向贸易。传递机制贸易关系容易形成三角闭合也就是“伙伴的伙伴也是伙伴”。偏好依附机制贸易大国更容易吸引新的贸易连接形成核心—边缘结构。ERGM可以同时估计这些机制的相对重要性。因为ERGM标准形式处理的是二元网络所以通常要把加权网络二值化。比如设定一个贸易额阈值超过阈值就视为存在贸易关系。阈值的选择很关键我一般取所有贸易额的中位数或75%分位数这样既能保留核心贸易联系又不会让网络过于稠密。如果阈值太低网络密度接近1ERGM拟合会出现退化如果阈值太高网络过于稀疏很多统计量趋近0估计不准确。4.2 R语言拟合与结果解读R里做ERGM主要用到statnet套件。刚才我们用igraph建了网络现在需要转换成statnet的network对象。我会用intergraph包完成转换library(statnet) library(intergraph) net_bin - asNetwork(g_bin) # g_bin是二值化后的igraph对象 model - ergm(net_bin ~ edges mutual gwesp(0.25) gwidegree(0.3, fixedTRUE) gwodegree(0.3, fixedTRUE)) summary(model)这里edges对应网络密度mutual对应互易性gwesp是加权共享伙伴项捕捉传递性趋势gwidegree和gwodegree分别捕捉入度和出度的中心化趋势。固定衰减参数是为了让模型更容易收敛我通常用0.25到0.5之间的值。拟合完成后看Estimate列和Pr(|z|)列。Estimate是对数比值比所以exp(Estimate)才是解释为优势比。比如mutual的Estimate如果是1.8exp(1.8)≈6.05说明在控制其他结构的情况下存在反向贸易连接的概率是随机网络下的6倍。这在贸易网络里几乎总是显著的。gwesp的正系数同样说明贸易关系容易形成闭合三角。如果gwidegree为负则说明高入度国家的聚集程度低于随机网络可能存在“少数国家接收多数连接”的星型结构。ERGM结果解释要非常小心因为多个统计项之间存在共线性。我通常跑两三个简化模型逐步加入互易项、传递项、度数项然后比较AIC/BIC。模型收敛失败时可以增加MCMC样本数比如model2 - ergm(net_bin ~ edges mutual, control control.ergm(MCMC.samplesize 10000))还可以调整衰变参数。如果仍不收敛检查网络密度是否过高或者某个统计项是否常数。记住ERGM不是黑箱每加一个项都要有业务上的理由不然就是过度拟合。5. 常见问题与实操心得5.1 数据清洗与网络对象构建的坑我踩过最大的坑是自环前面已经提过。另一个坑是节点编码不一致同一国家的ISO代码可能在不同文件中格式不同导致网络被拆成多个孤立节点。强烈建议在数据合并后做一次快速验证比如检查边数量nrow(unique(edges_clean[, c(from, to)]))这个数应该等于网络中的边数如果小于summary(g)中的边数说明存在重复边没清理干净。如果大于说明边表有冗余。还有一个容易忽略的问题节点属性比如GDP、地区分类和网络节点名的匹配顺序。igraph中V(g)$name的顺序并不总是和数据框的顺序一致要用match()或join的方式合并千万不要用cbind。我就吃过这种亏导致社区颜色和实际节点错位画出来的图误导了自己好几天。5.2 算法参数调整的经验自旋玻璃算法的结果对spins参数异常敏感建议多跑几次取稳定解。我在2018年贸易网络上跑的时候spins8和spins6的结果差异很大最后我结合经济常识选择了一个与全球经济板块分布较为吻合的划分。另外如果网络很大先用简化后的核心网络比如剔除贸易额低于100万美元的小边可以减少计算负担。ERGM更是重量级100个节点左右的网络都可能跑很久。我一般在跑之前先用描述性指标看看数据如果网络太密密度0.6ERGM拟合数值可能不稳定需要先稀疏化比如删掉权重最低的20%的边。对于大网络还可以考虑使用ergm包的多线程版本或采样MCMC。我在实际项目中甚至把网络按地区拆开分别建模反而更有解释力。5.3 结果可视化小技巧网络可视化对代码要求不高但能极大提升分析质量。我通常把社区检测结果映射到节点颜色节点大小与总贸易额成正比边的透明度与贸易额成正比。可以用igraph的plot()快速看整体也可以用ggraph包做更精致的图。特别提醒边数很多时直接画会变成一团毛线最好先过滤掉权重低的边只保留Top 20%的边。这样既保持了核心结构又让视觉清晰。对于ERGM的结果可以用gof()函数做拟合优度诊断画出模型预测的度数分布、边共享伙伴分布与实际网络的对比图。如果拟合效果好模型模拟出的网络应当能复现出实际网络的结构特征。6. 最后的补充加权网络带来的额外发现在TradeNetworkDistributionsAnalysis项目中我最惊喜的不是高互易性或高传递性而是分类性在二值网络和加权网络之间的差异。二值网络显示正同配但一旦用强度加权度重新计算同配系数可能明显下降甚至反转。这说明贸易强国的连接强度远比连接数量重要。如果你也想分析贸易网络建议不要只看二值化的度也试试用strength加权度重新计算全部指标往往能发现不一样的故事。另一个经验是网络分析真正考验人的不是算法运行而是每一步选择的业务解释。密度高不代表贸易体系稳定互易性高也不代表贸易公平。指标永远只是工具最终要回到经济现实中去验证。希望这篇博文能帮你少走一点弯路把R网络分析真正用到自己的数据上。本文还有配套的精品资源点击获取 
长河君 · 2026/9/2 2:16:00 · 阅读 1.2k R语言网络分析实战：从2018全球贸易数据到ERGM建模
简介这份针对2018年国际贸易网络分布分析的R语言资源面向经济学、社会学及复杂网络研究者解决从原始贸易数据到网络指标计算与模型解读的完整流程落地问题。资源包内仅有1个R脚本却集中覆盖网络密度、平均路径长度、传递性、互易性、分类性、自旋玻璃社区检测、结构等效性以及指数随机图模型ERGM等多类方法实现压缩包大小仅7KB轻量易用。脚本按分析模块清晰组织可对照理论定义分段运行并拆解输出既能直观理解密度、路径长度等基础指标在全球贸易格局中的含义也可借助社区检测识别区域贸易板块再通过ERGM探究贸易关系形成的影响机制。当前已有233人学习适合具备一定R基础并希望结合真实贸易数据进阶网络分析的中高级用户。 2018年全球国家间贸易数据用R跑一遍网络分析。这篇博文围绕TradeNetworkDistributionsAnalysis项目把进出口网络的密度、平均路径长度、传递性、互易性、分类性、自旋玻璃社区检测、结构等效性和ERGM模型的完整流程写出来。适合正在学网络分析的R用户以及想从网络视角看国际贸易的研究者。我会按实际操作的顺序把每一步的思路、R代码和踩过的坑都摊开讲。1. 项目整体设计与网络构建1.1 为什么把贸易关系看成网络国际贸易数据天然是一张有向加权网络每个国家或地区是节点国家之间的出口和进口流是边。传统统计只能回答“谁和谁贸易多”而网络分析能回答更深一层的问题这个贸易体系是紧密还是松散是否存在核心—边缘结构贸易关系是如何形成的对于2018年全球贸易可以用公开贸易流数据比如UN Comtrade或世界银行WITS构建邻接矩阵行是出口国列是进口国单元格是贸易额。这里有一个关键选择网络必须是有向的因为A国对B国出口与B国对A国出口是两条不同方向的边只有同时记录方向才能正确计算互易性、入度/出度等指标。同时贸易额是权重边权重取出口额还是进口额取决于你关注哪个视角。我做项目时直接采用出口矩阵因为各国上报的出口数据通常比进口数据更准确。2018年是个不错的分析窗口全球贸易总量处于周期波动中区域贸易协定密集网络结构相对稳定又有足够的公开数据支撑。分析结果既能反映当年的贸易格局也能作为后续年份对比的基线。1.2 R语言中构建网络对象在R里最常用的网络分析包是igraph和statnet。我的做法是先用igraph做指标计算和社区检测再用statnet里的network对象做ERGM。igraph建网很简单如果手头有edge list三行代码就能得到网络对象library(igraph) edges - read.csv(trade_edges_2018.csv, header TRUE) g - graph_from_data_frame(edges, directed TRUE)这里的edges至少包含三列from、to、weight分别代表出口国、进口国和贸易额。from和to最好先转成字符型否则graph_from_data_frame会做隐式转换可能导致后续匹配节点属性时出错。建完网络后我习惯立刻检查一下基本规模summary(g)这个输出会列出节点数、边数、是否有向、是否有权重。节点数一般在150到200之间取决于你纳入多少国家。如果某些小经济体数据缺失严重建议直接过滤掉避免网络被散点拖累。网络构建是后面所有分析的地基最容易掉坑的是重复边和缺失值。如果同一个出口国-进口国组合出现多次必须提前聚合否则网络会包含平行边。igraph中有simplify()可以合并但更稳妥的是在数据处理阶段用dplyr聚合library(dplyr) edges_agg - edges %% group_by(from, to) %% summarise(weight sum(weight, na.rm TRUE), .groups drop)缺失值要么删掉要么用0填充。0表示没有贸易关系而NA会导致后续计算直接报错。另外有些数据源会把“地区”或“未分类”行也包含进来这些行不属于真正的国家节点必须提前剔除。还有一个特别隐蔽的坑自环。有些出口数据会包含国家对自身的经济体报告如果不加过滤graph_from_data_frame会创建自环自环会影响密度、路径长度、传递性等多个指标。处理方法是直接过滤edges_clean - edges_agg[edges_agg$from ! edges_agg$to, ]2. 基础网络指标密度、路径、传递与互易2.1 网络密度与平均路径长度网络密度是所有实际存在的边数除以最大可能边数。在有向网络中最大可能边数是n*(n-1)因为贸易网络一般不考虑国家内部的自我贸易。R中直接用edge_density(g)计算。密度越接近1说明国家之间的贸易联系越紧密。2018年全球贸易网络的密度大概在0.4~0.6之间具体取决于纳入的国家数量和数据过滤阈值。如果网络太稀疏后面的社区检测和ERGM都会不稳定。平均路径长度average path length是衡量网络“小世界”性质的核心指标它表示任意两个可达节点之间的最短路径的平均步数。在R里mean_distance(g)对于有向图如果存在不可达节点对mean_distance默认会忽略这些对并且会在结果中带一个告警。你可以显式设置unconnectedTRUE来确认行为。很多人在这一步会忽略一个重要问题贸易网络中存在少数孤立或只出不进的国家这些国家会拉低网络连通性。建议先找出最大连通分量在核心网络上计算路径指标否则平均路径长度会被严重低估或计算失败。贸易网络的路径长度通常很短2到4步之间。这说明国际贸易高度全球化A国通过B国间接与C国贸易是很常见的。路径长度还可以结合最短路径的中介性betweenness来识别“中间人国家”这类分析对供应链研究特别有用。2.2 传递性与互易性的业务含义传递性transitivity也叫聚类系数衡量的是“朋友的朋友也是朋友”的程度。在贸易网络中如果A国大量出口到B国B国又大量出口到C国那么A国和C国之间也倾向于有贸易往来这样会形成三角形关系。R中transitivity(g)返回全局传递系数transitivity(g, typelocal)返回每个节点的局部聚类系数。贸易网络的传递性通常较高说明贸易关系存在区域集聚效应例如地理邻近或共同协定国家之间更容易形成紧密的三角贸易圈。互易性reciprocity是有向网络特有的指标它表示双向边的比例网络中既存在A→B也存在B→A的边占所有有向边的比例。R中reciprocity(g)默认计算“互易边数/总边数”的比率。贸易网络的互易性往往很高因为大多数国家之间都是既出口又进口只是贸易额不同。如果使用二值化网络只看边是否存在互易性可能超过0.7如果使用加权网络定义互易性会更复杂需要比较两个方向的权重。igraph里的reciprocity()只支持无权网络如果你需要加权互易性建议自己写一个简单函数计算两个方向权重的镜像程度。我常用的做法是# 加权互易性双向一致边权重占总权重比例 w - as_adjacency_matrix(g, attr weight, sparse FALSE) recip_weight - sum(pmin(w, t(w))) / sum(w)这种加权指标可以更真实地反映贸易的对称性而非仅仅看“有没有”双向关系。2.3 分类性核心—边缘结构的直接证据分类性assortativity也叫同配性衡量的是“喜欢和相似的人在一起”的程度。在网络分析里通常用度—度同配系数判断高贸易额国家是否倾向于连接高贸易额国家。R中assortativity_degree(g)取值范围是[-1,1]正值表示同配负值表示异配。做2018年贸易网络时我发现二值网络的分类性通常是正的说明少数贸易大国之间形成紧密的核心—核心连接而小国主要连接大国形成一个典型的“核心—边缘”结构。不过分类性只依赖度即连接数量而贸易网络更重要的是边的权重。如果你用强度strength替代度即把每个国家的总出口额和总进口额作为“重要性”指标也许会得到不一样的结果。我建议同时计算两个版本一个是基于度的分类性另一个是基于强度的分类性。实际项目中强度分类性可能比度分类性更能揭示资本和商品流向的集中趋势。3. 社区检测与结构等效性3.1 自旋玻璃社区检测算法原理社区检测是网络分析的重头戏目标是找出网络中联系紧密的国家组团。我选了自旋玻璃Spin Glass算法理由是它基于Potts自旋模型能够处理带权重的网络并且通过模拟退火找到全局最优的社区划分。在igraph里sg - cluster_spinglass(g, weights E(g)$weight, spins 8)运行后用membership(sg)取出每个国家的社区标签。这个算法的一个明显优势是允许节点带权重可以充分利用贸易额信息。但它有几个非常实际的前提条件第一算法要求网络是连通的所以要先删掉孤立节点或者只在最大的连通分量上运行。第二spins参数指定允许的“自旋”数也就是你预期组团的个数。设置得过多会导致社区过度分裂过少又会把不同社区强行合并。我的经验是先设定spins8再根据结果调整。如果网络规模较大还可以调大start.temp让模拟退火有更充分的搜索过程。跑完以后用V(g)$community - membership(sg)把社区信息保存在网络对象里。2018年的贸易网络通常能分出一批高度抱团的区域经济集团但也可能混着一些全球性贸易大国因为它们和很多区域都有联系社区归属会比较模糊。这时候不要怀疑算法反而应该把这看作全球贸易“多重区域化”的证据。另外自旋玻璃算法结果受随机种子影响比较大最好设置set.seed()固定种子并多跑几次看社区划分的稳定性。如果两次运行结果差异很大说明网络本身没有清晰的社区结构这时候即使有输出也不能过度解读。3.2 结构等效性分析结构等效性Structural Equivalence度量的是两个节点在网络中扮演的角色是否相同。简单来说如果两个国家的贸易伙伴集合高度重叠即使它们之间没有直接贸易也可以认为它们在网络中的“位置”等效。计算结构等效性通常有两种方法一是计算节点邻接向量的欧氏距离或余弦相似度二是用sna包中的sedegree或correlation距离。igraph中自带similarity()函数搭配methodjaccard可以计算基于共同邻居的相似度。但要注意这个函数默认把网络当成无权网络且只计算有边连接的相似性对于没有共同邻居的节点会直接给0。更严谨的做法是自建一个矩阵把每个节点的出边和入边分别提取出来组成一个2k维向量然后计算两两之间的欧氏距离adj - as_adjacency_matrix(g, sparse FALSE) node_vec - cbind(adj, t(adj)) # 每个节点对应一个2n维向量 dist_mat - as.matrix(dist(node_vec, method euclidean))距离越小说明两个国家在整个网络中的结构位置越接近。比如一个内陆小国和一个同体量的岛国可能彼此没有直接大量贸易但它们都主要依赖某个大国进行进出口那么它们的结构等效性就很高在供应链中扮演的角色类似。实际操作中我还会把结构等效性和社区检测结果放在一起画热图能直观看到哪些国家在贸易网络中扮演相似的“桥接型”或“边缘型”角色。这种描述性分析能为后续ERGM的节点属性选择提供灵感比如如果你发现某些国家结构等效可能意味着存在相似的经济体量或区域特色。4. 指数随机图模型ERGM4.1 ERGM到底在拟合什么前面算的密度、互易、传递、分类性都是描述统计只能告诉我们网络“是什么样”而指数随机图模型ERGM解决的则是“为什么网络会呈现这个样子”。ERGM把网络视为随机变量通过一组网络统计量边数、互易边数、三角形数量、度数分布等来估计网络形成的概率分布。你可以把它类比成逻辑回归在网络中每一条可能的边都有一个出现概率而这个概率由各种局部结构项共同决定。对于贸易网络最值得考察的生成机制有三类互易机制国家之间倾向于相互开放市场而不是单向贸易。传递机制贸易关系容易形成三角闭合也就是“伙伴的伙伴也是伙伴”。偏好依附机制贸易大国更容易吸引新的贸易连接形成核心—边缘结构。ERGM可以同时估计这些机制的相对重要性。因为ERGM标准形式处理的是二元网络所以通常要把加权网络二值化。比如设定一个贸易额阈值超过阈值就视为存在贸易关系。阈值的选择很关键我一般取所有贸易额的中位数或75%分位数这样既能保留核心贸易联系又不会让网络过于稠密。如果阈值太低网络密度接近1ERGM拟合会出现退化如果阈值太高网络过于稀疏很多统计量趋近0估计不准确。4.2 R语言拟合与结果解读R里做ERGM主要用到statnet套件。刚才我们用igraph建了网络现在需要转换成statnet的network对象。我会用intergraph包完成转换library(statnet) library(intergraph) net_bin - asNetwork(g_bin) # g_bin是二值化后的igraph对象 model - ergm(net_bin ~ edges mutual gwesp(0.25) gwidegree(0.3, fixedTRUE) gwodegree(0.3, fixedTRUE)) summary(model)这里edges对应网络密度mutual对应互易性gwesp是加权共享伙伴项捕捉传递性趋势gwidegree和gwodegree分别捕捉入度和出度的中心化趋势。固定衰减参数是为了让模型更容易收敛我通常用0.25到0.5之间的值。拟合完成后看Estimate列和Pr(|z|)列。Estimate是对数比值比所以exp(Estimate)才是解释为优势比。比如mutual的Estimate如果是1.8exp(1.8)≈6.05说明在控制其他结构的情况下存在反向贸易连接的概率是随机网络下的6倍。这在贸易网络里几乎总是显著的。gwesp的正系数同样说明贸易关系容易形成闭合三角。如果gwidegree为负则说明高入度国家的聚集程度低于随机网络可能存在“少数国家接收多数连接”的星型结构。ERGM结果解释要非常小心因为多个统计项之间存在共线性。我通常跑两三个简化模型逐步加入互易项、传递项、度数项然后比较AIC/BIC。模型收敛失败时可以增加MCMC样本数比如model2 - ergm(net_bin ~ edges mutual, control control.ergm(MCMC.samplesize 10000))还可以调整衰变参数。如果仍不收敛检查网络密度是否过高或者某个统计项是否常数。记住ERGM不是黑箱每加一个项都要有业务上的理由不然就是过度拟合。5. 常见问题与实操心得5.1 数据清洗与网络对象构建的坑我踩过最大的坑是自环前面已经提过。另一个坑是节点编码不一致同一国家的ISO代码可能在不同文件中格式不同导致网络被拆成多个孤立节点。强烈建议在数据合并后做一次快速验证比如检查边数量nrow(unique(edges_clean[, c(from, to)]))这个数应该等于网络中的边数如果小于summary(g)中的边数说明存在重复边没清理干净。如果大于说明边表有冗余。还有一个容易忽略的问题节点属性比如GDP、地区分类和网络节点名的匹配顺序。igraph中V(g)$name的顺序并不总是和数据框的顺序一致要用match()或join的方式合并千万不要用cbind。我就吃过这种亏导致社区颜色和实际节点错位画出来的图误导了自己好几天。5.2 算法参数调整的经验自旋玻璃算法的结果对spins参数异常敏感建议多跑几次取稳定解。我在2018年贸易网络上跑的时候spins8和spins6的结果差异很大最后我结合经济常识选择了一个与全球经济板块分布较为吻合的划分。另外如果网络很大先用简化后的核心网络比如剔除贸易额低于100万美元的小边可以减少计算负担。ERGM更是重量级100个节点左右的网络都可能跑很久。我一般在跑之前先用描述性指标看看数据如果网络太密密度0.6ERGM拟合数值可能不稳定需要先稀疏化比如删掉权重最低的20%的边。对于大网络还可以考虑使用ergm包的多线程版本或采样MCMC。我在实际项目中甚至把网络按地区拆开分别建模反而更有解释力。5.3 结果可视化小技巧网络可视化对代码要求不高但能极大提升分析质量。我通常把社区检测结果映射到节点颜色节点大小与总贸易额成正比边的透明度与贸易额成正比。可以用igraph的plot()快速看整体也可以用ggraph包做更精致的图。特别提醒边数很多时直接画会变成一团毛线最好先过滤掉权重低的边只保留Top 20%的边。这样既保持了核心结构又让视觉清晰。对于ERGM的结果可以用gof()函数做拟合优度诊断画出模型预测的度数分布、边共享伙伴分布与实际网络的对比图。如果拟合效果好模型模拟出的网络应当能复现出实际网络的结构特征。6. 最后的补充加权网络带来的额外发现在TradeNetworkDistributionsAnalysis项目中我最惊喜的不是高互易性或高传递性而是分类性在二值网络和加权网络之间的差异。二值网络显示正同配但一旦用强度加权度重新计算同配系数可能明显下降甚至反转。这说明贸易强国的连接强度远比连接数量重要。如果你也想分析贸易网络建议不要只看二值化的度也试试用strength加权度重新计算全部指标往往能发现不一样的故事。另一个经验是网络分析真正考验人的不是算法运行而是每一步选择的业务解释。密度高不代表贸易体系稳定互易性高也不代表贸易公平。指标永远只是工具最终要回到经济现实中去验证。希望这篇博文能帮你少走一点弯路把R网络分析真正用到自己的数据上。本文还有配套的精品资源点击获取