如何使用 R 计算二次分配过程？

Question

我有 3 个网络，我想用二次分配程序来衡量结构等价性。我已经执行了这段代码但是得到了这个错误。任何人都可以帮助并建议我走对路吗？

 rm(list=ls())
    getwd()
    library(sna)
    dat=read.csv(file.choose(),header=TRUE)
    dat1=read.csv(file.choose(),header=TRUE) 
    dat2=read.csv(file.choose(),header=TRUE) 

    el=as.matrix(dat)
    g=graph.edgelist(el,directed=FALSE)

    el=as.matrix(dat1)
    g1=graph.edgelist(el,directed=FALSE)

    el=as.matrix(dat2)
    g2=graph.edgelist(el,directed=FALSE)
    #Perform qap tests of graph correlation
     q.12<-qaptest(g,gcor,g1=1,g2=2)
Error in FUN(X[[i]], ...) : 
  as.sociomatrix.sna input must be an adjacency matrix/array, network, or list.
> q.13<-qaptest(g,gcor,g1=1,g2=3)
Error in FUN(X[[i]], ...) : 
  as.sociomatrix.sna input must be an adjacency matrix/array, network, or list.
> 
> #Examine the results
> summary(q.12)
Error in summary(q.12) : object 'q.12' not found
> plot(q.12)
Error in plot(q.12) : object 'q.12' not found
> summary(q.13)
Error in summary(q.13) : object 'q.13' not found
> plot(q.13)
Error in plot(q.13) : object 'q.13' not found

Answer 1

两个网络的节点大小应该相同。意味着如果一个网络有 50 个节点，另一个网络应该有相同的节点集。但是，没有边可以不同

如何使用 R 计算二次分配过程？

How to compute Quadratic Assignment Procedure using R?

r

sna