NOTE MLwiN macro code generated by MLPowSim
CODE 4 1 b41 c1002
SPLIT c1011 c1002 G31 G34
SPLIT c1012 c1002 G32 G34
NOTE calculate IGLS interval coverage
NED 0.975000 b42
JOIN "N-XC Fact1" b30 "N-XC Fact1"
JOIN "N-level 1" b21 "N-level 1"
JOIN "N-XC Fact2" b22 "N-XC Fact2"
CALC c1005= G31[1] - b42*G32[1]
CALC c1006= G31[1] + b42*G32[1]
CALC c1007 = (c1005 > 0 & c1006 > 0) | (c1005 < 0 & c1006 < 0)
AVER c1007 b202 b203 b204
JOIN "zpow0" b203 "zpow0"
CALC b204 = (b203)*(1-b203)/b41
CALC b205 = b203-b42*sqrt(b204)
JOIN "zlow0" b205 "zlow0"
CALC b205 = b203+b42*sqrt(b204)
JOIN "zupp0" b205 "zupp0"
NOTE calculate IGLS SE method
AVER G32[1] b202 b203 b204 b205
CALC b206= b203+b42*b205
CALC b207= b203-b42*b205
CALC b203=0.500000/b203
CALC b203 = (-1)*b203+b42
CALC b206=0.500000/b206
CALC b206 = (-1)*b206+b42
CALC b207=0.500000/b207
CALC b207 = (-1)*b207+b42
NPRO b203 b204
JOIN "spow0" b204 "spow0"
NPRO b206 b204
JOIN "slow0" b204 "slow0"
NPRO b207 b204
JOIN "supp0" b204 "supp0"
NOTE MCMC method
AVER G33[1] b202 b203 b204
JOIN "mpow0" b203 "mpow0"
CALC b204 = (b203)*(1-b203)/b41
CALC b205 = b203-b42*sqrt(b204)
JOIN "mlow0" b205 "mlow0"
CALC b205 = b203+b42*sqrt(b204)
JOIN "mupp0" b205 "mupp0"