/*
README: IN THESE TABLES WE USE IVIR8 (F.E. GRAVITY - version2) AS THE MAIN INSTRUMENT FOR THE MSH
*/

clear
set more off

cd "/Users/francesc/Desktop/revision_JIE
use data/work/maindata_130221.dta

log using log_131115.txt,replace

**********************************************************************************

ren pwt_openc tsh
ren fbsh msh
ren hk_shfb msh_hk

gen lnivtradesh2=ln(ivtradesh2)
** generate high openness variable
gen mostopen=(tsh>0.80)


gen msh_open=mostopen*msh
gen msh_nonopen=(1-mostopen) *msh

gen ivir3_open=mostopen*ivir3
gen ivir3_nonopen=(1-mostopen)*ivir3



gen ln_rely=ln(rel_y2000)
gen ln_relk=ln(rel_kya2000)
gen ln_relh=ln(rel_h2000)
gen ln_rela=ln(rel_a2000)

gen g70_00_rgdppop=(1/30)*ln(pwt_rgdppop/y1970)
gen lnrgdpop70=ln(y1970)

gen netrate=(msh-ertotal)


*** combining trade and immigration openness
** first principal component 
factor tsh msh, pcf factors(1)
predict princ_comp_open

factor ivir3 ivtradesh2, pcf factors(1)
predict princ_comp_inst


*************************************************************************************************
******** TABLE 1: DESCRIPTIVES
*************************************************************************************************

sum sample_fr pwt_rgdppop tsh msh pwt_ropen ertotal msh_hk iqual2 index_divm index_divt lnpop lnarea dis_equ tropicland euro1900 putt_migshare putt_statehist edes1975 gini d9_d1 ivtradesh1 ivtradesh2 ivtradesh3 ivtradesh9 ivir1 ivir2 ivir3 ivir9 
sum tsh ivtradesh*
sum msh ivir*

*scatter g_tsh tsh



*************************************************************************************************
******** TABLE 2: gravity regressions 
*************************************************************************************************

*(see outreg files in gravity dofile)





*************************************************************************************************
******** TABLE 3: TRADE OPENNESS AND INCOME PER PERSON
*************************************************************************************************

* 0) FR spec, linear pred as in FR, FR sample
*ivreg2 lngdppop (tsh = ivtradesh1) lnpop lnarea if sample_fr==1, robust ffirst /*FR sample income per person */
*outreg2 tsh lnpop lnarea dis_equ tropicland using table3.xls, dec(2) se bracket nolabel replace

* 1) FR spec, linear pred, FR sample
ivreg2 lngdppop (tsh = ivtradesh2) lnpop lnarea if sample_fr==1, robust ffirst /*FR sample income per person */
outreg2 tsh lnpop lnarea dis_equ tropicland using tables/table3.xls, dec(2) se bracket nolabel replace

*2) full sample
ivreg2 lngdppop (tsh = ivtradesh2) lnpop lnarea, robust	ffirst /*Full sample */
outreg2 tsh lnpop lnarea dis_equ tropicland using tables/table3.xls, dec(2) se bracket nolabel append

*3)   including distance to the equator
ivreg2 lngdppop (tsh = ivtradesh2) lnpop lnarea dis_equ if sample_fr==1, robust ffirst
outreg2 tsh lnpop lnarea dis_equ tropicland using tables/table3.xls, dec(2) se bracket nolabel append

*4) including All geography and colonial controls
ivreg2 lngdppop (tsh = ivtradesh2) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh lnpop lnarea dis_equ tropicland using tables/table3.xls, dec(2) se bracket nolabel append

**5) OLS
reg lngdppop tsh lnpop lnarea dis_equ d_* othercon1-othercon7 rich4 f_brit f_french if sample_fr==1, robust
outreg2 tsh lnpop lnarea dis_equ tropicland using tables/table3.xls, dec(2) se bracket nolabel append




*************************************************************************************************
******** TABLE 4: IMMIGRATION, TRADE AND INCOME
*************************************************************************************************

*1) FR sample
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) lnpop lnarea  if sample_fr==1, robust ffirst
outreg2 msh tsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel replace

*2) dist. equator
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) lnpop lnarea dis_equ if sample_fr==1, robust ffirst
outreg2 msh tsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel append

*3) Full sample
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) lnpop lnarea dis_equ, robust ffirst
outreg2 msh tsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel append

**4) FE predictor
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir8) lnpop lnarea dis_equ, robust ffirst
outreg2 msh tsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel append


**5 With controls for geo, resources, climate and colonial past
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 msh tsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel append
rivtest, ci

**6 REDUCED FORM REGRESSION (Direct regression)
regress lngdppop ivtradesh2 ivir3  lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1 , robust
outreg2 msh tsh lnpop lnarea dis_equ ivtradesh2 ivir3 using tables/table4.xls, dec(2) se bracket nolabel append


**7 OLS
regress lngdppop tsh msh  lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust
outreg2 msh tsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel append


**8 ONLY MSH 
ivreg2 lngdppop (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 msh tsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel append
rivtest, ci

*9 VALUE-ADDED TRADE SHARE and MSH
ivreg2 lngdppop (g_vatsh msh = ivtradesh2 ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 msh g_vatsh lnpop lnarea dis_equ using tables/table4.xls, dec(2) se bracket nolabel append


/*
*** CODE TO PRODUCE THE Angrist and Pishcke f-STAT for msh (for future reference)******************
local q=2
** first 
reg tsh ivtradesh2 ivir3 lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 f_brit f_french if sample_fr==1, robust
predict double tsh_1

reg msh tsh_1 lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 f_brit f_french if sample_fr==1, robust
predict double e, resid

** F-test for first stage of msh
reg e ivtradesh2 ivir3 lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 f_brit f_french if sample_fr==1, robust
testparm ivtradesh2 ivir3

dis r(F)*(`q'/(`q'-1))
** prob>F
dis 1-F(`q-1', 102, r(F)*(`q'/(`q'-1)))
*/




*************************************************************************************************
********* Table 5: ROBUSTNESS: MSH-HK, NET MIGRATION, HETEROGENEOUS EFFECTS
*************************************************************************************************


* 1: mig in hum cap, basic (no trade)
ivreg2 lngdppop (msh_hk = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel replace

* mig in hum cap, basic with institutions (no trade)
*ivreg2 lngdppop (msh_hk iqual2 = ivir3 edes1975 dis_equ) lnpop lnarea d_* tropicland othercon1-othercon7 oilres  if sample_fr==1, robust ffirst
*outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append

* 2 EMIGRATION
ivreg2  lngdppop (msh =  ivir3) ertotal lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres  rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append

* EMIGRATION and institutions
*ivreg2 lngdppop (msh iqual2 =  ivir3 edes1975 dis_equ) ertotal lnpop lnarea  d_* tropicland othercon1-othercon7 oilres  if sample_fr==1, robust ffirst
*outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append

* 3 NET IMMIGRATION
ivreg2  lngdppop (netrate =  ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres  rich4 f_brit f_french if sample_fr==1, robust  ffirst
outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append

* NET IMMIGRATION and institutions
*ivreg2 lngdppop (netrate iqual2 =  ivir3 edes1975 dis_equ) lnpop lnarea  d_* tropicland othercon1-othercon7 oilres   if sample_fr==1, robust ffirst
*outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append


********* Table B: HETEROGENEITY EFFECTS BY features of  RECEIVING COUNTRY

centile putt_statehist, centile(50)
gen high_statehist=0
replace high_statehist=1 if putt_statehist>.515

centile tsh, centile(50)
gen high_tsh=0
replace high_tsh=1 if tsh>0.80

centile ln_relh, centile(50)
gen high_relh=0
replace high_relh=1 if ln_relh>-0.4311

tab high_statehist
tab high_tsh
tab high_relh

gen msh_high_statehist = msh*high_statehist
gen msh_high_tsh = msh*high_tsh
gen msh_high_relh = msh*high_relh

gen msh_low_statehist = msh*(1-high_statehist)
gen msh_low_tsh = msh*(1-high_tsh)
gen msh_low_relh = msh*(1-high_relh)

gen ivir8_high_statehist = ivir8*high_statehist
gen ivir8_high_tsh = ivir8*high_tsh
gen ivir8_high_relh = ivir8*high_relh

gen ivir8_low_statehist = ivir8*(1-high_statehist)
gen ivir8_low_tsh = ivir8*(1-high_tsh)
gen ivir8_low_relh = ivir8*(1-high_relh)


gen ivir3_high_statehist = ivir3*high_statehist
gen ivir3_high_tsh = ivir3*high_tsh
gen ivir3_high_relh = ivir3*high_relh

gen ivir3_low_statehist = ivir3*(1-high_statehist)
gen ivir3_low_tsh = ivir3*(1-high_tsh)
gen ivir3_low_relh = ivir3*(1-high_relh)

/*
**4 BY STATEHIST
ivreg2 lngdppop (msh_high_statehist msh_low_statehist = ivir3_high_statehist ivir3_low_statehist) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append
test msh_high_statehist msh_low_statehist
*/

**4 BY OPENNESS TO TRADE
ivreg2 lngdppop (msh_high_tsh msh_low_tsh = ivir3_high_tsh ivir3_low_tsh) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french euro1900 if sample_fr==1, robust ffirst
outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append
test msh_high_tsh=msh_low_tsh

**5 BY RELH (AS DEFINED IN THE HALL AND JONES DECOMPOSITION, H/L)
ivreg2 lngdppop (msh_high_relh msh_low_relh = ivir3_high_relh ivir3_low_relh) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french euro1900 if sample_fr==1, robust ffirst
outreg2 princ_comp_open iqual2 msh* tsh ertotal netrate dis_equ using tables/table5.xls, dec(2) se bracket nolabel append
test msh_high_relh=msh_low_relh





*************************************************************************************************
********* Table 6: INSTITUTIONAL QUALITY
**** WHEN WE INCLUDE IQUAL AS REGRESSOR WE DO NOT INCLUDE THE COLONIAL ORIGIN CONTROLS
*************************************************************************************************




*1 - ENDOG INSTITUTIONS a la AC (euro descendents 1975)
ivreg2 lngdppop (msh iqual2 = ivir3 edes1975 dis_equ) lnpop lnarea  d_* othercon1-othercon7 oilres if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel replace

/*
***this is weaker and doesn't add anything so no point in reporting it.
** 2 FE PREDICTOR
ivreg2 lngdppop (msh iqual2 = ivir8 edes1975 dis_equ) lnpop lnarea d_* tropicland othercon1-othercon7 oilres  if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append
*/

* 2 CONTROL FOR DETERMINANTS OF INSTITUTIONS
ivreg2 lngdppop (msh = ivir3) lnpop lnarea dis_equ d_* edes1975 tropicland othercon1-othercon7 oilres if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append

*3 EVERYBODY NOW
ivreg2 lngdppop (tsh msh iqual2 = ivtradesh2 ivir3 edes1975 dis_equ) lnpop lnarea  d_* tropicland othercon1-othercon7 oilres if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append

*4 PW controls - no more controls
ivreg2 lngdppop (msh = ivir3) putt_statehist putt_migshare lnpop lnarea if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append

*5 PUTTERMAN-WEIL
ivreg2 lngdppop (msh iqual2 = ivir3 edes1975 dis_equ) putt_statehist putt_migshare lnpop lnarea  d_* tropicland othercon1-othercon7 oilres if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append

*6 PUTTERMAN-WEILL without institutions but with distance from equator
ivreg2 lngdppop (msh = ivir3) dis_equ putt_statehist putt_migshare lnpop lnarea d_* tropicland othercon1-othercon7 oilres  if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append

*7 IQUAL as dep. var.
ivreg2 iqual2 (msh = ivir3) putt_statehist putt_migshare lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 oilres  f_brit f_french if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append

*8 ALCALA AND CICCONE
ivreg2 lngdpemp (lnpwtropen iqual2 = lnivtradesh2 edes1975 dis_equ) lnpop lnarea if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append

*9 ALCALA AND CICCONE with controls and immigration share
ivreg2 lngdpemp (lnpwtropen msh iqual2 = lnivtradesh2 edes1975 dis_equ ivir3) lnpop lnarea  d_* tropicland othercon1-othercon7 oilres   if sample_fr==1, robust ffirst
outreg2 msh msh_hk iqual2 ace_avexpr ace_xcon1990sj lnpwtropen tsh dis_equ edes1975 putt_statehist putt_migshare using tables/table6.xls, dec(2) se bracket nolabel append




*************************************************************************************************
********* Table 7: Hall and Jones Decomposition.
*************************************************************************************************

** top panel: msh
ivreg2 ln_rely (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_frenc oilres if sample_fr==1, robust ffirst
outreg2 msh iqual2 lnpop lnarea using tables/table7.xls, dec(2) se bracket nolabel replace

ivreg2 ln_relk (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french oilres if sample_fr==1, robust ffirst
outreg2 msh iqual2 lnpop lnarea using tables/table7.xls, dec(2) se bracket nolabel append

ivreg2 ln_relh (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french oilres if sample_fr==1, robust ffirst
outreg2 msh  iqual2 lnpop lnarea using tables/table7.xls, dec(2) se bracket nolabel append

ivreg2 ln_rela (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french oilres if sample_fr==1, robust ffirst
outreg2 msh  iqual2 lnpop lnarea using tables/table7.xls, dec(2) se bracket nolabel append

*now with endogenous institutions
ivreg2 ln_rely (msh iqual2 = ivir3 edes1975 dis_equ) lnpop lnarea d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french oilres if sample_fr==1, robust ffirst
outreg2 msh  iqual2 lnpop lnarea using tables/table7.xls, dec(2) se bracket nolabel append

ivreg2 ln_rela (msh iqual2 = ivir3 edes1975 dis_equ) lnpop lnarea d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french oilres if sample_fr==1, robust ffirst
outreg2 msh  iqual2 lnpop lnarea using tables/table7.xls, dec(2) se bracket nolabel append


*********Inequality

replace gini=gini/100

** Gini \ MSH
ivreg2  gini (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 msh tsh princ_comp_open iqual2 using tables/table7.xls, dec(2) se bracket nolabel append

** 90-10 wage gap \ MSH TSH
ivreg2  d9_d1 (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 msh tsh princ_comp_open iqual2 using tables/table7.xls, dec(2) se bracket nolabel append



*************************************************************************************************
********* Table: 8 ethnic and linguistic FRACTIONALIZATION
*************************************************************************************************


**Exogenous diversity in trade and migration flows
**1 BASIC no diversity
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel replace

***2 Basic only migration diversity
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) index_divm  lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append

**3
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) index_divm index_divt lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append

**4
ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) index_divm3 index_divt2 lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append


**Fractionalization  and religion

ivreg2 frac_etn (tsh msh = ivtradesh2 ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append

ivreg2 frac_lang (tsh msh = ivtradesh2 ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append


*ivreg2 lngdppop (tsh msh = ivtradesh2 ivir3) frac_etn frac_lang catho80 muslim80 protmg80 lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
*outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append

ivreg2 lngdppop (msh = ivir3) frac_etn frac_lang catho80 muslim80 protmg80 lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append

*ivreg2 lngdppop (msh = ivir3) catho80 muslim80 protmg80 lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 oilres rich4 f_brit f_french if sample_fr==1, robust ffirst
*outreg2 tsh msh index_divm index_divt index_divm3 index_divt2 frac_etn frac_lang catho80 muslim80 protmg80 using tables/table8.xls, dec(2) se bracket nolabel append





*************************************************************************************************
**** TABLE 9: INNOVATION AND INEQUALITY
** INNOVATION MEASURED AS PATENTs PER MILLION ***** *** granted patents ****
*************************************************************************************************

*** patent variables ***

gen no_grant_per_million=1000*(no_grant/pwt_pop)
rename no_appl no_appl_per_million

gen lnpat_a=ln(no_appl_per_million)
gen lnpat_g=ln(no_grant_per_million)

gen pat_g=exp(lnpat_g)
gen levelpat=pat_g*pwt_pop
gen lnpat=ln(levelpat)

corr lnpat_a lnpat_g


**********Patenting

*1 Patents per cap \ MSH
ivreg2 pat_g (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 oilres f_brit f_french if sample_fr==1, robust
outreg2 msh tsh princ_comp_open iqual2 using tables/table9.xls, dec(2) se bracket nolabel replace

*2 Ln Patents per cap \ MSH
ivreg2 lnpat_g (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 oilres f_brit f_french if sample_fr==1, robust
outreg2 msh iqual2 tsh princ_comp_open iqual2 using tables/table9.xls, dec(2) se bracket nolabel append

*3 Ln Patents per cap \ MSH IQ
ivreg2 lnpat_g (msh iqual2 = ivir3 edes1975 dis_equ) lnpop lnarea  d_* othercon1-othercon7 oilres if sample_fr==1, robust
outreg2 msh iqual2 tsh  princ_comp_open iqual2 using tables/table9.xls, dec(2) se bracket nolabel append

*4 Ln Patents per cap \ TSH IQ
ivreg2 lnpat_g (tsh iqual2 = ivtradesh2 edes1975 dis_equ) lnpop lnarea  d_* othercon1-othercon7 oilres if sample_fr==1, robust
outreg2 msh iqual2 tsh  princ_comp_open iqual2 using tables/table9.xls, dec(2) se bracket nolabel append

*5 Ln Patents \ MSH
ivreg2 lnpat (msh = ivir3) lnpop lnarea dis_equ d_* tropicland othercon1-othercon7 rich4 oilres f_brit f_french if sample_fr==1, robust
outreg2 msh iqual2 tsh  princ_comp_open iqual2 using tables/table9.xls, dec(2) se bracket nolabel append

*6 Ln Patents \ MSH IQ
ivreg2 lnpat (msh iqual2 = ivir3 edes1975 dis_equ) lnpop lnarea  d_* othercon1-othercon7 oilres if sample_fr==1, robust
outreg2 msh iqual2 tsh  princ_comp_open iqual2 using tables/table9.xls, dec(2) se bracket nolabel append



*erase tables/table2_glin.txt
*erase tables/table2_gpml.txt
erase tables/table3.txt
erase tables/table4.txt
erase tables/table5.txt
erase tables/table6.txt
erase tables/table7.txt
erase tables/table8.txt
erase tables/table9.txt


******************THE MAGNITUDE OF THE EFFECT
drop if lngdppop==.
keep iso lngdppop msh iqual2
sum lngdppop, det
sort lngdppop

log close