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 1 40
MARK 0
SEED 1
ERASE C994-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 -0.140000 C998
JOIN C996 1.051000 C996
LINK 10 G40
LINK 1 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 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 G28[1] "mpow0"
NAME G29[1] "mlow0"
NAME G30[1] "mupp0"
CALC b41 = 1000
LOOP b23 20 500 20
    OBEY setup.txt
ENDL
MARK 1