/************ 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 STI9*/ /****sélection des variables*/ use Resultats/TousImpots/Temp/temp_tous_impots_$annee, replace keep if age>=18 & age<=122 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 k_cn* k_isf replace k_cn=0.001*uniform() if k_cn==0 save Resultats/TousImpots/Temp/temp2_tous_impots_$annee, replace /*****fractiles très riches*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=k_cn [w=pondv], nq(100) keep if py==100 drop py xtile py=k_cn [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=k_cn (mean) k_cn k_isf impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk [w=pondv], by(py) save Resultats/TousImpots/TableSTI9_$annee, replace /*****fractiles moyens riches*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=k_cn [w=pondv], nq(100) replace py=py-1 keep if py>=95 & py<=98 collapse (rawsum) pondv (min) seuil=k_cn (mean) k_cn k_isf impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk [w=pondv], by(py) append using Resultats/TousImpots/TableSTI9_$annee save Resultats/TousImpots/TableSTI9_$annee, replace /*****fractiles déciles*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=k_cn [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=k_cn (mean) k_cn k_isf impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk [w=pondv], by(py) append using Resultats/TousImpots/TableSTI9_$annee save Resultats/TousImpots/TableSTI9_$annee, replace /*****fractiles pauvres-moyens-riches*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear xtile py=k_cn [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=k_cn (mean) k_cn k_isf impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk [w=pondv], by(py) append using Resultats/TousImpots/TableSTI9_$annee save Resultats/TousImpots/TableSTI9_$annee, replace /*****moyennes*/ use Resultats/TousImpots/Temp/temp2_tous_impots_$annee, clear collapse (rawsum) pondv (min) seuil=k_cn (mean) k_cn k_isf impota_prim iraf ika iiaf csts is tf isf bouclier isf2 dmtg tvaaf tpaf thf cs ts cs_contr csk [w=pondv] append using Resultats/TousImpots/TableSTI9_$annee replace pondv=pondv/1000 replace k_isf=k_isf/k_cn replace impota_prim=impota_prim/k_cn replace iraf=iraf/k_cn replace ika=ika/k_cn replace iiaf=iiaf/k_cn replace csts=csts/k_cn replace is=is/k_cn replace tf=tf/k_cn replace isf=isf/k_cn replace bouclier=bouclier/k_cn replace isf2=isf2/k_cn replace dmtg=dmtg/k_cn replace tvaaf=tvaaf/k_cn replace tpaf=tpaf/k_cn replace thf=thf/k_cn replace cs=cs/k_cn replace ts=ts/k_cn replace cs_cont=cs_cont/k_cn replace csk=csk/k_cn save Resultats/TousImpots/TableSTI9_$annee, replace