program lprob2

lincom `0' + _cons
local LN_O = ln(r(estimate))
local SE_LN_O = r(se)/r(estimate)
local CIL = exp(`LN_O' - 1.96*`SE_LN_0')
local CIU = exp(`LN_O' + 1.96*`SE_LN_0')
local PROB = r(estimate)/(1+r(estimate))
local CIL = `CIL'/(1 + `CIL')
local CIU = `CIU'/(1 + `CIU')
display " "
display " Probability: " `PROB' "     95% CI: " `CIL' " - " `CIU'

end