version 7.0
use implementation, clear
*step2
set type double
list
*step3
blogit aff tot g e gei, or
lrtest, saving(0)
blogit aff tot g e, or
lrtest
return list
*step4
display (1-nchi2(r(df), r(chi2)*500/10^8, invchi2(r(df), (1-0.05))))*100
*step5
display round(10^8*(npnchi2(r(df), invchi2(r(df), (1-0.05)), (1-(80 /100)))/ r(chi2)), 1)