kathy_5900
2008-10-04, 09:24
[问题] 假定有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次.
请教高手,如何实现?
[分析] 以下是本人模拟1次的编程,但不知如何实现50次模拟,也想到做循环,但不知如何实现,请版主及其他高手们指点啊!! 急~~~~~~
>> n=2000;
Rho=[1 0.25 0 -0.25;...
0.25 1 0 0;...
0 0 1 0;...
-0.25 0 0 1];
Z=mvnrnd([0 0 0 0],Rho,n);
U=normcdf(Z,0,1);
n=1;
X=[binoinv(U(:,1),n,0.3) binoinv(U(:,2),n,0.3) binoinv(U(:,3),n,0.3)...
binoinv(U(:,4),n,0.3) ];
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次.
请教高手,如何实现?
[分析] 以下是本人模拟1次的编程,但不知如何实现50次模拟,也想到做循环,但不知如何实现,请版主及其他高手们指点啊!! 急~~~~~~
>> n=2000;
Rho=[1 0.25 0 -0.25;...
0.25 1 0 0;...
0 0 1 0;...
-0.25 0 0 1];
Z=mvnrnd([0 0 0 0],Rho,n);
U=normcdf(Z,0,1);
n=1;
X=[binoinv(U(:,1),n,0.3) binoinv(U(:,2),n,0.3) binoinv(U(:,3),n,0.3)...
binoinv(U(:,4),n,0.3) ];