-
Notifications
You must be signed in to change notification settings - Fork 0
/
Décomposition de la variance.do
238 lines (134 loc) · 7.77 KB
/
Décomposition de la variance.do
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
if "`c(username)'" =="guillaumedaudin" {
global dir_baseline_results "~/Documents/Recherche/2013 -- Trade Costs -- local/results/baseline"
global dir_referee1 "~/Documents/Recherche/2013 -- Trade Costs -- local/results/referee1"
global dir "~/Documents/Recherche/2013 -- Trade Costs -- local"
global dir_comparaison "~/Documents/Recherche/2013 -- Trade Costs -- local/results/comparaisons_various"
global dir_temp ~/Downloads/temp_stata
global dir_results "~/Documents/Recherche/2013 -- Trade Costs -- local/results"
global dir_redaction "~/Répertoires Git/trade_costs_git/redaction/JEGeo/revision_JEGeo/revised_article"
global dir_git "~/Répertoires Git/trade_costs_git/"
}
*** Juillet 2020: Lise, tout sur mon OneDrive
/* Fixe Lise P112*/
if "`c(hostname)'" =="LAB0271A" {
* baseline results sur hummels_tra dans son intégralité
global dir_baseline_results "C:\Users\lpatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results\baseline"
* résultats selon méthode référé 1
global dir_referee1 "C:\Users\lpatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results\referee1"
* stocker la comparaison des résultats
global dir_comparaison "C:\Users\lpatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results\referee1\comparaison_baseline_referee1"
/* Il me manque pour faire méthode 2 en IV
- IV_referee1_panel/results_estimTC_`year'_sitc2_3_`mode'.dta
- IV_ref1_y/results_estimTC_`year'_sitc2_3_`mode'.dta
*/
global dir_temp "C:\Users\lpatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\temp"
global dir "C:\Users\lpatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs"
global dir_results "C:\Users\lpatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results"
}
/* Nouveau portable Lise */
if "`c(hostname)'" =="MSOP112C" {
* baseline results sur hummels_tra dans son intégralité
global dir_baseline_results "C:\Users\Ipatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results\baseline"
* résultats selon méthode référé 1
global dir_referee1 "C:\Users\Ipatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results\referee1"
* stocker la comparaison des résultats
global dir_comparaison "C:\Users\Ipatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results\referee1\comparaison_baseline_referee1"
/* Il me manque pour faire méthode 2 en IV
- IV_referee1_panel/results_estimTC_`year'_sitc2_3_`mode'.dta
- IV_ref1_y/results_estimTC_`year'_sitc2_3_`mode'.dta
*/
global dir_temp "C:\Users\Ipatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\temp"
global dir "C:\Users\Ipatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs"
global dir_results "C:\Users\Ipatureau\OneDrive - Université Paris-Dauphine\Université Paris-Dauphine\trade_costs\results"
}
set more off
capture program drop decompo_var
program decompo_var
egen sd_tot=sd(prix_trsp), by(year mode)
gen var_tot=sd_tot^2
foreach dim in pays produit secteur {
if "`dim'"=="pays" egen sd_intra_`dim'=sd(prix_trsp), by(year mode iso_o)
if "`dim'"=="produit" egen sd_intra_`dim'=sd(prix_trsp), by(year mode sitc2)
if "`dim'"=="secteur" egen sd_intra_`dim'=sd(prix_trsp), by(year mode sector)
gen var_intra_`dim'=sd_intra_`dim'^2
egen var_intra_`dim'_moy=mean(var_intra_`dim'), by(year mode)
gen share_var_inter_`dim' = 1-var_intra_`dim'_moy/var_tot
}
bys year mode: keep if _n==1
foreach mode in air ves {
if "`mode'"== "air" local title_graph "Air"
if "`mode'"== "ves" local title_graph "Ves"
twoway (line share_var_inter_produit year) (line share_var_inter_secteur year) (line share_var_inter_pays year) /*(line product_country_variance year)*/ ///
if mode=="`mode'", ///
legend(label(1 "Share of between-product variance") label(3 "Share of between-country variance") label(2 "Share of between-sector variance") /*label(3 "Share of between product x country variance")*/ ///
rows(3)) ///
title("`mode'") name("`mode'", replace) ///
scheme(s1mono)
}
graph combine air ves, rows(1) scheme(s1mono)
end
use "$dir/data/hummels_tra.dta", clear
gen sector = substr(sitc2,1,3)
decompo_var
graph export "$dir_results/Decomposition_variance/Décomposition de la variance à la mimine_brut.png", replace
keep year mode sd_tot-share_var_inter_secteur
export delimited "$dir_results/Decomposition_variance/Décomposition de la variance à la mimine_brut.csv", replace
use "$dir/data/hummels_tra.dta", clear
gen sector = substr(sitc2,1,3)
bys sector: egen c_95_prix_trsp2 = pctile(prix_trsp2),p(95)
bys sector: egen c_05_prix_trsp2 = pctile(prix_trsp2),p(05)
drop if prix_trsp2 < c_05_prix_trsp2 | prix_trsp2 > c_95_prix_trsp2
decompo_var
graph export "$dir_results/Decomposition_variance/Décomposition de la variance à la mimine_ss_val_ext.png", replace
graph export "$dir_redaction/Décomposition de la variance à la mimine_ss_val_ext.png"
keep year mode sd_tot-share_var_inter_secteur
export delimited "$dir_results/Decomposition_variance/Décomposition de la variance à la mimine_ss_val_ext.csv", replace
/*Décidemment, je n'arrive pas à réconcilier les différentes méthodes.... :(
*******Test
use "$dir/data/hummels_tra.dta", clear
encode iso_o, gen(iso_o_num)
encode sitc2, gen(sitc2_num)
areg prix_trsp if mode=="ves" & year==1974, absorb(sitc2_num)
**Méthode xtreg
foreach mode in air ves {
use "$dir/data/hummels_tra.dta", clear
keep if mode=="`mode'"
encode iso_o, gen(iso_o_num)
encode sitc2, gen(sitc2_num)
gen iso_sitc = iso_o+sitc2
egen iso_sitc_num = group(iso_sitc)
xtreg prix_trsp if year==1974, i(sitc2_num)
local rho_sitc=e(rho)
xtreg prix_trsp if year==1974, i(iso_o_num)
local rho_iso=e(rho)
capture xtreg prix_trsp if year==1974, i(iso_sitc_num)
local rho_iso_sitc=e(rho)
matrix ana_var = (1974,`rho_sitc',`rho_iso',`rho_iso_sitc')
mat colnames ana_var = year product_variance country_variance product_country_variance
foreach y of num 1975(1)2013 {
xtreg prix_trsp if year==`y', i(sitc2_num)
local rho_sitc=e(rho)
xtreg prix_trsp if year==`y', i(iso_o_num)
local rho_iso=e(rho)
capture xtreg prix_trsp if year==`y', i(iso_sitc_num)
local rho_iso_sitc=e(rho)
matrix A = (`y',`rho_sitc',`rho_iso',`rho_iso_sitc')
matrix ana_var = ana_var \ A
}
matrix list ana_var
clear
svmat ana_var, names (col)
twoway (line product_variance year) (line country_variance year) /*(line product_country_variance year)*/, ///
legend(label(1 "Share of between product variance") label(2 "Share of between country variance") /*label(3 "Share of between product x country variance")*/ ///
rows(2)) ///
title("`mode'") name("`mode'", replace)
gen mode = "`mode'"
save "$dir/results/Décomposition de la variance/Décomposition de la variance à la xtreg_`mode'.dta", replace
}
use "$dir/results/Décomposition de la variance/Décomposition de la variance à la xtreg_air.dta", clear
append using "$dir/results/Décomposition de la variance/Décomposition de la variance à la xtreg_ves.dta", replace
erase "$dir/results/Décomposition de la variance/Décomposition de la variance à la xtreg_air.dta"
erase "$dir/results/Décomposition de la variance/Décomposition de la variance à la xtreg_ves.dta"
graph combine air ves, rows(2)
graph export "$dir/results/Décomposition de la variance/Décomposition de la variance à la xtreg.png", replace
export delimited "$dir/results/Décomposition de la variance/Décomposition de la variance à la xtreg.csv", replace