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, setup.txt and analyse.txt
INIT 5 5000 5000 4 40
MARK 0
SEED 1
ERASE C994-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 -0.228000 C998
JOIN C998 0.262000 C998
JOIN C998 0.191000 C998
JOIN C998 0.123000 C998
JOIN C996 0.155000 C996
JOIN C996 0.839000 C996
LINK 10 G40
LINK 2 G21
LINK 4 G22
LINK 4 G23
LINK 4 G24
LINK 4 G25
LINK 4 G26
LINK 4 G27
LINK 4 G31
LINK 4 G32
LINK 4 G33
LINK 2 G34
LINK 0 G40
NAME G21[1] "N-level 1"
NAME G21[2] "N-level 2"
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 G22[2] "zpow1"
NAME G23[2] "zlow1"
NAME G24[2] "zupp1"
NAME G25[2] "spow1"
NAME G26[2] "slow1"
NAME G27[2] "supp1"
NAME G22[3] "zpow2"
NAME G23[3] "zlow2"
NAME G24[3] "zupp2"
NAME G25[3] "spow2"
NAME G26[3] "slow2"
NAME G27[3] "supp2"
NAME G22[4] "zpow3"
NAME G23[4] "zlow3"
NAME G24[4] "zupp3"
NAME G25[4] "spow3"
NAME G26[4] "slow3"
NAME G27[4] "supp3"
CALC b41 = 1000
LOOP b22 20 300 20
    OBEY simu2.txt
ENDL
MARK 1