[求助]复杂数据模拟编程的问题

[kathy_5900]

假定有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(P|Y=1)=(-0.6+0.25x1+025x2+0.5(x1x2)),

想要产生样本量n=100的样本,共模拟50次. 请教高手,如何实现?

[kathy_5900]

模型有误已修改,请版主及各位高手予指点,急!!!!