Skip to content

islands

Daria-Maltseva edited this page Jul 14, 2020 · 2 revisions

Islands

AA 4th collaboration network (Ct')

In Pajek, I generated islands for Ct' network [5,100] and got 2062 for SNS and 1415 for SNA (our) datasets. I extracted them (Network + Partition - Extract Subnetwork [1-*] and saved.
Visualization:

setwd("C:/Mail.Ru Cloud/ANR HSE/ANR Projects/SNA Vlado Batagelj/Comparison/Analyses/WAl")
p <- read.csv("./AAlislands.clu",header=FALSE,skip=2)$V1
t <- table(p)
head(t)
tail(t)
t <- as.vector(table(p))
t1 <- table(t)
tl <- t1[2:53] # there was a problem that the value "4" appeared (with the value 1 - as there would be 1 island of 4 nodes); I did not take it. 
tl # checked 

setwd("C:/Mail.Ru Cloud/ANR HSE/ANR Projects/SNA Vlado Batagelj/Comparison/Analyses/WAbm")
p2 <- read.csv("./AAbmIslands.clu",header=FALSE,skip=2)$V1
t2 <- table(p2)
head(t2)
tail(t2)
t2 <- as.vector(table(p2))
tbm <- table(t2)
tbm # checked 

length(tl)
length(tbm)

x <- as.integer(names(tl))
length(x)
x1 <- as.integer(names(tbm))
length(x1)

yt <- c(1,5, 10, 30, 50,100,300, 500); yl <- c("1","5", "10", "30", "50","100","300", "500")
plot(x,tl,log="xy", cex=1,pch=16, col='blue', main="N of islands and nodes", yaxt="n", xlab = "N of nodes", ylab = "N of islands")
axis(side=2,at=yt, labels=yl)
points(x1,tbm,cex=1,pch=16, col='red')
legend("topright", 
       legend = c("SNS", "SNA"), 
       col = c('blue', 'red'), 
       pch = c(16,16), 
       bty = "n", 
       pt.cex = 2, 
       cex = 0.9, 
       text.col = "black", 
       horiz = F , 
       #inset = c(0.1, 0.1)
)

KK' network

In Pajek, I generated islands for KK' network [2,100] and got 426 islands for SNS(1), 292 for SNS(2), and 270 for SNA (our) datasets. I extracted them (Network + Partition - Extract Subnetwork [1-*] and saved.
Visualization:

setwd("C:/Mail.Ru Cloud/ANR HSE/ANR Projects/SNA Vlado Batagelj/Comparison/Analyses/WKl")
p <- read.csv("./KKl_norm_islands.clu",header=FALSE,skip=2)$V1
t <- table(p)
head(t)
tail(t)
t <- as.vector(table(p))
tl <- table(t)
tl #checked, there was the same case as above -- 1 island of 1 node; did not take it. 
tl <- tl[2:7]
tl 

setwd("C:/Mail.Ru Cloud/ANR HSE/ANR Projects/SNA Vlado Batagelj/Comparison/Analyses/WKls")
pls <- read.csv("./KKls_norm_islands.clu",header=FALSE,skip=2)$V1
tls <- table(pls)
tls
tls <- as.vector(table(pls))
tls <- table(tls)
tls 

setwd("C:/Mail.Ru Cloud/ANR HSE/ANR Projects/SNA Vlado Batagelj/Comparison/Analyses/WKbm")
pbm <- read.csv("./KKbm_norm_islands.clu",header=FALSE,skip=2)$V1
tbm <- table(pbm)
head(tbm)
tail(tbm)
tbm <- as.vector(table(pbm))
tbm <- table(tbm)
tbm

x1 <- as.integer(names(tl))
length(x1)
x2 <- as.integer(names(tls))
length(x2)
x3 <- as.integer(names(tbm))
length(x3)

yt <- c(1,10, 100, 300); yl <- c("1","10", "100", "300")
plot(x1,tl,log="xy", cex=1.2,pch=16, col='blue', main="N of islands and nodes", yaxt="n", xlab = "N of nodes", ylab = "N of islands")
axis(side=2,at=yt, labels=yl)
points(x2,tls,cex=1.2,pch=16, col='black')
points(x3,tbm,cex=1.2,pch=16, col='red')
legend("topright", 
       legend = c("SNS(1)", "SNS(2)", "SNA"), 
       col = c('blue', 'black', 'red'), 
       pch = c(16,16), 
       bty = "n", 
       pt.cex = 2, 
       cex = 0.9, 
       text.col = "black", 
       horiz = F , 
       #inset = c(0.1, 0.1)
)