数据模拟编程问题

shiyiming

假定有4个离散变量（x1~x4）均服从二项分布，这四个变量间相关系数矩阵如下:
        X1        X2        X3        x4
X1        1        0.25        0        -0.25
X2        0.25        1        0        0
X3        0        0        1        0
X4        -0.25        0        0        1
因变量Y(取值为0或1)与自变量(x1~x4)之间满足以下logistic模型:
logit(Pr(Y=1))= -0.6+0.25x1+025x2+0.5(x1x2)

想要产生n=100的样本,共模拟50次,生成50个SAS数据集.
请教高手,如何实现?