/R/R: network randomization test igprah stakoverflow

R: network randomization test igprah stakoverflow

# NETWORK RANDOMIZATION TEST
# TEST IF FEATURES ARE UNUSUAL

# Create Network
g <- graph_from_data_frame(edges, directed=FALSE)
gd = edge_density(g)

# Network Average path length
g.apl <-mean_distance(g, directed=FALSE)

# Create one random graph with the same number of nodes and edges as g
g.random <- erdos.renyi.game(n = gorder(g), p.or.m = gd, type = "gnp")

# Generate 1000 random graphs
gl <- vector('list',1000)

for(i in 1:1000){
  gl[[i]] <- erdos.renyi.game(n = gorder(g), p.or.m = gd, type = "gnp")
}

# Calculate average path length of 1000 random graphs
gl.apls <- unlist(lapply(gl, mean_distance, directed = FALSE))

# Plot the distribution of average path lengths
hist(gl.apls, xlim = range(c(1.5, 10)))
abline(v = g.apl, col = "red", lty = 3, lwd = 2)

# Calculate the proportion of graphs with an average path length lower than our observed
mean(gl.apls < g.apl)

Source: stackoverflow.com

Similar codes

0 vote125 views1 answer
Pdlyr mutate string extract
R
0 vote126 views1 answer
R squared regression in r with ggplot
R
0 vote136 views1 answer
How to enter character string with quotes in r
R
0 vote118 views1 answer
Correlation matrix using factors r
R
0 vote127 views1 answer
Ways to examine model in r
R
0 vote137 views1 answer
R change column value conditionally
R
0 vote122 views2 answers
How to round all numeric column types in r
R
0 vote119 views1 answer
Replace na in a column with values from another df
R
1 vote124 views1 answer
Empty plot in r
R
0 vote123 views1 answer
Calculating RMSE, MAE, MSE, Rsquared manually in R
R