Corrected bugs in J=4 scripts

main
Corentin Choisy 12 months ago
parent 04a51d794d
commit 59d753d8e3

@ -931,16 +931,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

@ -931,16 +931,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

@ -931,16 +931,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