NOTE MLwiN macro code generated by MLPowSim
NOTE This is outer code to be run directly in MLwiN
NOTE You will also need 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.001000 C998
JOIN C998 0.595000 C998
JOIN C996 0.648000 C996
LINK 10 G40
LINK 1 G21
LINK 2 G22
LINK 2 G23
LINK 2 G24
LINK 2 G25
LINK 2 G26
LINK 2 G27
LINK 2 G31
LINK 2 G32
LINK 2 G33
LINK 1 G34
LINK 0 G40
NAME G21[1] "Samplesize"
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"
CALC b41 = 1000
LOOP b23 5 50 5
    OBEY setup.txt
ENDL
MARK 1