NOTE MLwiN macro code generated by MLPowSim
CODE 2 1 b41 c1002
SPLIT c1011 c1002 G31 G34
SPLIT c1012 c1002 G32 G34
NOTE calculate IGLS interval coverage
NED 0.975000 b42
JOIN "Samplesize" b23 "Samplesize"
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.405500)/b203
CALC b203 = b203+b42
CALC b206=(-0.405500)/b206
CALC b206 = b206+b42
CALC b207=(-0.405500)/b207
CALC b207 = b207+b42
NPRO b203 b204
JOIN "spow0" b204 "spow0"
NPRO b206 b204
JOIN "slow0" b204 "slow0"
NPRO b207 b204
JOIN "supp0" b204 "supp0"