孤立頂点が消滅する閾値に関する実験

ソースコード

library(igraph)

Warning: package 'igraph' was built under R version 4.3.3

n=100 #頂点数
N=400 #サンプル数
cs=seq(1,1.8,0.1) #1から1.8まで0.1刻み

ratio=c()
for (c in cs){
  count=0
  for (i in 1:N){
    p=c*log(n)/n
    g=sample_gnp(n,p)
    if(diameter(g)<=log(n)){
      count=count+1
    }
  }
  ratio=c(ratio,count/N)  
}

グラフ

barplot(ratio,names.arg = cs)