*This files replicates the results given in the online appendix with dummies for every 200 years, Table E7 *Remember to change the paths as appropriate use "M:\Nov2014\Plov\Resub_JDE_CVS\Stata\ReplicationFiles\Europe3.dta", clear gen _Iyear__600_700 = _Iyear_600 + _Iyear_700 gen _Iyear__800_900 = _Iyear_800 + _Iyear_900 gen _Iyear__1000_1100 = _Iyear_1000 + _Iyear_1100 gen _Iyear__1200_1300 = _Iyear_1200 + _Iyear_1300 *Gen interactions: forvalues k=600(200)1300{ gen Gleyqualitywheat55_`k'= _Iyear__`k'*(qualitybarleyshareover55+GoodQGleyShare) } foreach VAR in CoastDistKm RomanHeritage Equator{ forvalues k=600(200)1300{ gen `VAR'_`k'= _Iyear__`k'*`VAR' } } forvalues k=600(200)1300{ gen qualitywheat55luv`k'= _Iyear__`k'*qualitywheatshareover55 } forvalues k=600(200)1300{ gen qualitybarley40luv`k'= _Iyear__`k'*qualitybarleyshareover40 } forvalues k=600(200)1300{ gen qualitybarley55luv`k'= _Iyear__`k'*qualitybarleyshareover55 } forvalues k=600(200)1300{ gen qualitybarley70luv`k'= _Iyear__`k'*qualitybarleyshareover70 } forvalues k=600(200)1300{ gen ln_CoastDistKm`k'= ln(1+CoastDistKm)*_Iyear__`k' } forvalues k=600(200)1300{ gen RomanRobustHeritage_`k'= _Iyear__`k'*RomanHeritageRobust } forvalues k=600(200)1300{ gen RiverDist_`k'= _Iyear__`k'*AvgDistRiverNew } forvalues k=600(200)1300{ gen IronD_`k'= _Iyear__`k'*IronDummy } gen ln_Avg_MT = ln(Avg_MT) gen qualityluvisolwheat55_I_1000 = qualitywheatshareover55*i_1000 gen qualityluvisolbarley55_I_1000 = qualitybarleyshareover55*i_1000 gen ln_qualityluvisolbarley55_I_1000 =ln(1+qualityluvisolbarley55_I_1000 ) gen ln_CityCountDens = ln(1+CityCountDens) gen Avg_MT_sq = Avg_MT^2 ************************************************** ********************* Table E7 ******************* ************************************************** *Column 1: Urbanization, NonFlex, No controls: xi: areg CityCountDens qualityluvisolbarley55_I_1000 i.year if nutslevel ==2 & undefinedsoilshare < 0.2 & year < 1400, a(nuts_code) cluster(nuts_code) outreg2 using M:\Nov2014\Plov\Resub_JDE_CVS\Stata\T3, excel dec(5) adjr2 append *Column 2: Urbanization, NonFlex, Controls: xi: areg CityCountDens qualityluvisolbarley55_I_1000 RomanHeritage_* Unis CoastDistKm_* Avg_MT i.year if nutslevel ==2 & undefinedsoilshare < 0.2 & year < 1400, a(nuts_code) cluster(nuts_code) outreg2 using M:\Nov2014\Plov\Resub_JDE_CVS\Stata\T3, excel dec(5) adjr2 append *Column 3: Urbanization, Flex, No controls: xi: areg CityCountDens qualitybarley55* i.year if nutslevel ==2 & undefinedsoilshare < 0.2 & year < 1400, a(nuts_code) cluster(nuts_code) outreg2 using M:\Nov2014\Plov\Resub_JDE_CVS\Stata\Years_200, excel dec(5) adjr2 append *Column 4: Urbanization, Flex, Controls: xi: areg CityCountDens qualitybarley55* RomanHeritage_* Unis CoastDistKm_* Avg_MT i.year if nutslevel ==2 & undefinedsoilshare < 0.2 & year < 1400,a(nuts_code) cluster(nuts_code) outreg2 using M:\Nov2014\Plov\Resub_JDE_CVS\Stata\Years_200, excel dec(5) adjr2 append