Dear Stata experts.

using
Code:
webuse auto
// Binary outcome dominance analysis with factor varaible (e.g., Azen & Traxel, 2009)
Code:
domin foreign trunk weight, reg(logit) fitstat(e(r2_p)) sets((i.rep78))
// That order gives the ranking of variables based on their importance, from there, calculating the bootstrapping as in the help file is by mean of

Code:
bootstrap, reps(500): domin foreign trunk weight, reg(logit) fitstat(e(r2_p)) sets((i.rep78))
// following each method i end up having same results which is not what i am looking for results

General dominance statistics: Logistic regression
Number of obs = 59
Overall Fit Statistic = 0.5229

| Dominance Standardized Ranking
foreign | Stat. Domin. Stat.
------------+------------------------------------------------------------------------
trunk | 0.0834 0.1595 3
weight | 0.2611 0.4994 1
set1 | 0.1783 0.3411 2


My question is:
- Using bootstrapping, how would it be possible to calculate the mean for dominance value for bootstrapped samples for the 500 samples.
- How to calculate the a “reproducibility” value,which represents the proportion of bootstrapped samples that the given dominance pattern is produced. The closer the reproducibility value is to 1.00, the greater the robustness of the dominance results.

Much appreciate your help.