NOTE MLwiN macro code generated by MLPowSim CODE 6 1 b41 c30 CODE 13 1 b41 c31 SPLIT c1011 c30 c51-c56 SPLIT c1012 c31 c101-c113 NOTE calculate IGLS interval coverage NED 0.975000 b42 JOIN c209 b21 c209 JOIN c210 b22 c210 CALC c101 = c101*(c101>0) - c101*(c101<0) CALC c101 = sqrt(c101) CALC c200= c51 - b42*c101 CALC c201= c200>0 AVER c201 b202 b203 b204 JOIN c211 b203 c211 CALC b204 = (b203)*(1-b203)/b41 CALC b205 = b203-b42*sqrt(b204) JOIN c251 b205 c251 CALC b205 = b203+b42*sqrt(b204) JOIN c271 b205 c271 CALC c103 = c103*(c103>0) - c103*(c103<0) CALC c103 = sqrt(c103) CALC c200= c52 - b42*c103 CALC c201= c200>0 AVER c201 b202 b203 b204 JOIN c212 b203 c212 CALC b204 = (b203)*(1-b203)/b41 CALC b205 = b203-b42*sqrt(b204) JOIN c252 b205 c252 CALC b205 = b203+b42*sqrt(b204) JOIN c272 b205 c272 NOTE calculate IGLS SE method AVER c101 b202 b203 b204 b205 CALC b206= b203+b42*b205 CALC b207= b203-b42*b205 CALC b203=0.100000/b203 CALC b203 = (-1)*b203+b42 CALC b206=0.100000/b206 CALC b206 = (-1)*b206+b42 CALC b207=0.100000/b207 CALC b207 = (-1)*b207+b42 NPRO b203 b204 JOIN c231 b204 c231 NPRO b206 b204 JOIN c291 b204 c291 NPRO b207 b204 JOIN c311 b204 c311 AVER c103 b202 b203 b204 b205 CALC b206= b203+b42*b205 CALC b207= b203-b42*b205 CALC b203=0.600000/b203 CALC b203 = (-1)*b203+b42 CALC b206=0.600000/b206 CALC b206 = (-1)*b206+b42 CALC b207=0.600000/b207 CALC b207 = (-1)*b207+b42 NPRO b203 b204 JOIN c232 b204 c232 NPRO b206 b204 JOIN c292 b204 c292 NPRO b207 b204 JOIN c312 b204 c312