cd "C:\Users\alejo\Dropbox\Teaching\EC2009\JEcEd\Data\"
*cd "C:\Users\sbrp073\Dropbox\Teaching\EC2009\JEcEd\Data"

capture log close
log using JEEcEd.txt, replace

use Colombia-2017-full-data.dta, replace

* Keeping only manufacturing firms
keep if a0==1

* Determining how many manufacturing firms are there in the data
tab a0

*** Firm heterogeneity: examples

gen log10emp = log10(l1)
twoway (kdensity log10emp if a4b==15) (kdensity log10emp if a4b==17 | a4b==18 | a4b==19) (kdensity log10emp if a4b==24 | a4b==25)


*** Indicators of global engagement: exporting, importing, being foreign-owned

* Create a dummy variable called 'export' = 1 if an establishment exports some of its output (directly or indirectly) and 0 otherwise
gen pctexp = d3b + d3c
gen export=.
replace export=0 
replace export=1 if pctexp>0 

* Create a dummy variable called 'import' = 1 if an establishment imports some of its inputs and 0 otherwise
gen import=.
replace import=1 if d12b>0 & d12b~=.
replace import=0 if d12b==0 & d12b~=-9

* Create a dummy variable called 'foreign' = 1 if the % of establishment owned by foreigners is at least 10% and 0 otherwise
gen b2 = b2a + b2b + b2c + b2d
gen foreign=.
replace foreign=1 if b2b>=10
replace foreign=0 if b2b<10 & b2~=.

tab export
tab import
tab foreign

tab export import
tab export foreign
tab import foreign

* Calculate observations per industry
egen obsind = count(id), by(a4b)

* plotting incidence of exporting, importing and beign foreign-owned in sectors with more than 20 establishments:
graph hbar (mean) export import foreign if obsind>=20, over(a4b)

* Plotting export intensity distribution:
hist pctexp if pctexp>0, bin(10) frac xtitle(export intensity)


*** Constructing performance indicators:
* log employment
gen logemp = log(l1)

* log skill-intensity (share of permanent non-production workers in total permanent workers)
gen skillint = l3b/(l3a + l3b)
replace skillint=. if skillint<0 | skillint>1
gen logskillint = log(skillint)

* log sales
gen logsales = log(d2)

* log sales per worker
gen salespw = d2/l1
replace salespw=. if salespw<0
gen logsales_pw = log(d2/l1) 

* log capital per worker
gen kpw = n6a/l1
gen logk_pw = log(kpw)

* log female employment
gen femempint = (l5a + l5b)/l1
replace femempint=. if femempint<0 | femempint>1
gen logfemempint = log(femempint)

* log R&D expenditure/sales
gen rd = h9
replace rd=. if rd<0
gen logrd = log(rd)

* export regressions
reg logemp export,r
reg logsales export, r
reg logsales_pw export,r
reg logk_pw export,r
reg logskillint export,r
reg logfemempint export,r
reg logrd export,r

* import regressions
reg logemp import,r
reg logsales import, r
reg logsales_pw import,r
reg logk_pw import,r
reg logskillint import,r
reg logfemempint import,r
reg logrd import,r

* foreign-owned regressions
reg logemp foreign,r
reg logsales foreign, r
reg logsales_pw foreign,r
reg logk_pw foreign,r
reg logskillint foreign,r
reg logfemempint foreign,r
reg logrd foreign,r