NOTE MLwiN macro code generated by MLPowSim
CODE 6 1 b41 c1002
SPLIT c1011 c1002 G31 G34
SPLIT c1012 c1002 G32 G34
NOTE calculate IGLS interval coverage
NED 0.975000 b42
JOIN "N-level 3" b30 "N-level 3"
JOIN "N-level 1" b21 "N-level 1"
JOIN "N-level 2" b22 "N-level 2"
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"
CALC c1005= G31[2] - b42*G32[2]
CALC c1006= G31[2] + b42*G32[2]
CALC c1007 = (c1005 > 0 & c1006 > 0) | (c1005 < 0 & c1006 < 0)
AVER c1007 b202 b203 b204
JOIN "zpow1" b203 "zpow1"
CALC b204 = (b203)*(1-b203)/b41
CALC b205 = b203-b42*sqrt(b204)
JOIN "zlow1" b205 "zlow1"
CALC b205 = b203+b42*sqrt(b204)
JOIN "zupp1" b205 "zupp1"
CALC c1005= G31[3] - b42*G32[3]
CALC c1006= G31[3] + b42*G32[3]
CALC c1007 = (c1005 > 0 & c1006 > 0) | (c1005 < 0 & c1006 < 0)
AVER c1007 b202 b203 b204
JOIN "zpow2" b203 "zpow2"
CALC b204 = (b203)*(1-b203)/b41
CALC b205 = b203-b42*sqrt(b204)
JOIN "zlow2" b205 "zlow2"
CALC b205 = b203+b42*sqrt(b204)
JOIN "zupp2" b205 "zupp2"
NOTE calculate IGLS SE method
AVER G32[1] b202 b203 b204 b205
CALC b206= b203+b42*b205
CALC b207= b203-b42*b205
CALC b203=21.500000/b203
CALC b203 = (-1)*b203+b42
CALC b206=21.500000/b206
CALC b206 = (-1)*b206+b42
CALC b207=21.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"
AVER G32[2] b202 b203 b204 b205
CALC b206= b203+b42*b205
CALC b207= b203-b42*b205
CALC b203=3.000000/b203
CALC b203 = (-1)*b203+b42
CALC b206=3.000000/b206
CALC b206 = (-1)*b206+b42
CALC b207=3.000000/b207
CALC b207 = (-1)*b207+b42
NPRO b203 b204
JOIN "spow1" b204 "spow1"
NPRO b206 b204
JOIN "slow1" b204 "slow1"
NPRO b207 b204
JOIN "supp1" b204 "supp1"
AVER G32[3] b202 b203 b204 b205
CALC b206= b203+b42*b205
CALC b207= b203-b42*b205
CALC b203=(-6.000000)/b203
CALC b203 = b203+b42
CALC b206=(-6.000000)/b206
CALC b206 = b206+b42
CALC b207=(-6.000000)/b207
CALC b207 = b207+b42
NPRO b203 b204
JOIN "spow2" b204 "spow2"
NPRO b206 b204
JOIN "slow2" b204 "slow2"
NPRO b207 b204
JOIN "supp2" b204 "supp2"