看这里!
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)