**** TEM QUE ARRUMAR O ENDEREÇO PARA SALVAR O ARQUIVO DE RESULTADOS
log using ".....\indicadores_desigualdade.log", replace


*** ABRIR O BANCO DE DADOS... É O ARQUIVO DTA!!!
use ".....\xxxxxxxx.dta", clear


replace V4614 = . if V4614> 400000
replace V4621 = . if V4621> 400000



*** cálculo do "gini"

sort V4621
gen dom_i = _n if V4621~=.
egen tot_domic = max(dom_i)
gen pi = dom_i/tot_domic

egen tot_renda = sum(V4621)
gen double parc_renda = V4621/tot_renda

gen double fi = parc_renda if _n==1
replace fi = fi[_n-1]+parc_renda if _n>1 & V4621~=.

gen termo_gini = fi if _n==1
replace termo_gini = fi + fi[_n-1] if _n > 1 & V4621~=.

egen soma_termo_gini = sum(termo_gini)

gen gini = 1 - (soma_termo_gini/tot_domic)



*** cálculo dos percentiles da renda
sort V4621
pctile pct = V4621, nq(10) genp(decil)
list decil pct in 1/10

*** calculo da renda dos 50% mais pobres e 10% mais ricos
/*

     +--------------+
     | decil    pct |
     |--------------|
  1. |    10    111 |
  2. |    20    183 |
  3. |    30    250 |
  4. |    40    323 |
  5. |    50    404 |
     |--------------|
  6. |    60    524 |
  7. |    70    709 |
  8. |    80    800 |
  9. |    90   1274 |
 10. |     .      . |
     +--------------+
*/

egen soma_50_pobres = sum (V4621) if V4621<=404
gen p_renda_50_pobres = soma_50_pobres/tot_renda

egen soma_10_ricos = sum (V4621) if V4621>=1274
gen p_renda_10_ricos = soma_10_ricos/tot_renda

sum p_renda_50_pobres p_renda_10_ricos



**** parada de pen

sort V4621
xtile decil_domic = V4621, nq(10) 
sort decil_domic
by decil_domic: egen renda_media_decil = mean(V4621)

tab renda_media_decil 

twoway (line renda_media_decil decil_domic, sort), ytitle(renda media) ylabel(, labels) xtitle(decimos da renda) xlabel(#10, labels) title(Parada de Pen)

twoway (bar renda_media_decil decil_domic, sort), ytitle(renda media) ylabel(, labels) xtitle(decimos da renda) xlabel(#10, labels) title(Parada de Pen)



*** curva de lorenz

twoway (line fi pi, sort) (line pi pi, sort), ytitle(renda acumulada) ylabel(, labels) xtitle(pop acumulada) xlabel(#10, labels) title(Curva de Lorenz)


**** Theil - T
sort V4621
egen dom_i_excl_zero = count(V4621) if V4621~=0
egen tot_domic_excl_zero = max(dom_i_excl_zero)

gen double termo_theil_T = parc_renda*log(parc_renda*tot_domic_excl_zero)
egen theil_T = sum(termo_theil_T)
tab theil_T

***** Theil - L

gen double termo_theil_L_1 = log(1/(parc_renda*tot_domic_excl_zero))
egen double termo_theil_L_2 = sum(termo_theil_L_1)
gen theil_L = (1/tot_domic_excl_zero)*termo_theil_L_2
tab theil_L


inequal7 V4621



** urbana ou rural

gen urbana =1 if V4105<=3
replace urbana =0 if V4105>3


inequal7 V4621 if urbana==1
inequal7 V4621 if urbana==0

*** quando calculamos os indices para as areas urbanas e rurais, estamos
*** calculando a desigualdade para cada grupo.... e não a desigualdade entre os grupos


*ssc inequal7