Logit变换之后如何求p置信区间

如题，我在计算一个概率时，对其进行了对数变换，
logit(p/1-p)=α+βx＋u
现在我已经通过线性回归计算出了上式的各个系数以及logit(p/1-p)的预测值，想问一下各位，如何求解p的预测值和置信区间呢？？？如果有Stata程序更好，多谢~

看这里!

https://www.stata.com/support/faqs/statistics/prediction-confidence-intervals/

或遵循此代码:

predict error, stdp

generate lb = xb - invnormal(0.975)*error

generate ub = xb + invnormal(0.975)*error

generate plb = invlogit(lb)

generate pub = invlogit(ub)