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 3 40
MARK 0
SEED 1
ERASE C994-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 -0.226000 C998
JOIN C998 0.257000 C998
JOIN C998 0.146000 C998
JOIN C996 0.156000 C996
JOIN C996 0.839000 C996
NOTE put MVN variances for predictors in model
JOIN c994 0.120000 c994
JOIN c994 0.000000 c994
JOIN c994 0.000000 c994
JOIN c995 0.125000 c995
JOIN c995 0.045000 c995
JOIN c995 0.249000 c995
LINK 10 G40
LINK 2 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 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"
CALC b41 = 1000
LOOP b22 20 300 20
    OBEY simu2.txt
ENDL
MARK 1