r - 如何从两个二分网络制作一个有向三分网络？

Question

对不起，因为我觉得这个疑问应该更简单，但我找不到满意的答案。

我有两个定量二分网络（显示 AB 和 BC 之间的生态关系）。我的问题是我不知道如何加入两者以建立一个有向的定量三方网络（如典型的食物网）。每个二分网络的 B 级具有相同的顶点组成（在数量和顶点名称方面）。此外，A 级和 C 级之间不能交互。因此，最终三方网络的顺序和方向应该是C->B->A

有什么建议么？

感谢您的关注！

Answer 1

这个怎么样：

准备一些数据（我希望你会为我做饭:)）

library(igraph)
set.seed(666)
# herbivore-plant
m_hp <- matrix( rbinom(12, 10, p=0.2), 4, 3)
dimnames(m_hp) <- list(
  consuming=paste0("h", seq(1, nrow(m_hp))),
  consumed=paste0("p", seq(1, ncol(m_hp)))
)
# carnivore-herbivore
m_ch <- matrix( rbinom(20, 10, p=0.2), 5, 4)
dimnames(m_ch) <- list(
  consuming=paste0("c", seq(1, nrow(m_ch))),
  consumed=paste0("h", seq(1, ncol(m_ch)))
)

...所以它看起来像你的（我想）：

m_hp

##          consumed
## consuming p1 p2 p3
##        h1  3  1  0
##        h2  1  3  1
##        h3  5  5  3
##        h4  1  2  0

m_ch

##          consumed
## consuming h1 h2 h3 h4
##        c1  0  4  0  2
##        c2  1  2  1  2
##        c3  1  2  3  3
##        c4  3  5  1  2
##        c5  0  2  0  2

现在通过边缘列表将它们变成 igraph 对象

el_hp <- as.data.frame(as.table(m_hp), stringsAsFactors = FALSE)
el_ch <- as.data.frame(as.table(m_ch), stringsAsFactors = FALSE)
el <- rbind(el_hp, el_ch)
g <- graph.data.frame( el[el$Freq != 0 , ]  )
V(g)$type <- substr(V(g)$name, 1, 1)

联合网络的邻接矩阵：

get.adjacency(g, sparse=FALSE, attr="Freq")

##    h1 h2 h3 h4 c2 c3 c4 c1 c5 p1 p2 p3
## h1  0  0  0  0  0  0  0  0  0  3  1  0
## h2  0  0  0  0  0  0  0  0  0  1  3  1
## h3  0  0  0  0  0  0  0  0  0  5  5  3
## h4  0  0  0  0  0  0  0  0  0  1  2  0
## c2  1  2  1  2  0  0  0  0  0  0  0  0
## c3  1  2  3  3  0  0  0  0  0  0  0  0
## c4  3  5  1  2  0  0  0  0  0  0  0  0
## c1  0  4  0  2  0  0  0  0  0  0  0  0
## c5  0  2  0  2  0  0  0  0  0  0  0  0
## p1  0  0  0  0  0  0  0  0  0  0  0  0
## p2  0  0  0  0  0  0  0  0  0  0  0  0
## p3  0  0  0  0  0  0  0  0  0  0  0  0

以图形方式

t <- match(V(g)$type, c("p", "h", "c") )
plot(g, vertex.color=t )

甚至

l <- layout_with_fr(g, miny=t, maxy=t  )
plot(g, vertex.color=t, layout=l, edge.width=E(g)$Freq)