2016-05-15 179 views
3

简而言之,我想根据相关强度绘制边缘,删除不重要的值后。我能做到这一点的正相关性对与edge.betweeness,但不幸的是不是否定:igraph不适用edge.width负相关系数

data <- matrix(rnorm(100),10,10)  
colnames(data) <- LETTERS[1:10] 

library(Hmisc) 
cor1 <- rcorr(data) 
diag(cor1$r) <- 0 

library(igraph) 

#####Example 1: 
graph <- graph.adjacency(cor1$r, weighted=TRUE, mode="lower") 
graph <- delete.edges(graph, E(graph)[ weight < 0.3 & weight > -0.3]) 
#####trying to pass edge weights to edge.width 
plot.igraph(graph, vertex.size=20, edge.width=E(graph)$weight) 
###edge.width=E(graph)$weight is ignored 

#####Example 2: 
graph <- graph.adjacency(cor1$r, weighted=TRUE, mode="lower") 
graph <- delete.edges(graph, E(graph)[ weight < 0.3]) #omitting the 2nd condition 
E(graph)$weight <- edge.betweenness(graph) #apparently required 
plot.igraph(graph, vertex.size=20, edge.width=E(graph)$weight) 
####this does work, but only for positive correlation coefficients 

#####Example 3: 
graph <- graph.adjacency(cor1$r, weighted=TRUE, mode="lower") 
graph <- delete.edges(graph, E(graph)[ weight < 0.3 & weight > -0.3]) 
E(graph)$weight <- edge.betweenness(graph) 
#####gives error: Error in .Call("R_igraph_edge_betweenness", graph, directed, weights, : 
#################At centrality.c:2046 : Weight vector must be non-negative, Invalid value 

那么,如何传递的负相关性值edge.width

回答

4

您可以使用有色边缘来表示负相关和正相关,并使用edge.width来表示相关的大小。在下面的例子中,我进行了如下修改:

  1. 改变图形对象的名称g1,因为是 在的igraph包的功能。
  2. 为简洁起见,请使用边缘删除条件的绝对值。
  3. edge.width参数更改为abs(E(g1)$weight)*8。绝对值确保权重始终为正值。 乘以8只会使边缘宽度变大。
  4. 添加edge.color参数到颜色顶点蓝色为 正相关,红色为负相关。

#####Example 3: 
g1 <- graph.adjacency(cor1$r, weighted=TRUE, mode="lower") 
g1 <- delete.edges(g1, E(g1)[ abs(weight) < 0.3 ]) 
plot.igraph(g1, vertex.size=20, edge.width=abs(E(g1)$weight)*8, 
     edge.color=ifelse(cor1$r > 0, "blue","red")) 

enter image description here

+0

非常感谢你。卓见。 – nouse