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 4 40
MARK 0
SEED 1
ERASE C994-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 -0.167000 C998
JOIN C998 0.166000 C998
JOIN C998 0.165000 C998
JOIN C998 0.560000 C998
JOIN C996 0.081000 C996
JOIN C996 0.562000 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 c994 0.020000 c994
JOIN c994 0.000000 c994
JOIN c994 0.902000 c994
JOIN c995 0.125000 c995
JOIN c995 0.045000 c995
JOIN c995 0.249000 c995
JOIN c995 0.013000 c995
JOIN c995 -0.006000 c995
JOIN c995 0.116000 c995
LINK 10 G40
LINK 2 G21
LINK 4 G22
LINK 4 G23
LINK 4 G24
LINK 4 G25
LINK 4 G26
LINK 4 G27
LINK 4 G31
LINK 4 G32
LINK 4 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"
NAME G22[4] "zpow3"
NAME G23[4] "zlow3"
NAME G24[4] "zupp3"
NAME G25[4] "spow3"
NAME G26[4] "slow3"
NAME G27[4] "supp3"
CALC b41 = 1000
LOOP b22 10 150 10
    OBEY simu2.txt
ENDL
MARK 1