NOTE MLwiN macro code generated by MLPowSim
NOTE b23 - number of units
ERASE c1011 c1012
GENErate 1 b23 c1
PUT b23 1 c4
PUT b23 1 c5
NAME c1 'l1id' c4 'cons' c5 'resp'
CLEAr
RESP c5
IDEN 1 c1
EXPL 1 c4
SETV 1 c4
ERROR 0
BATCH 1
ERASE G33
LOOP b40 1 b41
    PICK 1 c998 b51
    EDIT 1 c1098 b51
    PICK 1 c996 b51
    EDIT 1 c1096 b51
    PUT b23 1 c5
    SIMU c5
    START
    NOTE We use zeroes for residual starting values
    MCRS 1
    MCMC 0 1000 1 5.8 50 10 1 1 1 1 1 1
    MCMC 1 5001 1 c1090 c1091 c1003 c1004 1 1
    PUPN c1003 c1004
    CODE 2 1 5001 c1002
    SPLIT c1090 c1002 G31 G34
    SORT 1 G31[1] G31[1] G31[1] G31[1]
    PICK 4875 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
PAUSE 1