在R中使用Statnet包导入加权邻接矩阵

Question

我试图使用Unix服务器上的二分和statnet套件(包括网络、sna和其他几个包)来分析R中2模式网络的加权1模式投影。使用二部和矩阵代数的组合，投影工作得很好，但是当我试图将值矩阵作为加权网络对象导入时，使用下面的代码，我似乎放松了原始矩阵中的值。

MNDocnet<-as.network(MNDocmatrix,matrix.type="adjacency",directed=FALSE, hyper=FALSE, loops=TRUE, multiple=FALSE, bipartite = FALSE, ignore.eval=FALSE, names.eval="patients")

谢谢你能提供的任何帮助。

Answer 1

如果没有数据结构，很难准确地知道，但在我看来，这种语法是正确的。下面是一个例子

制作样本输入数据

> adjmat<-matrix(c(0,1,2,3,0,4,5,6,0),ncol=3)
> adjmat
     [,1] [,2] [,3]
[1,]    0    3    5
[2,]    1    0    6
[3,]    2    4    0

将矩阵转换为网络对象

> test<-as.network(adjmat,matrix.type='adjacency',ignore.eval=FALSE,names.eval='sample')

打印名为“sample”的属性的边缘值

> test%e%'sample'
[1] 1 2 3 4 5 6

请注意，如果要将其转换回值矩阵，则需要给出提供值的属性的名称：

> as.matrix(test)
  1 2 3
1 0 1 1
2 1 0 1
3 1 1 0

与

> as.matrix(test,attrname='sample')
  1 2 3
1 0 3 5
2 1 0 6
3 2 4 0