计算R中随机网络的r平方值

Question

I created a set of random networks in R using the erdos-renyi algorithm. However, I cannot find a function to calculate the R-squared (Pearson correlation coefficient) value of the graph. Fro example, if I use the function erdos renyi game in the igraph package to make my network

erdos.renyi.game(281, 447, type = "gnm",m = 465, directed = FALSE)

How would I be able to calculate the R-squared value for this?

Answer 1

library(igraph)
library(intergraph)
library(sna)

# create the graphs
gs <- replicate(2,  
                sample_gnm(281, 447,m = 465, directed = FALSE),
                simplify = F) 
# sample_gnm is preferred to erdos.renyi.game now

# turn them into network objects
nets <- lapply(gs, asNetwork)

Now use gcor to get the correlation and qaptest to test the significance.

gcor(nets, g1 = 1, g2 =2)
[1] 0.003537025
qaptest(nets, gcor, g1 = 1, g2 = 2)

QAP Test Results

Estimated p-values:
    p(f(perm) >= f(d)): 0.17 
p(f(perm) <= f(d)): 0.834