test_measures.R

In MattiaGirardi1997/Understanding-Complex-Networks:

library(data.table)
library(igraph)
library(DescTools)
library(ggplot2)
data <- fread("removed_loops/output/master_measures_removed_loops.csv")
i <- 1
net <- fread(sprintf("removed_loops/data/%s.csv", data[i, Name]))
g <- graph_from_data_frame(net, directed = F)
Gini(degree(g))
Gini(degree.distribution(g))
plot(master_data[NetworkDomain == "Social"]$GiniDegreeDistribution)
plot(master_data[NetworkDomain == "Non-Social"]$GiniDegreeDistribution)

plot(degree.distribution(g))


eb <- edge.betweenness.community(g)

plot(eb, g)


data <- fread("removed_loops/output/master_measures_removed_loops.csv")
for(i in 1:nrow(data)){
  name <- data[i, Name]
  net <- graph_from_data_frame(fread(sprintf("removed_loops/data/%s.csv", data[i, Name])), directed = F)
  GiniDegreeDistribution <- Gini(degree_distribution(net))
  if(i == 1){
    res <- data.table(data[i, 1:4], GiniDegreeDistribution)
  } else {
    res <- rbind(res, data.table(data[i, 1:4], GiniDegreeDistribution))
  }
}

data$GiniDegreeDistribution <- res$GiniDegreeDistribution
data <- master_data[-(103:111)]
write.table(data, file = "removed_loops/output/master_measures_removed_loops.csv",
            sep = ",", row.names = F)

### Constraint
data <- fread("removed_loops/output/removed_loops_table.csv")
for(i in 1:nrow(data)){
  name <- data[i, Name]
  net <- fread(sprintf("removed_loops/data/%s.csv", data[i, Name]))
  net <- graph_from_data_frame(net, directed = F)
  MeanConstraint <- mean(constraint(net))
  GiniConstraint <- Gini(constraint(net))
  if(i == 1){
    write.table(data.table(data[i, 1:4], MeanConstraint, GiniConstraint),
                file = "removed_loops/output/constraint.csv",
                row.names = F, sep = ",")
  } else {
    write.table(data.table(data[i, 1:4], MeanConstraint, GiniConstraint),
                file = "removed_loops/output/constraint.csv",
                row.names = F, sep = ",", col.names = F, append = T)
  }
  rm(net, MeanConstraint, GiniConstraint)
}

### Modularity
data <- fread("removed_loops/output/removed_loops_table.csv")
for(i in 1:nrow(data)){
  name <- data[i, Name]
  net <- fread(sprintf("removed_loops/data/%s.csv", data[i, Name]))
  net <- graph_from_data_frame(net, directed = F)
  Modularity <- modularity(edge.betweenness.community(net))
  if(i == 1){
    write.table(data.table(data[i, 1:4], Modularity),
                file = "removed_loops/output/modularity.csv",
                row.names = F, sep = ",")
  } else {
    write.table(data.table(data[i, 1:4], Modularity),
                file = "removed_loops/output/modularity.csv",
                row.names = F, sep = ",", col.names = F, append = T)
  }
  rm(net, MeanConstraint, GiniConstraint)
}





tab <- data[which(!data$Name %in% master_data$Name), c("Name", "NetworkDomain")]
res[res$Name %in% tab$Name]

data$Name %in% master_data$Name