r - 異質周辺確率を持つバイナリ変数の相関

Question

異なる変数が異なる限界確率を持つ 50 の相関バイナリ変数を作成したいと思います。

まず、相関行列を作成します。

cor.mat=matrix(.9,nrow=50,ncol=50)
cor.mat[,9:11]=.1 
cor.mat[9:11,]=.1  
diag(cor.mat)=1 

次に、「rmvbin」を使用してデータを生成します。

library(bindata)
marg=rep(c(0.4,0.6),c(25,25))
a<-rmvbin(100, margprob=marg, bincorr=cor.mat)

ただし、次のエラーが表示されます。

Error in commonprob2sigma(commonprob, simulvals) : 
  Matrix commonprob not admissible.

また、margprob の代わりに commonprob を試しました。

common=rep(c(0.4,0.6),c(25,25))
a<-rmvbin(100, commonprob=common, bincorr=cor.mat)

しかし、私は得る:

Error in if (n != dim(commonprob)[2]) { : argument is of length zero
> 

私は何を間違っていますか？