1

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

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

有什么建议么?

感谢您的关注!

4

1 回答 1

4

这个怎么样:

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

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)

在此处输入图像描述

于 2015-10-22T12:49:56.787 回答