set more off cap log close log using results_doda_et_al_2016, replace text *use full_dataset_with_5DVs.dta, replace use data_doda_et_al_2016.dta,replace gene lnemp=log(employees) list organisationname icbsectordescription year d_p d_s d_t if s_5==1 replace d_t=0 if s_5==1 replace d_p=1 if s_5==1 list organisationname icbsectordescription year d_p d_s d_t if s_5==1 list organisationname icbsectordescription year d_p d_s d_t if s_9==1 replace d_t=0 if s_9==1 replace d_s=1 if s_9==1 list organisationname icbsectordescription year d_p d_s d_t if s_9==1 list organisationname icbsectordescription year d_p d_s d_t if s_10==1 replace d_t=0 if s_10==1 replace d_s=1 if s_10==1 list organisationname icbsectordescription year d_p d_s d_t if s_10==1 list organisationname icbsectordescription year d_p d_s d_t if s_3==1 replace d_t=0 if s_3==1 replace d_s=1 if s_3==1 list organisationname icbsectordescription year d_p d_s d_t if s_3==1 * description of the dataset sort year id by year: sum lnco2 lni* lnoprev d_* sum p0*_d m0*_d dm0*_d d0*_d d_* y_* ************************************************* ***generate dummies for JAPAN and Australia******************** ************************************************* gen d_JA=0 replace d_JA=1 if countrycode =="JA" replace d_eu27=0 if countrycode =="AU" gen d_AU=0 replace d_AU=1 if countrycode =="AU" sort id year ************************************************ *************Outliers********* **************************************** drop if organisationname=="OneSteel" & year==2007 ********************************************************** *******generate growth rate of the log Co2/oprev********** ********************************************************** by id: gen ln_CO2gr=(lni_co2oprev-l.lni_co2oprev) by id: gen CO2gr=(i_co2oprev-l.i_co2oprev)/l.i_co2oprev by id: gen EMPLgr=(employees-l.employees)/l.employees by id: gen ln_EMPLgr=(lnemp-l.lnemp) ***Outliers***** drop if CO2gr>9 & CO2gr!=. ***Descriptive statistics**** sum lnco2 lnoprev lni_co2oprev ln_CO2gr lnemp p0*_s m0*_s dm0*_s d0*_s hist lni_co2oprev, normal bin(30) graph save CO2Distribution,replace gen count=1 preserve sort year by year: egen TOTAL=total(count) tab year TOTAL graph bar TOTAL, over(year) graph save numberfirms,replace egen firmsbysector = total(count), by (icbsectordescription) label var firmsbysector "Count of reporting firms" sort firmsbysector graph hbar firmsbysector, over(icbsectordescription, sort(1) descending) graph save firmsbysector,replace restore keep if year==2010 ***Table 4: Continuous variables for CMPs************** **** Column (I): All CMPs and Firms *********************** ******************************************************** *preserve regress ln_CO2gr /// p0*_s m0*_s dm0*_s d0*_s /// d_s d_t d_eu27 d_JA d_AU /// lnemp y_2 y_3 y_4, vce (robust) *outreg2 using test1.doc, sideway label dec(3) replace **** Column (II): Restriction on CO2 Measurement*********************** ***** Measurement Index > 10 (Average)**************** ******************************************************** keep if m_measure>10 *keep if m01>1 regress ln_CO2gr /// p0*_s dm0*_s d0*_s /// d_s d_t d_eu27 d_JA d_AU /// lnemp y_2 y_3 y_4, vce (robust) *outreg2 using test2.doc, sideway label replace dec(3) **Column (III): drop Banks and Financial Services and Column (II) restriction ************************************************************ preserve drop if icbsectordescription=="8350 Banks" | icbsectordescription=="8770 Financial Services" |icbsectordescription=="8530 Nonlife Insurance"|icbsectordescription=="8570 Life Insurance"|icbsectordescription=="8630 Real Estate Investment & Services"|icbsectordescription=="8670 Real Estate Investment Trusts" regress ln_CO2gr /// p0*_s dm0*_s d0*_s /// d_s d_t d_eu27 d_JA d_AU /// lnemp y_2 y_3 y_4, vce (robust) *outreg2 using test2.doc, sideway label dec(3) restore preserve *****Column (IV): Only firms in the most dirty sectors and Column (II) restriction ************************************************************** keep if sectordirtyness==1| sectordirtyness==2 regress ln_CO2gr /// p0*_s dm0*_s d0*_s /// d_s d_t d_eu27 d_JA d_AU /// lnemp y_2 y_3 y_4, vce (robust) *outreg2 using test2.doc, sideway label dec(3) restore preserve ********Column (V): Only Firms in the ETS sectors and Column (II) restriction keep if eu_ets==1 regress ln_CO2gr /// p0*_s dm0*_s d0*_s /// d_s d_t /// lnemp y_2 y_3 y_4, vce (robust) *outreg2 using test2.doc, sideway label dec(3) restore log close