Рассчитать значение r-квадрата для случайной сети в R

Question

Я создал множество случайных сетей в R, используя алгоритм erdos-renyi. Однако я не могу найти функцию для вычисления значения R-квадрата (коэффициента корреляции Пирсона) графика. Например, если я использую игру erdos renyi в пакете igraph, чтобы создать свою сеть

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