とある身内のご要望にお応えして、stochastic block model (SBM)を生成するコードを載せます。 Poisson次数分布の、一番基本的なSBMを生成します。 Affinity matrix (p_rs)の要素を陽に指定して、任意のブロック構造を作れます。 cinとcoutは、rescaleされた…

速度を比べっこ。 ランダムグラフのedgelist(links)を作ってみて、「各頂点の隣接頂点の配列」を要素に持つ、入れ子配列(nbsimpleとnb)を２通りで作ってみる。 links = round(Int64,10000*rand(50000,2)+1) links = vcat(links,hcat(links[:,2],links[:,1]))…

Indicator行列 例えば、各要素は１〜３番目のどれかが１、他は０になっていて、 要素１：１番目が１、他は０ 要素２：３番目が１、他は０ 要素３：２番目が１、他は０ 要素４：１番目が１、他は０ というとき、 1 0 0 0 0 1 0 1 0 1 0 0 という行列。 作り方…

どこにでも書いてある、DFSで連結成分を出してくるという操作のサンプル。 特にJuliaだから特殊ということはないはず。 枝の集合linksが与えられたときに、連結成分の頂点CCnodesを返す。 Stackでやる方法 function dfs_stack(nb,root) visited = Int64[] st…