モンテカルロパーミュテーション

比較する

Rには固定周辺度数を引数に分割表を乱数発生する関数r2dtable()がある(こちら) Patefield's algorithmを用いる 地道にやるのとPatefield's でやるのとの違いはどれくらい? ほとんど違わない Nr<-2 Nc<-2 M<-matrix(sample(1:100,Nr*Nc,replace=TRUE),Nr,Nc)…