Good morning,

Could you, please, help me find the formula used by Stata to compute predicted probabilities after fitting a logistic regression model?
I tried to guess it from the logit formula but my results (p1 or p2) are different from those obtained by Stata (pr).

Thanks in advance!
--
version 17

Code:
* Example generated by -dataex-. For more info, type help dataex
clear
input byte age double(vitesserisk irrespect sensation danger) byte accidenté
21                  4  4.833333333333333  4.972222222222222  4.103174603174604 1
21 3.2000000000000006 3.1666666666666665  4.472222222222222 3.6111111111111107 0
18  4.333333333333333 2.6666666666666665  4.083333333333333  3.492063492063492 1
24  4.733333333333333 2.3333333333333335  4.111111111111111   2.76984126984127 0
22  5.066666666666666 3.6666666666666665  3.027777777777778 1.5317460317460316 1
23                4.2 2.6666666666666665  2.027777777777778  2.261904761904762 1
22                4.8  4.833333333333333 3.1111111111111107 2.2698412698412698 1
20  4.866666666666666                  3 2.5833333333333335 1.8333333333333333 1
20  4.866666666666666                4.5 2.4444444444444446 2.0555555555555554 1
24  4.733333333333333 2.3333333333333335 2.5833333333333335 2.4603174603174605 1
end
label values accidenté labels3
label def labels3 0 "Non", modify
label def labels3 1 "Oui", modify

Code:
logistic accidenté age vitesserisk irrespect sensation danger, coef
predict pr
predict lp, xb
gen logitp=1.131809 + (-0.278859*age) + (0.4179795*vitesserisk) + (-0.2565562*irrespect) + (-0.0087661*sensation) + (-0.4366206*danger)
list pr lp logitp age Éducation revenu emploi vitesserisk irrespect sensation danger accidenté in 1/5
gen p1=1/(1+exp(-logitp))
gen p2=100/(1+exp(-logitp))
list pr lp logitp p1 p2 age vitesserisk irrespect sensation danger accidenté in 1/5