在R语言ABM中构建指定规模的对称Agent社交网络矩阵
实现方法与R包推荐
推荐R包:igraph
igraph是R中专注于网络分析的工具包,内置多种生成指定度序列无向图的函数,能完美满足你的需求——生成对称、无自环且每个Agent连接数在1-10之间的邻接矩阵。
具体实现步骤
- 安装并加载包
install.packages("igraph") library(igraph)
- 生成符合要求的度序列
先随机生成100个1-10之间的整数作为每个Agent的目标连接数,注意无向图要求所有节点度的总和为偶数,若总和为奇数则微调其中一个节点的度(保持在1-10范围内):
set.seed(123) # 可选,确保结果可复现 degrees <- sample(1:10, 100, replace = TRUE) # 确保度总和为偶数 while(sum(degrees) %% 2 != 0) { # 随机选一个未达最大度的节点,将其度加1 idx <- sample(which(degrees != 10), 1) degrees[idx] <- degrees[idx] + 1 }
- 生成无向图并转换为邻接矩阵
使用sample_degseq函数生成符合度序列的简单无向图(无自环、无多重边),再转换为对称邻接矩阵:
# 生成无向图,method="vl"是Viger-Latouche算法,专门生成简单图 g <- sample_degseq(degrees, method = "vl") # 转换为稠密邻接矩阵 adj_matrix <- as_adjacency_matrix(g, type = "both", sparse = FALSE)
- 验证结果
通过行和检查每个Agent的连接数是否符合要求:
rowSums(adj_matrix) # 结果应与degrees一致
为什么不用你原来的代码?
你原来的随机生成方法无法控制每个节点的连接数,而通过度序列生成的方式能精准约束每个Agent的社交网络规模,同时保证对称性和无自环的约束条件。
内容的提问来源于stack exchange,提问作者clara
相关产品推荐
相关产品推荐

