r - 相関バイナリ変数を生成するための限界確率を計算する方法

Question

それぞれが相互に相関する長さ 100 のバイナリ変数を 50 個作成したいとします。

指定された rho を使用して相関行列を作成します。

cor.mat <- matrix(0.2,nrow=50, ncol=50)
diag(cor.mat) <- 1

次に使用しますrmvbin：

library(bindata)
rmvbin<-rmvbin(100, margprob=x, bincorr=cor.mat)

ただし、margprob 引数の計算方法がわかりません。誰か助けてくれませんか？

各行と列の確率の合計のベクトルである必要がありますか?

Answer 1

margprob残りの変数とは関係なく、任意の 1 つのバイナリ変数が 1 である確率の単純な繰り返しベクトルであるべきです。この値を呼び出しますp。変数が同一分布であると仮定すると (相関行列が当てはまると思われます)、 margprob=rep(p,50).

相関行列を使用して限界確率を決定することはできないため、各行と列の確率の合計のベクトルであってはなりません。確率変数の限界確率が何であるかを理解するのに問題がある場合は、問題についてより多くのコンテキストを提供する必要があり、 math.stackexchange.comにより適した質問になります。