2011-06-30から1日間の記事一覧

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

またしてもこちらから。 行列Mの縦横の和が等しい行列をランダムに生成して、モンテカルロ法でMに偏りがあるかないかを調べたい件。 r2dtableという関数を使ってみる。 set.seed(100) r2dtable(1,rowSums(M),colSums(M)) 要素は0か1にしたい。さてどうしたも…