st: Bootstrapped Std Errors for Concentration Index Decomposition

[Abekah Nkrumah <ankrumah@gmail.com>]

over their determinants. I have calculated the concentration indicies
together with their standard errors via the code below. My difficulty
however is how to calculated bootstrapped standard errors for the  for
the components of the decomposition (especially the contributions and
concentration index of the X's). Please can anyone help?

egen raw_rank=rank(wealth_index_dec), unique
sort raw_rank
quietly sum sample_weight
gen wi= sample_weight/r(sum)
gen cusum= sum(wi)
gen wj=cusum[_n-1]
replace wj=0 if wj==.
gen rank=wj+0.5*wi

gen weight12=sample_weight
replace weight12=weight12*10^6
replace weight12= int(sample_weight)

qui sum  neghaz [fw=weight12]
scalar mean= r(mean)
cor  neghaz rank [fw=weight12], c
sca c= (2/mean)*r(cov_12)
sca list c

gen y= neghaz

global X ............................................

qui regr y $X [pw=weight12],vce(cluster hh_number)
sum y [aw=weight12]
sca m_y=r(mean)

foreach x of global X {
     qui  {
            sca b_`x' = _b[`x']
            corr rank `x' [aw=weight12], c
            sca cov_`x' = r(cov_12)
            sum `x' [aw=weight12]
            sca elas_`x' = (b_`x'*r(mean))/m_y
            sca CI_`x' = 2*cov_`x'/r(mean)
            sca con_`x' = elas_`x'*CI_`x'
            sca prcnt_`x' = con_`x'/CI*100
      }
      di "`x' elasticity:", elas_`x'
      di "`x' concentration index:", CI_`x'
      di "`x' contribution:", con_`x'
      di "`x' percentage contribution:", prcnt_`x'
}


Please can you help in terms of how to calculate the bootstrapped
Standard errors for the contributions from the decomposition?

Thank you very much

Gordon
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/