計算R中隨機網絡的r平方值

Question

我使用erdos-renyi算法在R中創建了一組隨機網絡。 但是，我找不到函數來計算圖的R平方（皮爾森相關系數）值。 例如，如果我使用igraph軟件包中的功能erdos renyi game來建立我的網絡

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

我將如何為此計算R平方值？

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)

現在，使用gcor獲得相關性，並使用qaptest檢驗其重要性。

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