상관 관계가있는 이진 데이터를 R?

Question

지정된 파이 계수를 가진 2 진수 데이터의 2 벡터를 원한다면 어떻게 R을 사용하여 시뮬레이션 할 수 있습니까?

예를 들어, 내가 어떻게 두 x 같은 벡터와 0.79

> x = c(1, 1, 0, 0, 1, 0, 1, 1, 1) 
> y = c(1, 1, 0, 0, 0, 0, 1, 1, 1) 
> cor(x,y) 
[1] 0.7905694 

Answer 1

패키지이 더 복잡한 이진 데이터를 생성하는 좋은이다 bindata의 효율적인 코르와 지정된 벡터 길이의 y를 만들 수 있습니다 상관 구조.

library(bindata) 

## Construct a binary correlation matrix 
rho <- 0.7905694 
m <- matrix(c(1,rho,rho,1), ncol=2) 

## Simulate 10000 x-y pairs, and check that they have the specified 
## correlation structure 
x <- rmvbin(1e5, margprob = c(0.5, 0.5), bincorr = m) 
cor(x) 
#   [,1]  [,2] 
# [1,] 1.0000000 0.7889613 
# [2,] 0.7889613 1.0000000