/************ 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 STI8*/ /****sélection des variables*/ use Resultats/TousImpots/Temp/temp_tous_impots_$annee, replace keep if age>=18 & age<=122 egen masse_cs_contr=total(cs_contr*pondv/1000000000) egen masse_yr_prim=total(yr_prim*pondv/1000000000) egen masse_yr_prima=total(yr_prima*pondv/1000000000) gen ik=is+tf+isf+dmtg gen ika=is+tf+isf-bouclier+dmtg replace irf=irf+ir_remp*masse_cs_contr/masse_yr_prim replace iraf=iraf+ira_remp*masse_cs_contr/masse_yr_prima replace tvaf=tvaf+tva_conso_remp*masse_cs_contr/masse_yr_prim replace tpf=tpf+tp_conso_remp*masse_cs_contr/masse_yr_prim replace thf=thf+th_remp*masse_cs_contr/masse_yr_prim gen iif=tvaf+tpf+thf replace tvaaf=tvaaf+tva_consoa_remp*masse_cs_contr/masse_yr_prima replace tpaf=tpaf+tp_consoa_remp*masse_cs_contr/masse_yr_prima gen iiaf=tvaaf+tpaf+thf gen csts=cs+ts-cs_contr+csk replace yr_prim=yr_prim*masse_cs_contr/masse_yr_prim replace yr_prima=yr_prima*masse_cs_contr/masse_yr_prima replace ya_prim=ya_prim-cs_contr replace ya_prima=ya_prima-cs_contr gen rfin_prim=rfin_cn*yk_prim/(yk_cn+0.00001) gen rfon_prim=rfon_cn*yk_prim/(yk_cn+0.00001) gen sal_prim=sal_superbrut_cn*ya_prim/(ya_cn+0.00001) gen nonsal_prim=nonsal_brut_cn*ya_prim/(ya_cn+0.00001) gen rfin_prima=rfin_cn*yk_prima/(yk_cn+0.00001) gen rfon_prima=rfon_cn*yk_prima/(yk_cn+0.00001) gen sal_prima=sal_superbrut_cn*ya_prima/(ya_cn+0.00001) gen nonsal_prima=nonsal_brut_cn*ya_prima/(ya_cn+0.00001) keep pondv y_sec impot_sec irf ik ika iif csts is tf isf bouclier isf2 dmtg tvaf tpf thf csts yk_prim rfin_prim rfon_prim ya_prim sal_prim nonsal_prim yr_prim y_seca impota_sec iraf iiaf tvaaf tpaf yk_prima rfin_prima rfon_prima ya_prima sal_prima nonsal_prima yr_prima save Resultats/TousImpots/Temp/temp3_tous_impots_$annee, replace /*****fractiles très riches*/ use Resultats/TousImpots/Temp/temp3_tous_impots_$annee, clear xtile py=y_sec [w=pondv], nq(100) keep if py==100 drop py xtile py=y_sec [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_sec (mean) y_sec impot_sec irf ik iif csts is tf isf bouclier isf2 dmtg tvaf tpf thf yk_prim rfin_prim rfon_prim ya_prim sal_prim nonsal_prim yr_prim [w=pondv], by(py) save Resultats/TousImpots/TableSTI8_$annee, replace /*****fractiles moyens riches*/ use Resultats/TousImpots/Temp/temp3_tous_impots_$annee, clear xtile py=y_sec [w=pondv], nq(100) replace py=py-1 keep if py>=95 & py<=98 collapse (rawsum) pondv (min) seuil=y_sec (mean) y_sec impot_sec irf ik iif csts is tf isf bouclier isf2 dmtg tvaf tpf thf yk_prim rfin_prim rfon_prim ya_prim sal_prim nonsal_prim yr_prim [w=pondv], by(py) append using Resultats/TousImpots/TableSTI8_$annee save Resultats/TousImpots/TableSTI8_$annee, replace /*****fractiles déciles*/ use Resultats/TousImpots/Temp/temp3_tous_impots_$annee, clear xtile py=y_sec [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_sec (mean) y_sec impot_sec irf ik iif csts is tf isf bouclier isf2 dmtg tvaf tpf thf yk_prim rfin_prim rfon_prim ya_prim sal_prim nonsal_prim yr_prim [w=pondv], by(py) append using Resultats/TousImpots/TableSTI8_$annee save Resultats/TousImpots/TableSTI8_$annee, replace /*****fractiles pauvres-moyens-riches*/ use Resultats/TousImpots/Temp/temp3_tous_impots_$annee, clear xtile py=y_sec [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_sec (mean) y_sec impot_sec irf ik iif csts is tf isf bouclier isf2 dmtg tvaf tpf thf yk_prim rfin_prim rfon_prim ya_prim sal_prim nonsal_prim yr_prim [w=pondv], by(py) append using Resultats/TousImpots/TableSTI8_$annee save Resultats/TousImpots/TableSTI8_$annee, replace /*****moyennes*/ use Resultats/TousImpots/Temp/temp3_tous_impots_$annee, clear collapse (rawsum) pondv (min) seuil=y_sec (mean) y_sec impot_sec irf ik iif csts is tf isf bouclier isf2 dmtg tvaf tpf thf yk_prim rfin_prim rfon_prim ya_prim sal_prim nonsal_prim yr_prim [w=pondv] append using Resultats/TousImpots/TableSTI8_$annee replace pondv=pondv/1000 replace impot_sec=impot_sec/y_sec replace irf=irf/y_sec replace ik=ik/y_sec replace iif=iif/y_sec replace csts=csts/y_sec replace is=is/y_sec replace tf=tf/y_sec replace isf=isf/y_sec replace bouclier=bouclier/y_sec replace isf2=isf2/y_sec replace dmtg=dmtg/y_sec replace tvaf=tvaf/y_sec replace tpf=tpf/y_sec replace thf=thf/y_sec replace yk_prim=yk_prim/y_sec replace rfin_prim=rfin_prim/y_sec replace rfon_prim=rfon_prim/y_sec replace ya_prim=ya_prim/y_sec replace sal_prim=sal_prim/y_sec replace nonsal_prim=nonsal_prim/y_sec replace yr_prim=yr_prim/y_sec save Resultats/TousImpots/TableSTI8_$annee, replace