周辺度数に応じて要素を振る

またしてもこちらから。
行列Mの縦横の和が等しい行列をランダムに生成して、モンテカルロ法でMに偏りがあるかないかを調べたい件。
r2dtableという関数を使ってみる。

set.seed(100)
r2dtable(1,rowSums(M),colSums(M))

要素は0か1にしたい。さてどうしたものやら。