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 2 40
MARK 0
SEED 1
ERASE C994-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 -0.405500 C998
JOIN C996 0.250000 C996
JOIN C996 0.000000 C996
LINK 10 G40
LINK 2 G21
LINK 1 G22
LINK 1 G23
LINK 1 G24
LINK 1 G25
LINK 1 G26
LINK 1 G27
LINK 1 G31
LINK 1 G32
LINK 1 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"
CALC b41 = 1000
LOOP b22 10 50 5
    OBEY simu2.txt
ENDL
MARK 1