cap log close

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


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

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


/*
LINHAS de POBREZA
						alimentar		pobreza			pobreza_alta
North 			
Metro Belem 				130,55		200,21			371,81
Urban 						124,34		194,49			402,32
Rural 						122,26		177,33			280,29
			
Nordeste			
Metro Fortaleza 			122,26		188,77			379,44
Metro Recife 				128,48		198,30			400,41
Metro Salvador 				130,55		205,93			432,83
Urban 						124,34		190,67			345,12
Rural						122,26		175,42			266,94
			
Southeast 			
Metro Rio De Janeiro 		128,48		204,02			482,40
Metro Sao Paul0 			134,70		219,27			579,65
Metro Belo Horizonte 		122,26		196,39			465,24
Urban 						132,62		207,83			461,43
Rural						120,19		184,95			352,75
			
l5 south 			
Metro Curitiba 				124,34		200,21			501,47
Metro Port0 Alegre 			132,62		211,65			488,12
Urban 						118,12		188,77			427,11
Rural 						113,97		171,61			297,45
			
l9 Center West 			
Brasilia 					128,48		207,83			577,74
Goiania Municipality 		122,26		196,39			469,06
Urban 						126,41		200,21			453,80
Rural						124,34		190,67			364,19
			
National 						126,41		196,39		419,48
Metropolitan 					128,48		202,11		469,06
Urban excluding Metropolitan 	126,41		196,39		417,57
Rural 							120,19		179,23		312,70


*/


*** cálculo da proporção de pobres "H"

gen z = 345.12

gen pobre_vuln = 1 if V4621<=z
replace pobre_vuln = 0 if V4621>z & V4621~=.

egen n_pobres = sum(pobre_vuln)
egen total_n = count(pobre_vuln)
* a função count conta o numero de casos nao missing da vr
gen H = n_pobres/total_n

*** cálculo da razao do hiato da renda "I"


gen hiato = z - V4621 if pobre_vuln==1
replace hiato = 0 if pobre_vuln==0

egen soma_hiato = sum(hiato)

gen pz = n_pobres * z

gen I = soma_hiato / pz


**** cálculo do HI (indice de Foster, Greer e Thorbeck quando alfa =1 )

gen HI = H * I


**** cálculo do indice de Sen

sort V4621 
gen ind = _n if pobre_vuln==1
gen peso = n_pobres + 1 - ind if pobre_vuln==1
gen hiato_pond = hiato * peso
egen soma_hiato_pond = sum (hiato_pond)  
gen primeiro_termo = 2 / ((n_pobres +1) * total_n * z)

gen sen = primeiro_termo * soma_hiato_pond

**** cálculo do índice de Foster, Greer e Thorbeck quando alfa =2

gen hiato_2 = (z - V4621)^2 if pobre_vuln==1
replace hiato_2 = 0 if pobre_vuln==0

egen soma_hiato_2 = sum(hiato_2)

gen den = (total_n * (z^2))

gen FGT = soma_hiato_2 / den

****

sum H I HI sen FGT


*****

*** temos que instalar e dps usar


*apoverty varlist [weight] [if exp] [in exp] [, line(#) varpl(varname) gen(newvarname) all h h2 apg pgr igr w s tak thon fgt1 fgt2 fgt3 fgt4 fgt5 fgt6 fgt7 fgt8 fgt9 chu1 chu2 chu3 chu4 chu5]


apoverty V4621, varpl(z) gen(pobre_vuln_1) h igr pgr s fgt3   



**** decomposição

** urbana ou rural

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

egen n_pobres_urb1 = sum(pobre_vuln) if urbana==1
egen n_pobres_urb = mean(n_pobres_urb1)
drop n_pobres_urb1

egen total_n_urb1 = count(pobre_vuln) if urbana==1
egen total_n_urb = mean(total_n_urb1)
drop total_n_urb1

gen H_urb = n_pobres_urb/total_n_urb

egen n_pobres_rur1 = sum(pobre_vuln) if urbana==0
egen n_pobres_rur = mean(n_pobres_rur1)
drop n_pobres_rur1

egen total_n_rur1 = count(pobre_vuln) if urbana==0
egen total_n_rur = mean(total_n_rur1)
drop total_n_rur1

gen H_rur = n_pobres_rur/total_n_rur

gen pi_urb = total_n_urb / total_n
gen pi_rur = total_n_rur / total_n

gen H_alt = (pi_urb*H_urb) + (pi_rur*H_rur)

sum H H_alt H_urb H_rur



*** cálculo da razao do hiato da renda "I" separadamente para as areas urbana e rural

egen soma_hiato_urb1 = sum(hiato) if urbana==1
egen soma_hiato_urb = mean(soma_hiato_urb1)
drop soma_hiato_urb1

gen pz_urb = n_pobres_urb * z

gen I_urb = soma_hiato_urb / pz_urb


egen soma_hiato_rur1 = sum(hiato) if urbana==0
egen soma_hiato_rur = mean(soma_hiato_rur1)

gen pz_rur = n_pobres_rur * z

gen I_rur = soma_hiato_rur / pz_rur


gen I_alt = (pi_urb*I_urb) + (pi_rur*I_rur)

sum I I_alt I_rur I_urb



**** cálculo do HI (indice de Foster, Greer e Thorbeck quando alfa =1 )
*** separadamente para as areas urbana e rural

gen HI_urb = H_urb * I_urb
gen HI_rur = H_rur * I_rur

gen HI_alt = (pi_urb*HI_urb) + (pi_rur*HI_rur)


sum HI HI_alt HI_urb HI_rur


**************

**** cálculo do índice de Foster, Greer e Thorbeck quando alfa =2
*** para as areas urbana e rural

*urbana
egen soma_hiato_2_urb1 = sum(hiato_2) if urbana==1
egen soma_hiato_2_urb = mean(soma_hiato_2_urb1)
drop soma_hiato_2_urb1

gen den_urb = (total_n_urb * (z^2))

gen FGT_urb = soma_hiato_2_urb / den_urb

*rural
egen soma_hiato_2_rur1 = sum(hiato_2) if urbana==0
egen soma_hiato_2_rur = mean(soma_hiato_2_rur1)
drop soma_hiato_2_rur1

gen den_rur = (total_n_rur * (z^2))

gen FGT_rur = soma_hiato_2_rur / den_rur

gen FGT_alt = (pi_urb*FGT_urb) + (pi_rur*FGT_rur)

sum FGT FGT_alt FGT_rur FGT_urb



*povdeco5 varname [ [fweights aweights] if exp in range], pline(#) varpl(zvar) [bygroup(groupvar)]


povdeco5 V4621, varpl(z) by(urbana)



/*        subgroup poverty 'share', S_k = v_k.FGT_k(a)/FGT(a), and
        subgroup poverty 'risk',  R_k = FGT_k(a)/FGT(a) = S_k/v_k.
*/