我正在用igraph绘制图形,我希望边缘具有不同的颜色,具体取决于它们所代表的连接的强度。我可以设置颜色,但我不能将它们与连接强度的值联系起来。 我目前的代码如下:
library(igraph)
library(raster)
library(ggplot2)
library(statnet)
library(qgraph)
connectivityMatrix <- as.matrix(read.table(file=myFile,sep='')))
coordinates <- as.matrix(read.table(file=coordinatesFile))
connectivityMatrix<-connectivityMatrix[1:833,1:833]
CM<-connectivityMatrix[subsetX,subsetY]
COORD<-coordinates[subset,]
net <- as.network(CM, matrix.type = "adjacency", directed = TRUE)
minX<-min(coordinates[,1])
maxX<-max(coordinates[,1])
minY<-min(coordinates[,2])
maxY<-max(coordinates[,2])
p<-plot(net, coord=COORD,xlim=c(minX,maxX),ylim=c(minY,maxY),edge.col=c('red','yellow','cyan','blue'),object.scale=0.005, vertex.col='dimgrey',edge.lwd=1)
在上面的代码中,有一种方法可以将使用edge.col指定的颜色与它们在CM中表示的值范围相关联吗?这样,对应于连接矩阵中的值0-x1的边缘将被绘制为红色,x1-x2被绘制为“黄色”,......和x3-x4被绘制为蓝色。 x1,x2,x3是范围限制,x4是CM的最大值。 有没有人知道如何做到这一点?是否可以添加包含边缘颜色和它们代表的值范围的图例?
提前感谢你,
答案 0 :(得分:6)
您可以使用colorRamp作为缩放功能。例如,请参阅下面的代码。
library(igraph)
#Create a random weighted graph
g = erdos.renyi.game(10,0.5)
E(g)$weight = runif(ecount(g))
#Color scaling function
c_scale <- colorRamp(c('red','yellow','cyan','blue'))
#Applying the color scale to edge weights.
#rgb method is to convert colors to a character vector.
E(g)$color = apply(c_scale(E(g)$weight), 1, function(x) rgb(x[1]/255,x[2]/255,x[3]/255) )
#plot using igraph
plot.igraph(g)