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 3 40
MARK 0
SEED 1
ERASE C994-C998
NOTE setup the values of beta, sigma2u, sigma2e etc.
JOIN C998 -0.103000 C998
JOIN C998 0.170000 C998
JOIN C998 0.591000 C998
JOIN C996 0.642000 C996
NOTE put MVN variances for predictors in model
JOIN c994 0.240000 c994
JOIN c994 0.026000 c994
JOIN c994 1.000000 c994
LINK 10 G40
LINK 1 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 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"
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 b23 50 1500 50
    OBEY setup.txt
ENDL
MARK 1