2014-02-12 52 views
3

计算相配我有一个数据集如下:在IGRAPH

set.seed(123) 
A = data.frame(rnorm(10),rnorm(10),rnorm(10),rnorm(10)) 

,然后用于IGRAPH包进行网络出了以下内容:

inv<-cor(t(A)) 
inv[inv<0.5] <- 0 
inv[inv==1] <- 0 
g1 <- graph.adjacency(inv, mode = "undirected", diag=FALSE, weighted=TRUE) 

我们计算G1的相配coefficeint ,

assortativity (g1, types1, types2 = NULL, directed = TRUE) 

我现在的问题是,我应该如何设置“types”,它在文档中说它是ver tex值。究竟是什么意思?我想计算任何5个顶点网络的协调性。谁能告诉我这是怎么完成的?

+0

....基于什么的assortativity?这个问题的答案可能有助于你理解为什么有这些类型... – user1320502

+0

例如,基于一组选定的顶点(网络底层中的特定节点不使用它们的顶点标签)?你能告诉我如何做到这一点吗? – user2258452

回答

4

所以我想你想要的标称版本的assortativity。例如。

V(g1)$foo <- sample(1:3, replace=TRUE, vcount(g1)) 
assortativity.nominal(g1, types=V(g1)$foo) 
# [1] -0.2270916 

类型必须是从1开始的整数。请参阅文档中的详细信息。

+0

对于有向图的连续变量,其中一个会如何?我有一个关联度(范围从0到1)作为两个向量的入边和出边,但我似乎无法将它们添加到args类型1和类型2。 –

+0

您需要对其进行分离,我很害怕。对于连续数量,AFAIK assortativity没有定义。虽然我可以想象,推广它并不难。 –