NOTE MLwiN macro code generated by MLPowSim
NOTE This is outer code to be run directly in MLwiN
NOTE You will also need simu2.txt, simu3.txt, genresp.txt, setup.txt and analyse.txt
INIT 5 5000 5000 3 40
MARK 0
SEED 1
ERASE C994-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 0.500000 C998
JOIN C996 1.200000 C996
JOIN C700 1.200000 C700
JOIN C996 0.400000 C996
JOIN C750 0.400000 C750
JOIN C996 8.000000 C996
LINK 10 G40
LINK 3 G21
LINK 1 G22
LINK 1 G23
LINK 1 G24
LINK 1 G25
LINK 1 G26
LINK 1 G27
LINK 1 G28
LINK 1 G29
LINK 1 G30
LINK 1 G31
LINK 1 G32
LINK 1 G33
LINK 3 G34
LINK 0 G40
NAME G21[3] "N-XC Fact1"
NAME G21[1] "N-level 1"
NAME G21[2] "N-XC Fact2"
NAME G22[1] "zpow0"
NAME G23[1] "zlow0"
NAME G24[1] "zupp0"
NAME G25[1] "spow0"
NAME G26[1] "slow0"
NAME G27[1] "supp0"
NAME G28[1] "mpow0"
NAME G29[1] "mlow0"
NAME G30[1] "mupp0"
CALC b41 = 1000
LOOP b30 20 100 20
    OBEY simu2.txt
ENDL
MARK 1