孤立頂点が消滅するしきい値に関する実験

ソースコード

\[ p = \frac{\log n}{n} \]

library(igraph)

n=100 #頂点数
N=400 #サンプル数
cs=seq(1,1.8,0.1) 

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

グラフ