-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathTableau comparaison referee1 & qy.do
338 lines (216 loc) · 9.86 KB
/
Tableau comparaison referee1 & qy.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
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
*ssc inst _gwtmean
*Programme fait à partir de "Comparaison.do
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_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_referee1_yearly/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_referee1_yearly/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
/*
************Comparaison de base
use "$dir/data/hummels_tra.dta", clear
contract year iso_o sitc2 mode
tab _freq
/*Cela confirme que year iso_o sitc2 mode sont les clefs du fichier*/
use "$dir/data/base_hs10_newyears.dta", clear
contract year iso_o mode hs dist_entry
tab _freq
*Ce n’est pas une clef unique ?! Donc il y a plusieurs consignements par produit/district dans base_hs10_newyears ?
use "$dir/data/hummels_tra.dta", clear
keep if year >=2005
merge 1:m year iso_o sitc2 mode using "$dir/data/base_hs10_newyears.dta", force
tab mode _merge
**Semble suggérer qu’il y a plus de choses dans HS10 que dans hummels_tra... (même au delà des cnt)
drop if sitc2==""
tab mode _merge
drop if mode=="cnt"
tab _merge
*****Donc tout le problème est bien lié à mode=="cnt" et aux sitc vides
generate sector = substr(sitc2,1,3)
codebook sector
contract year iso_o sector mode
describe
******************************************
***************** Pour vérifier que le merge se fait sur les bases d’orgine... c’est bon aussi
use "$dir/data/hummels_tra.dta", clear
rename sitc2 sector
drop if sector==""
assert strlen(sector)==5
replace sector = substr(sector,1,3)
keep if year >=2005
contract year iso_o sector mode
save temp_hummels_tra.dta, replace
use "$dir/data/base_hs10_newyears.dta", clear
rename sitc2 sector
drop if sector==""
assert strlen(sector)==5
replace sector = substr(sector,1,3)
drop if mode=="cnt"
contract year iso_o sector mode
merge 1:1 year sector iso_o mode using temp_hummels_tra.dta, force
erase temp_hummels_tra.dta
***********************
*/
** Faire tourner sur toutes les années /mode
/* method peut être
- "baseline" (nos benchmark results en s=3 digits, k=5 digits)
- "baseline10" (nos benchmark results en s=3 digits, k=10 digits)
- "baselinesamplereferee1" = notre methode sur le sample issu de la méthode du référé 1, en s=3, k=5 ou 10 (A ACTUALISER)
- "referee1" (methode OLS référé 1), s=3 k=10
- "IV_referee1_panel" (??)
- "IV_referee1_yearly" (??)
*/
******************************************************
******************************************************
do "$dir_git/Open_year_mode_method_model.do"
*****************************************************************************************
***on lance les programmes
*****************************************************************************************
*global method baseline
*global method referee1
******
******************Pour la table de comparaison
collect clear
capture program drop table_comparaison_part
program table_comparaison_part
args method
global method `method'
local time_span $time_span
local model $model
foreach mode in air ves {
capture erase $dir_temp/data_`model'_${method}_`mode'.dta
foreach year of num `time_span' {
open_year_mode_method_model `year' `mode' $method `model'
capture append using $dir_temp/data_`model'_${method}_`mode'.dta
save $dir_temp/data_`model'_${method}_`mode'.dta, replace
}
use $dir_temp/data_`model'_${method}_`mode'.dta, replace
drop if beta==.
egen value_year=total(val), by(year)
generate weight = val/value_year
generate N = 0
generate Nb_sectors = 0
generate Nb_partners = 0
generate Nb_pairs = 0
gen paire =iso_o+sector
gen value_tot = 0
foreach year of num `time_span' {
capture levelsof iso_o if year==`year'
replace Nb_partners=r(r) if year==`year'
capture levelsof sector if year==`year'
replace Nb_sectors=r(r) if year==`year'
capture levelsof paire if year==`year'
replace Nb_pairs=r(r) if year==`year'
egen N_`year'=count(beta), by(year)
replace N=N_`year' if year==`year'
egen value_`year'=total(val), by(year)
replace value_tot=value_`year' if year==`year'
drop N_`year'
}
collect, tags(var[N] mode[`mode'] digit[${method}]): /*
*/ sum N [aweight=weight]
collect, tags(var[Nb_sectors] mode[`mode'] digit[${method}]): /*
*/ sum Nb_sectors [aweight=weight]
collect, tags(var[Nb_partners] mode[`mode'] digit[${method}]): /*
*/ sum Nb_partners[aweight=weight]
collect, tags(var[Nb_pairs] mode[`mode'] digit[${method}]): /*
*/ sum Nb_pairs[aweight=weight]
collect, tags(var[value_tot] mode[`mode'] digit[${method}]): /*
*/ sum value_tot[aweight=weight]
collect, tags(var[beta] mode[`mode'] digit[${method}]) /*
*/ : sum beta[aweight=weight], det
save $dir_temp/data_`model'_${method}_`mode'.dta, replace
collect, tags(var[beta_${debut}] mode[`mode'] digit[${method}]) /*
*/ : sum beta[aweight=weight] if year==${debut}, det
collect, tags(var[beta_${fin}] mode[`mode'] digit[${method}]) /*
*/ : sum beta[aweight=weight] if year==${fin}, det
gen beta_pour_collapse = beta*weight
collapse (sum) beta_pour_collapse, by(year)
quietly gen ln_beta=ln(beta_pour_collapse)
collect _r_b, tags(mode[`mode'] digit[${method}]): regress ln_beta year /*[aweight=weight]*/
*if "`mode'"=="air" & "${method}"=="non_séparé_wgt" blink
collect layout (result[mean]#var[Nb_sectors Nb_partners Nb_pairs value_tot] /*
*/ var[beta]#result[mean p50 sd] (colname[year]#result) (var[beta_$debut beta_${fin}]#result[mean])) /*
*/ (mode#digit)
}
collect style cell, warn nformat (%3.1f)
collect style cell var[beta beta_$debut beta_$fin], warn nformat(%3.2f)
collect style cell var[Nb_pairs]#var[value_tot]#var[Nb_sectors]#var[Nb_partners], warn nformat(%9.0fc)
collect style cell colname[year], warn nformat(%4.3f)
collect style column, nodelimiter dups(center) position(top) width(asis)
collect label levels colname year "Time trend", modify
collect label levels result p50 "Median", modify
collect label levels var value_tot "Covered trade value", modify
collect label levels digit referee1 "Estimating $\beta_{i,s}$", modify
collect label levels digit non_séparé_qy "Price per qy", modify
collect label levels digit non_séparé_wgt "Price per kg", modify
end
**********Pour comparer baseline / beta referee1
global time_span 2005 (1) 2013
global debut 2005
global fin 2013
global model nlAetI
table_comparaison_part referee1
table_comparaison_part baseline
collect preview
collect export /*
*/ "$dir_git/redaction/JEGeo/revision_JEGeo/revised_article/Online_Appendix/Comp_baseline_referee1.tex", /*
*/ tableonly replace
collect clear
**********
********************Pour comparer baseline / wgt / qy
*
global time_span 2009 (1) 2019
global debut 2009
global fin 2019
global model nlAetI
/**Mais en fait, ce ne sont pas les bonnes régressions, car non-séparées
table_comparaison_part hs10_qy1_qy
table_comparaison_part hs10_qy1_wgt
table_comparaison_part baseline
*/
table_comparaison_part non_séparé_qy
table_comparaison_part non_séparé_wgt
table_comparaison_part baseline
collect preview
collect export /*
*/ "$dir_git/redaction/JEGeo/revision_JEGeo/revised_article/Online_Appendix/Comp_baseline_wgt_qy.tex", /*
*/ tableonly replace
collect clear