NOTE MLwiN macro code generated by MLPowSim
NOTE b30 - number of first factor units, b22 - number of second factor units
NOTE b21 - number of replications
CALC b23 = b30*b21*b22
CALC b24 = b21*b30
ERASE c1011 c1012
GENErate 1 b23 c1
CODE b22 b24 1 c2
CODE b30 b21 b22 c3
PUT b23 1 c4
PUT b23 1 c5
NAME c1 'l1id' c2 'xcid2' c3 'xcid1' c4 'cons' c5 'resp'
CLEAr
RESP c5
IDEN 3 c3
IDEN 2 c2
IDEN 1 c1
EXPL 1 c4
SETV 1 c4
SETV 2 c4
SETV 3 c4
ERROR 0
BATCH 1
XCLA 1
ERASE G33
LOOP b40 1 b41
    PICK 1 c998 b51
    EDIT 1 c1098 b51
    PICK 1 c996 b51
    EDIT 1 c1096 b51
    PICK 2 c996 b51
    EDIT 2 c1096 b51
    PICK 3 c996 b51
    EDIT 3 c1096 b51
    OBEY genresp.txt
    START
    NOTE Put in correct starting values for variances
    PICK 1 c996 b51
    EDIT 1 c1096 b51
    PICK 2 c996 b51
    EDIT 2 c1096 b51
    PICK 3 c996 b51
    EDIT 3 c1096 b51
    NOTE We use zeroes for residual starting values
    MCRS 1
    MCMC 0 5000 1 5.8 50 10 1 1 1 1 1 1
    MCMC 1 10000 1 c1090 c1091 c1003 c1004 1 1
    PUPN c1003 c1004
    CODE 4 1 10000 c1002
    SPLIT c1090 c1002 G31 G34
    SORT 1 G31[1] G31[1] G31[1] G31[1]
    PICK 250 G31[1] b77
    CALC b78= b77>0
    JOIN G33[1] b78 G33[1]
    JOIN c1098 c1096 c1011 c1011
    SEPICK c1001
    JOIN c1001 c1012 c1012
ENDL
OBEY analyse.txt
XCLA 0
PAUSE 1