if (FALSE) {
set.seed(42)
## test under H0
a <- b <- G <- list()
for(i in 1:10){
a[[i]] <- igraph::sample_gnp(50,0.5)
b[[i]] <- igraph::sample_gnp(50,0.5)
}
G <- list(a,b)
k1 <- fraiman.test(G)
k1
## test under H1
a <- b <- G <- list()
for(i in 1:10){
a[[i]] <- igraph::sample_gnp(50,0.5)
b[[i]] <- igraph::sample_gnp(50,0.6)
}
G <- list(a,b)
k2 <- fraiman.test(G)
k2
}
Run the code above in your browser using DataLab