library(bindata)
rmvbin(10, bincorr=cbind(c(1,0.5),c(0.5,1)), margprob=c(0.5,0.5))
[,1] [,2]
[1,] 1 1
[2,] 0 0
[3,] 0 1
[4,] 1 0
[5,] 0 0
[6,] 1 1
[7,] 0 0
[8,] 0 1
[9,] 1 1
[10,] 0 1
x12<-rmvbin(10, bincorr=cbind(c(1,0.5),c(0.5,1)), margprob=c(0.5,0.5))
> x1<-x12[,1]
> x2<-x12[,2]library(bindata)
rmvbin(10, bincorr=cbind(c(1,0