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, setup.txt and analyse.txt
INIT 5 5000 5000 3 40
MARK 0
SEED 1
ERASE C990-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 21.500000 C998
JOIN C998 3.000000 C998
JOIN C998 -6.000000 C998
JOIN C996 12.000000 C996
JOIN C996 2.500000 C996
JOIN C996 140.000000 C996
LINK 10 G40
LINK 3 G21
LINK 3 G22
LINK 3 G23
LINK 3 G24
LINK 3 G25
LINK 3 G26
LINK 3 G27
LINK 3 G31
LINK 3 G32
LINK 3 G33
LINK 3 G34
LINK 0 G40
NAME G21[3] "N-level 3"
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"
CALC b41 = 1000
LOOP b30 10 40 10
    OBEY simu2.txt
ENDL
MARK 1