Corrected bug in J=4 scenarios

main
Corentin Choisy 12 months ago
parent 85a30b71ae
commit 04a51d794d

@ -344,16 +344,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -400,16 +400,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
@ -449,16 +449,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -512,16 +512,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
@ -561,16 +561,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -987,16 +987,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
@ -1036,16 +1036,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -1099,16 +1099,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
@ -1148,16 +1148,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear

@ -344,16 +344,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -400,16 +400,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
@ -449,16 +449,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -512,16 +512,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
@ -561,16 +561,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -987,16 +987,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
@ -1036,16 +1036,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -1099,16 +1099,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
@ -1148,16 +1148,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear

@ -344,16 +344,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -400,16 +400,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
@ -449,16 +449,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -512,16 +512,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
@ -561,16 +561,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.3,0.3,0.3 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.3,0.3,0.3 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.3,0.3,0.3 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -987,16 +987,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
@ -1036,16 +1036,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@ -1099,16 +1099,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
@ -1148,16 +1148,16 @@ forvalues replication = 1/1000 {
mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84) mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
if `difi'==1 { if `difi'==1 {
mat D= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0) mat E= (0.5,0.5,0.5 \ 0,0,0 \ 0,0,0 \ 0,0,0)
} }
if `difi'==2 { if `difi'==2 {
mat D= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0) mat E= (0,0,0 \ 0.5,0.5,0.5 \ 0,0,0 \ 0,0,0)
} }
if `difi'==3 { if `difi'==3 {
mat D= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0) mat E= (0,0,0 \ 0,0,0 \ 0.5,0.5,0.5 \ 0,0,0)
} }
if `difi'==4 { if `difi'==4 {
mat D= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5) mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
} }
mat D=B+E mat D=B+E
qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear

Loading…
Cancel
Save