/************ simulation des taux globaux de prélèvements: calcul de statistiques descriptives*/ clear set mem 500m set more off cd "C:\piketty\RevolutionFiscale2010\Simulateur\" /******** appel des paramètres*/ do Programmes/do_appel_parametres /*************préparation du tableau STI5*/ /****sélection des variables*/ use Resultats/TousImpots/Temp/temp_tous_impots_$annee, replace keep if age>=18 & age<=65 keep if ya_prim>0 & tempspartiel==0 gen ik=is+tf+isf+dmtg gen ika=is+tf+isf-bouclier+dmtg gen iif=tvaf+tpf+thf gen iiaf=tvaaf+tpaf+thf gen csts=cs+ts+csk keep pondv y_prim impot_prim irf ik ika iif csts is tf isf bouclier isf2 dmtg tvaf tpf thf cs ts cs_contr csk yk_prim ya_prim y_prima impota_prim iraf iiaf tvaaf tpaf yk_prima ya_prima save Resultats/TousImpots/Temp/temp2_tous_impots_$annee, replace /*****fractiles très riches*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=y_prima [w=pondv], nq(100) keep if py==100 drop py xtile py=y_prima [w=pondv], nq(1000) replace py=py-1 replace py=0 if py>=0 & py<900 replace py=900 if py>=900 & py<990 replace py=990 if py>=990 & py<999 replace py=999 if py>=999 collapse (rawsum) pondv (min) seuil=y_prima (mean) y_prima impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk yk_prima ya_prima [w=pondv], by(py) save Resultats/TousImpots/TableSTI5_$annee, replace /*****fractiles moyens riches*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=y_prima [w=pondv], nq(100) replace py=py-1 keep if py>=95 & py<=98 collapse (rawsum) pondv (min) seuil=y_prima (mean) y_prima impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk yk_prima ya_prima [w=pondv], by(py) append using Resultats/TousImpots/TableSTI5_$annee save Resultats/TousImpots/TableSTI5_$annee, replace /*****fractiles déciles*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=y_prima [w=pondv], nq(100) replace py=py-1 replace py=10*int(py/10) if py<95 replace py=95 if py>=95 & py<99 replace py=99 if py>=99 collapse (rawsum) pondv (min) seuil=y_prima (mean) y_prima impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk yk_prima ya_prima [w=pondv], by(py) append using Resultats/TousImpots/TableSTI5_$annee save Resultats/TousImpots/TableSTI5_$annee, replace /*****fractiles pauvres-moyens-riches*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=y_prima [w=pondv], nq(100) replace py=py-1 replace py=0 if py<50 replace py=50 if py>=50 & py<90 replace py=90 if py>=90 collapse (rawsum) pondv (min) seuil=y_prima (mean) y_prima impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk yk_prima ya_prima [w=pondv], by(py) append using Resultats/TousImpots/TableSTI5_$annee save Resultats/TousImpots/TableSTI5_$annee, replace /*****moyennes*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear collapse (rawsum) pondv (min) seuil=y_prima (mean) y_prima impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk yk_prima ya_prima [w=pondv] append using Resultats/TousImpots/TableSTI5_$annee replace pondv=pondv/1000 replace impota_prim=impota_prim/y_prima replace iraf=iraf/y_prima replace ika=ika/y_prima replace iiaf=iiaf/y_prima replace csts=csts/y_prima replace is=is/y_prima replace tf=tf/y_prima replace isf=isf/y_prima replace bouclier=bouclier/y_prima replace isf2=isf2/y_prima replace dmtg=dmtg/y_prima replace tvaaf=tvaaf/y_prima replace tpaf=tpaf/y_prima replace thf=thf/y_prima replace cs=cs/y_prima replace ts=ts/y_prima replace cs_cont=cs_cont/y_prima replace csk=csk/y_prima replace yk_prima=yk_prima/y_prima replace ya_prima=ya_prima/y_prima save Resultats/TousImpots/TableSTI5_$annee, replace