/**************************************************************************** * * * Code for Arouna, Mihcler, and Lokossou (2020) * * * ****************************************************************************/ cap log close clear all cd "C:\Users\jdmichler\Dropbox\Contract Farming\JDE CF Data\" set pformat %5.4f log using "CF_results.smcl", replace *Load data from Aminou use Benin_panel_20200923, clear *Set up regression globals quietly tabulate arrond, generate(_a) global depvar rsize ryield psold inc_pc global depvar15 rsize_15 ryield_15 psold_15 inc_pc_15 global arond _a2- _a33 global hhvar15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 agric_15 rice_forma_15 assoc_15 global trt T1C T2C T3C T2T1 T3T1 T3T2 /**************************************************************************** * * * Table 1: Attributes of existing contracts * * * ****************************************************************************/ tab cont_know tab cont_adopt tab cont_write1 tab cont_price1 tab cont_qual1 tab cont_quant1 tab cont_train1 tab cont_cred1 /**************************************************************************** * * * Figure 3: Outcomes by treatment group * * * ****************************************************************************/ by typ_contract, sort: sum rsize ryield psold inc_pc if year == 2016 twoway (kdensity rsize if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) /// (kdensity rsize if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) /// (kdensity rsize if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) /// (kdensity rsize if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) /// ytitle("Density") xtitle("Rice area (ha)") /// xline(.7721364, lpattern(solid) lwidth(thin) lcolor(black)) /// text(.18 `=.7721364-.01' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(1.007895, lpattern(solid) lwidth(thin) lcolor(gs10)) /// text(.18 `=1.007895-.01' "T1 mean", color(gs10) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(.9035317, lpattern(solid) lwidth(thin) lcolor(sky)) /// text(.18 `=.9035317-.01' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(1.030595, lpattern(solid) lwidth(thin) lcolor(turquoise)) /// text(.18 `=1.030595+.01' "T3 mean", color(turquoise) j(left) size(vsmall) place(ne) orient(vertical))), /// legend(label(1 "Control [C]") label(2 "Price [T1]") /// label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3]") /// position(6) col(2) row(2)) saving(rsize, replace) graph export "rsize.emf", as(emf) replace twoway (kdensity ryield if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) /// (kdensity ryield if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) /// (kdensity ryield if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) /// (kdensity ryield if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) /// ytitle("Density") xtitle("Yield (kg/ha)") /// xline(1652.126, lpattern(solid) lwidth(thin) lcolor(black)) /// text(.00009 `=1652.126-10' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(2132.92, lpattern(solid) lwidth(thin) lcolor(gs10)) /// text(.00009 `=2132.92-10' "T1 mean", color(gs10) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(2036.185, lpattern(solid) lwidth(thin) lcolor(sky)) /// text(.00009 `=2036.185-10' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(2220.142, lpattern(solid) lwidth(thin) lcolor(turquoise)) /// text(.00009 `=2220.142-10' "T3 mean", color(turquoise) j(left) size(vsmall) place(nw) orient(vertical))), /// legend(label(1 "Control [C]") label(2 "Price [T1]") /// label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3] ") /// position(6) col(2) row(2)) saving(ryield, replace) graph export "ryield.emf", as(emf) replace twoway (kdensity psold if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) /// (kdensity psold if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) /// (kdensity psold if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) /// (kdensity psold if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) /// ytitle("Density") xtitle("Market participation (%)") /// xline(24.95936, lpattern(solid) lwidth(thin) lcolor(black)) /// text(.0045 `=24.95936-.5' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(51.13158, lpattern(solid) lwidth(thin) lcolor(gs10)) /// text(.0045 `=51.13158-.5' "T1 mean", color(gs10) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(57.25, lpattern(solid) lwidth(thin) lcolor(sky)) /// text(.0045 `=57.25-.5' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(67.15613, lpattern(solid) lwidth(thin) lcolor(turquoise)) /// text(.0045 `=67.15613-.5' "T3 mean", color(turquoise) j(left) size(vsmall) place(nw) orient(vertical))), /// legend(label(1 "Control [C]") label(2 "Price [T1]") /// label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3] ") /// position(6) col(2) row(2)) saving(psold, replace) graph export "psold.emf", as(emf) replace twoway (kdensity inc_pc if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) /// (kdensity inc_pc if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) /// (kdensity inc_pc if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) /// (kdensity inc_pc if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) /// ytitle("Density") xtitle("Income per capita (US$)") /// xline(265.2537, lpattern(solid) lwidth(thin) lcolor(black)) /// text(-.0001 `=265.2537-5' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(378.8236, lpattern(solid) lwidth(thin) lcolor(gs10)) /// text(-.0001 `=378.8236+5' "T1 mean", color(gs10) j(left) size(vsmall) place(ne) orient(vertical)) /// xline(370.9411, lpattern(solid) lwidth(thin) lcolor(sky)) /// text(-.0001 `=370.9411-5' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) /// xline(458.6367, lpattern(solid) lwidth(thin) lcolor(turquoise)) /// text(-.0001 `=458.6367+5' "T3 mean", color(turquoise) j(left) size(vsmall) place(ne) orient(vertical))), /// legend(label(1 "Control [C]") label(2 "Price [T1]") /// label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3] ") /// position(6) col(2) row(2)) saving(inc_pc, replace) graph export "inc_pc.emf", as(emf) replace gr combine rsize.gph ryield.gph psold.gph inc_pc.gph, col(2) iscale(.5) commonscheme graph export "outcomes.emf", as(emf) replace /**************************************************************************** * * * Table 2, B.1: Baseline characteristics and balance tests * * * ****************************************************************************/ by treated, sort: eststo: quietly estpost sum $depvar15 $hhvar15 esttab using sum_stat_treated.csv, replace cell(mean(fmt(3)) sd(par fmt(3))) p(3) label nodepvar by treatment, sort: eststo: quietly estpost sum $depvar15 $hhvar15 esttab using sum_stat_treatment.csv, replace cell(mean(fmt(3)) sd(par fmt(3))) p(3) label nodepvar foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' treated, vce(cl groupement) } // Compare T with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' T1C, vce(cl groupement) } // Compare T1 with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' T2C, vce(cl groupement) } // Compare T2 with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' T3C, vce(cl groupement) } // Compare T3 with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' T2T1, vce(cl groupement) } // Compare T2 with T1 foreach dvar of varlist $depvar15 $hhvar15{ eststo: quietly reg `dvar' T3T1, vce(cl groupement) } // Compare T3 with T1 foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' T3T2, vce(cl groupement) } // Compare T3 with T2 esttab, se nostar drop(_cons) matrix C = r(coefs) eststo clear local rnames : rownames C local models : coleq C local models : list uniq models local i 0 foreach name of local rnames { local ++i local j 0 capture matrix drop b capture matrix drop se foreach model of local models { local ++j matrix tmp = C[`i', 2*`j'-1] if tmp[1,1]<. { matrix colnames tmp = `model' matrix b = nullmat(b), tmp matrix tmp[1,1] = C[`i', 2*`j'] matrix se = nullmat(se), tmp } } ereturn post b quietly estadd matrix se eststo `name' } esttab using balance_new.csv, cells(b(star fmt(3)) se(par fmt(3))) /// noobs label replace /**************************************************************************** * * * Table 3: Baseline characteristics and balance tests * * * ****************************************************************************/ foreach dvar of varlist $trt { eststo: reg `dvar' $depvar15 $hhvar15 $arond, vce(cl groupement) test $depvar $hhvar } // Compare across all treatments esttab, se nostar drop(_cons) esttab using balance_ftest.csv, cells(b(star fmt(3)) se(par fmt(3))) /// noobs label replace /**************************************************************************** * * * Table 4: Tests for attrition bias * * * ****************************************************************************/ preserve gen attrited = 1 if typ_contract == 4 replace attrited = 0 if attrited == . by attrited, sort: eststo: quietly estpost sum treated $depvar15 $hhvar15 esttab using sum_stat_attrit.csv, replace cell(mean(fmt(3)) sd(par fmt(3))) p(3) label nodepvar foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' attrited if treated == 0 | treated == 1, vce(cl groupement) } // Compare T with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' attrited if treatment == 0 | treatment == 1, vce(cl groupement) } // Compare T1 with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' attrited if treatment == 0 | treatment == 2, vce(cl groupement) } // Compare T2 with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' attrited if treatment == 0 | treatment == 3, vce(cl groupement) } // Compare T3 with C foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' attrited if treatment == 1 | treatment == 2, vce(cl groupement) } // Compare T2 with T1 foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' attrited if treatment == 1 | treatment == 3, vce(cl groupement) } // Compare T3 with T1 foreach dvar of varlist $depvar15 $hhvar15 { eststo: quietly reg `dvar' attrited if treatment == 2 | treatment == 3, vce(cl groupement) } // Compare T3 with T2 esttab, se nostar drop(_cons) matrix C = r(coefs) eststo clear local rnames : rownames C local models : coleq C local models : list uniq models local i 0 foreach name of local rnames { local ++i local j 0 capture matrix drop b capture matrix drop se foreach model of local models { local ++j matrix tmp = C[`i', 2*`j'-1] if tmp[1,1]<. { matrix colnames tmp = `model' matrix b = nullmat(b), tmp matrix tmp[1,1] = C[`i', 2*`j'] matrix se = nullmat(se), tmp } } ereturn post b quietly estadd matrix se eststo `name' } esttab using balance_attrit.csv, cells(b(star fmt(3)) se(par fmt(3))) /// noobs label replace * attrition on treatment reg attrited treatment, vce(cl groupement) reg attrited treatment $arond, vce(cl groupement) reg attrited treatment $arond $hhvar15, vce(cl groupement) restore /**************************************************************************** * * * Table B.2: Coefficient of variation in outcomes * * * ****************************************************************************/ by typ_contract, sort: cv2 $depvar /**************************************************************************** * * * Table 5, B.3, B.4, & B.5: Treatment effect of contract * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Rice area */ /* OLS */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_1, replace): /// reg rsize treated $arond, vce(cl groupement) outreg2 using "CTa_jde.doc", /// title (Table 5a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_2, replace): /// reg rsize treated $hhvar15 $arond, vce(cl groupement) outreg2 using "CTa_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_3, replace): /// reg rsize treated rsize_15 $arond, vce(cl groupement) outreg2 using "CTa_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_4, replace): /// reg rsize treated rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "CTa_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: yield */ /* OLS */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_1, replace): /// reg ryield treated $arond, vce(cl groupement) outreg2 using "CTb_jde.doc", /// title (Table 5b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_2, replace): /// reg ryield treated $hhvar15 $arond, vce(cl groupement) outreg2 using "CTb_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_3, replace): /// reg ryield treated ryield_15 $arond, vce(cl groupement) outreg2 using "CTb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_4, replace): /// reg ryield treated ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "CTb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Market participation */ /* OLS */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_1, replace): /// reg psold treated $arond, vce(cl groupement) outreg2 using "CTc_jde.doc", /// title (Table 5c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_2, replace): /// reg psold treated $hhvar15 $arond, vce(cl groupement) outreg2 using "CTc_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_3, replace): /// reg psold treated psold_15 $arond, vce(cl groupement) outreg2 using "CTc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_4, replace): /// reg psold treated psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "CTc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel D: Income */ /* OLS */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_1, replace): /// reg inc_pc treated inc_pc_15 $arond, vce(cl groupement) outreg2 using "CTd_jde.doc", /// title (Table 5d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_2, replace): /// reg inc_pc treated $hhvar15 $arond, vce(cl groupement) outreg2 using "CTd_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_3, replace): /// reg inc_pc treated inc_pc_15 $arond, vce(cl groupement) outreg2 using "CTd_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_4, replace): /// reg inc_pc treated inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "CTd_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) *Mean dependent variable in control group sum rsize ryield psold inc_pc if treated == 0 & year == 2016 /* MHT */ mhtexp rsize ryield psold inc_pc, treatment(treated) restore /* Lee Bounds */ eststo: quietly leebounds rsize treated eststo: quietly leebounds ryield treated eststo: quietly leebounds psold treated eststo: quietly leebounds inc_pc treated esttab using Lee1.csv, cells(b(star fmt(3)) se(par fmt(3))) /// noobs label replace /**************************************************************************** * * * Table 6, 7, B.6, & B.7: Treatment effects of each contract characteristic * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Rice area */ /* OLS */ reg rsize i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPa_jde.doc", /// title (Table 6a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg rsize i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPa_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* ANCOVA */ reg rsize i.typ_contract rsize_15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPa_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg rsize i.typ_contract rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPa_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* Panel B: yield */ /* OLS */ reg ryield i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPb_jde.doc", /// title (Table 6b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg ryield i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPb_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* ANCOVA */ reg ryield i.typ_contract ryield_15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg ryield i.typ_contract ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* Panel C: Market participation */ /* OLS */ reg psold i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPc_jde.doc", /// title (Table 6c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg psold i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPc_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* ANCOVA */ reg psold i.typ_contract psold_15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg psold i.typ_contract psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* Panel D: Income */ /* OLS */ reg inc_pc i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPd_jde.doc", /// title (Table 6d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg inc_pc i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPd_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* ANCOVA */ reg inc_pc i.typ_contract inc_pc_15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPd_jde.doc", append ctitle (ANCOVA) /// drop( _a* inc_pc_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg inc_pc i.typ_contract inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPd_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) /* MHT */ mhtexp rsize ryield psold inc_pc, treatment(typ_contract) /* Randomization inference */ /* Panel A: Rice area */ /* OLS */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_11, replace): /// reg rsize i.treatment $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_rsize_12, replace): /// reg rsize i.treatment $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_13, replace): /// reg rsize i.treatment $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_21, replace): /// reg rsize i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_rsize_22, replace): /// reg rsize i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_23, replace): /// reg rsize i.treatment $hhvar15 $arond, vce(cl groupement) /* ANCOVA */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_31, replace): /// reg rsize i.treatment rsize_15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_rsize_32, replace): /// reg rsize i.treatment rsize_15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_33, replace): /// reg rsize i.treatment rsize_15 $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_41, replace): /// reg rsize i.treatment rsize_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_rsize_42, replace): /// reg rsize i.treatment rsize_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_43, replace): /// reg rsize i.treatment rsize_15 $hhvar15 $arond, vce(cl groupement) /* Panel B: yield */ /* OLS */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_11, replace): /// reg ryield i.treatment $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_12, replace): /// reg ryield i.treatment $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_13, replace): /// reg ryield i.treatment $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_21, replace): /// reg ryield i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_22, replace): /// reg ryield i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_23, replace): /// reg ryield i.treatment $hhvar15 $arond, vce(cl groupement) /* ANCOVA */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_31, replace): /// reg ryield i.treatment ryield_15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_32, replace): /// reg ryield i.treatment ryield_15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_33, replace): /// reg ryield i.treatment ryield_15 $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_41, replace): /// reg ryield i.treatment ryield_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_42, replace): /// reg ryield i.treatment ryield_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_43, replace): /// reg ryield i.treatment ryield_15 $hhvar15 $arond, vce(cl groupement) /* Panel C: Market participation */ /* OLS */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_11, replace): /// reg psold i.treatment $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_12, replace): /// reg psold i.treatment $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_13, replace): /// reg psold i.treatment $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_21, replace): /// reg psold i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_22, replace): /// reg psold i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_23, replace): /// reg psold i.treatment $hhvar15 $arond, vce(cl groupement) /* ANCOVA */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_31, replace): /// reg psold i.treatment psold_15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_32, replace): /// reg psold i.treatment psold_15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_33, replace): /// reg psold i.treatment psold_15 $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_41, replace): /// reg psold i.treatment psold_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_42, replace): /// reg psold i.treatment psold_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_43, replace): /// reg psold i.treatment psold_15 $hhvar15 $arond, vce(cl groupement) /* Panel D: Income */ /* OLS */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_11, replace): /// reg inc_pc i.treatment $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_12, replace): /// reg inc_pc i.treatment $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_13, replace): /// reg inc_pc i.treatment $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_21, replace): /// reg inc_pc i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_22, replace): /// reg inc_pc i.treatment $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_23, replace): /// reg inc_pc i.treatment $hhvar15 $arond, vce(cl groupement) /* ANCOVA */ ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_31, replace): /// reg inc_pc i.treatment inc_pc_15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_32, replace): /// reg inc_pc i.treatment inc_pc_15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_33, replace): /// reg inc_pc i.treatment inc_pc_15 $arond, vce(cl groupement) ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_41, replace): /// reg inc_pc i.treatment inc_pc_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_42, replace): /// reg inc_pc i.treatment inc_pc_15 $hhvar15 $arond, vce(cl groupement) ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) /// reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_43, replace): /// reg inc_pc i.treatment inc_pc_15 $hhvar15 $arond, vce(cl groupement) restore /**************************************************************************** * * * Table 8: Heterogeneity of treatment effects of contract * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Contract Indicator */ /* Column 1: Rice area */ reg rsize i.treated##c.hhsize_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", /// title (Table 8: Treatment effects on Rice area) ctitle (hhsize) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg rsize i.treated##c.age_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", ctitle (age) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg rsize i.treated##i.sex_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", ctitle (sex) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg rsize i.treated##c.rice_exp_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", ctitle (exp) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg rsize i.treated##i.primeduc_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", ctitle (edu) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg rsize i.treated##i.agric_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", ctitle (activity) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg rsize i.treated##i.rice_forma_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", ctitle (train) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg rsize i.treated##i.assoc_15 rsize_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_area_jde.doc", ctitle (assoc) /// drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Column 2: yield */ reg ryield i.treated##c.hhsize_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", /// title (Table 8: Treatment effects on Rice yield) ctitle (hhsize) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg ryield i.treated##c.age_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", ctitle (age) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg ryield i.treated##i.sex_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", ctitle (sex) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg ryield i.treated##c.rice_exp_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", ctitle (exp) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg ryield i.treated##i.primeduc_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", ctitle (edu) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg ryield i.treated##i.agric_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", ctitle (activity) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg ryield i.treated##i.rice_forma_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", ctitle (train) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg ryield i.treated##i.assoc_15 ryield_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_yield_jde.doc", ctitle (assoc) /// drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Column 3: Market participation */ reg psold i.treated##c.hhsize_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", /// title (Table 8: Treatment effects on Market participation) ctitle (hhsize) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg psold i.treated##c.age_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", ctitle (age) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg psold i.treated##i.sex_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", ctitle (sex) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg psold i.treated##c.rice_exp_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", ctitle (exp) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg psold i.treated##i.primeduc_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", ctitle (edu) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg psold i.treated##i.agric_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", ctitle (activity) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg psold i.treated##i.rice_forma_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", ctitle (train) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg psold i.treated##i.assoc_15 psold_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_market_jde.doc", ctitle (assoc) /// drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Column 4: Income */ reg inc_pc i.treated##c.hhsize_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", /// title (Table 8: Treatment effects on Income) ctitle (hhsize) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg inc_pc i.treated##c.age_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", ctitle (age) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg inc_pc i.treated##i.sex_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", ctitle (sex) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg inc_pc i.treated##c.rice_exp_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", ctitle (exp) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg inc_pc i.treated##i.primeduc_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", ctitle (edu) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg inc_pc i.treated##i.agric_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", ctitle (activity) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg inc_pc i.treated##i.rice_forma_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", ctitle (train) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg inc_pc i.treated##i.assoc_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) outreg2 using "h_inc_jde.doc", ctitle (assoc) /// drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 /// agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /**************************************************************************** * * * Figures 4-6: Heterogeneity of treatment effects of contract * * * ****************************************************************************/ /* Contract Types */ /* Column 1: Rice area */ qui reg rsize i.typ_contract##c.hhsize_15 rsize_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(hhsize_15=(0(2)20)) marginsplot, title("") ytitle("Rice area (ha)") xtitle("Household size (n)") /// legend(pos(6) col(2) row(2)) saving(rsize_1, replace) graph export "rsize_1.emf", as(emf) replace qui reg rsize i.typ_contract##c.rice_exp_15 rsize_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(rice_exp_15=(0(5)30)) marginsplot, title("") ytitle("Rice area (ha)") xtitle("Exp. producing rice (years)") /// legend(pos(6) col(2) row(2)) saving(rsize_4, replace) graph export "rsize_4.emf", as(emf) replace qui reg rsize i.typ_contract##i.rice_forma_15 rsize_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, over(rice_forma_15) marginsplot, title("") ytitle("Rice area (ha)") xtitle("Farming contracts") /// xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(rsize_7, replace) graph export "rsize_7.emf", as(emf) replace /* Column 2: yield */ qui reg ryield i.typ_contract##c.hhsize_15 ryield_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(hhsize_15=(0(2)20)) marginsplot, title("") ytitle("Yield (kg/ha)") xtitle("Household size (n)") /// legend(pos(6) col(2) row(2)) saving(ryield_1, replace) graph export "ryield_1.emf", as(emf) replace qui reg ryield i.typ_contract##c.rice_exp_15 ryield_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(rice_exp_15=(0(5)30)) marginsplot, title("") ytitle("Yield (kg/ha)") xtitle("Exp. producing rice (years)") /// legend(pos(6) col(2) row(2)) saving(ryield_4, replace) graph export "ryield_4.emf", as(emf) replace qui reg ryield i.typ_contract##i.rice_forma_15 ryield_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, over(rice_forma_15) marginsplot, title("") ytitle("Yield (kg/ha)") xtitle("Farming contracts") /// xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(ryield_7, replace) graph export "ryield_7.emf", as(emf) replace /* Column 3: market participation */ qui reg psold i.typ_contract##c.hhsize_15 psold_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(hhsize_15=(0(2)20)) marginsplot, title("") ytitle("Market participation (%)") xtitle("Household size (n)") /// legend(pos(6) col(2) row(2)) saving(psold_1, replace) graph export "psold_1.emf", as(emf) replace qui reg psold i.typ_contract##c.rice_exp_15 psold_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(rice_exp_15=(0(5)30)) marginsplot, title("") ytitle("Market participation (%)") xtitle("Exp. producing rice (years)") /// legend(pos(6) col(2) row(2)) saving(psold_4, replace) graph export "psold_4.emf", as(emf) replace qui reg psold i.typ_contract##i.rice_forma_15 psold_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, over(rice_forma_15) marginsplot, title("") ytitle("Market participation (%)") xtitle("Farming contracts") /// xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(psold_7, replace) graph export "psold_7.emf", as(emf) replace /* Column 4: income */ qui reg inc_pc i.typ_contract##c.hhsize_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(hhsize_15=(0(2)20)) marginsplot, title("") ytitle("Income per capita (US$)") xtitle("Household size (n)") /// legend(pos(6) col(2)) saving(inc_pc_1, replace) graph export "inc_pc_1.emf", as(emf) replace qui reg inc_pc i.typ_contract##c.rice_exp_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, at(rice_exp_15=(0(5)30)) marginsplot, title("") ytitle("Income per capita (US$)") xtitle("Exp. producing rice (years)") /// legend(pos(6) col(2)) saving(inc_pc_4, replace) graph export "inc_pc_4.emf", as(emf) replace qui reg inc_pc i.typ_contract##i.rice_forma_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement) qui margins typ_contract, over(rice_forma_15) marginsplot, title("") ytitle("Income per capita (US$)") xtitle("Farming contracts") /// xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(inc_pc_7, replace) graph export "inc_pc_7.emf", as(emf) replace *Combine selected graphs gr combine rsize_1.gph ryield_1.gph psold_1.gph inc_pc_1.gph, col(2) iscale(.5) commonscheme graph export "hetero_hhsize.emf", as(emf) replace gr combine rsize_4.gph ryield_4.gph psold_4.gph inc_pc_4.gph, col(2) iscale(.5) commonscheme graph export "hetero_exp.emf", as(emf) replace gr combine rsize_7.gph ryield_7.gph psold_7.gph inc_pc_7.gph, col(2) iscale(.5) commonscheme graph export "hetero_train.emf", as(emf) replace restore /**************************************************************************** * * * Table B.3-B.3: Treatment effect on input use * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Table B.3: Treated */ reg lns treated lns15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTin_jde.doc", /// title (Table B.3: Treatment effects on Inputs) ctitle (ln(seed)) /// drop(_a* lns15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg lnf treated lnf15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTin_jde.doc", append ctitle (ln(fert)) /// drop(_a* lnf15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg lnp treated lnp15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTin_jde.doc", append ctitle (ln(pest)) /// drop(_a* lnp15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg lnh treated lnh15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTin_jde.doc", append ctitle (ln(herb)) /// drop(_a* lnh15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) reg lnl treated lnl15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTin_jde.doc", append ctitle (ln(labr)) /// drop(_a* lnl15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) *Mean dependent variable in control group sum lns lnf lnp lnh lnl if treated == 0 & year == 2016 /* Table B.4: Treatment Types */ reg lns i.typ_contract lns15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPin_jde.doc", /// title (Table B.4: Treatment effects on Inputs) ctitle (ln(seed)) /// drop(_a* lns15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg lnf i.typ_contract lnf15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPin_jde.doc", append ctitle (ln(fert)) /// drop(_a* lnf15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg lnp i.typ_contract lnp15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPin_jde.doc", append ctitle (ln(pest)) /// drop(_a* lnp15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg lnh i.typ_contract lnh15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPin_jde.doc", append ctitle (ln(herb)) /// drop(_a* lnh15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) reg lnl i.typ_contract lnl15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPin_jde.doc", append ctitle (ln(labr)) /// drop(_a* lnl15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) /// (1.typ_contract=3.typ_contract), common cons mtest(b) restore /**************************************************************************** * * * Table B.10: Treatment effect of price and quantity guarantee * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 drop if T1C == . /* Panel A: Rice area */ /* OLS */ reg rsize T1C $arond if year == 2016, vce(cl groupement) outreg2 using "CT1a_jde.doc", /// title (Table B.10a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T1C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1a_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg rsize T1C rsize_15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T1C rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: yield */ /* OLS */ reg ryield T1C $arond if year == 2016, vce(cl groupement) outreg2 using "CT1b_jde.doc", /// title (Table B.10b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T1C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1b_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg ryield T1C ryield_15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T1C ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Market participation */ /* OLS */ reg psold T1C $arond if year == 2016, vce(cl groupement) outreg2 using "CT1c_jde.doc", /// title (Table B.10c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T1C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1c_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg psold T1C psold_15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T1C psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel D: Income */ /* OLS */ reg inc_pc T1C $arond if year == 2016, vce(cl groupement) outreg2 using "CT1d_jde.doc", /// title (Table B.10d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No, /// Household FE, No, Year FE, No) reg inc_pc T1C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1d_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg inc_pc T1C inc_pc_15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T1C inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CT1d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) restore /**************************************************************************** * * * Table B.11: Treatment effect of price and extension training * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 drop if T2C == . /* Panel A: Rice area */ /* OLS */ reg rsize T2C $arond if year == 2016, vce(cl groupement) outreg2 using "T2Ca_jde.doc", /// title (Table B.11a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T2C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Ca_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg rsize T2C rsize_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Ca_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T2C rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Ca_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: yield */ /* OLS */ reg ryield T2C $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cb_jde.doc", /// title (Table B.11b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T2C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cb_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg ryield T2C ryield_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T2C ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Market participation */ /* OLS */ reg psold T2C $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cc_jde.doc", /// title (Table B.11c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T2C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cc_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg psold T2C psold_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T2C psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel D: Income */ /* OLS */ reg inc_pc T2C $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cd_jde.doc", /// title (Table B.11d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T2C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cd_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg inc_pc T2C inc_pc_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cd_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T2C inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2Cd_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) restore /**************************************************************************** * * * Table B.12: Treatment effect of price, extension training, and loans * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 drop if T3C == . /* Panel A: Rice area */ /* OLS */ reg rsize T3C $arond if year == 2016, vce(cl groupement) outreg2 using "T3Ca_jde.doc", /// title (Table B.12a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T3C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Ca_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg rsize T3C rsize_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Ca_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T3C rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Ca_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: yield */ /* OLS */ reg ryield T3C $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cb_jde.doc", /// title (Table B.12b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T3C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cb_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg ryield T3C ryield_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T3C ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cb_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Market participation */ /* OLS */ reg psold T3C $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cc_jde.doc", /// title (Table B.12c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T3C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cc_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg psold T3C psold_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T3C psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cc_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel D: Income */ /* OLS */ reg inc_pc T3C $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cd_jde.doc", /// title (Table B.12d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T3C $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cd_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg inc_pc T3C inc_pc_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cd_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T3C inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T3Cd_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) restore /**************************************************************************** * * * Table B.13: Treatment effects of extension training * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 drop if T2T1 == . /* Panel A: Rice area */ /* OLS */ reg rsize T2T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2a_jde.doc", /// title (Table B.13a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2a_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg rsize T2T1 rsize_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T2T1 rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: yield */ /* OLS */ reg ryield T2T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2b_jde.doc", /// title (Table B.13b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2b_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg ryield T2T1 ryield_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T2T1 ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Market participation */ /* OLS */ reg psold T2T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2c_jde.doc", /// title (Table B.13c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2c_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg psold T2T1 psold_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T2T1 psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel D: Income */ /* OLS */ reg inc_pc T2T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2d_jde.doc", /// title (Table B.13d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2d_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg inc_pc T2T1 inc_pc_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T2T1 inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T2d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) sum rsize ryield psold inc_pc if treatment == 1 & year == 2016 restore /**************************************************************************** * * * Table B.14: Treatment effects of input loans * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 drop if T3T2 == . /* Panel A: Rice area */ /* OLS */ reg rsize T3T2 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3a_jde.doc", /// title (Table B.14a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3a_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg rsize T3T2 rsize_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T3T2 rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Yield */ /* OLS */ reg ryield T3T2 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3b_jde.doc", /// title (Table B.14b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3b_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc ) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg ryield T3T2 ryield_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T3T2 ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Market participation */ /* OLS */ reg psold T3T2 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3c_jde.doc", /// title (Table B.14c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3c_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg psold T3T2 psold_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T3T2 psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel D: Income */ /* OLS */ reg inc_pc T3T2 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3d_jde.doc", /// title (Table B.14d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3d_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg inc_pc T3T2 inc_pc_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T3T2 inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T2T3d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) sum rsize ryield psold inc_pc if treatment == 2 & year == 2016 restore /**************************************************************************** * * * Table B.15: Treatment effects of extension and loans * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 drop if T3T1 == . /* Panel A: Rice area */ /* OLS */ reg rsize T3T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3a_jde.doc", /// title (Table B.15a: Treatment effects on Rice area) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3a_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg rsize T3T1 rsize_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg rsize T3T1 rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3a_jde.doc", append ctitle (ANCOVA) /// drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Yield */ /* OLS */ reg ryield T3T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3b_jde.doc", /// title (Table B.15b: Treatment effects on yield) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3b_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg ryield T3T1 ryield_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg ryield T3T1 ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3b_jde.doc", append ctitle (ANCOVA) /// drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Market participation */ /* OLS */ reg psold T3T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3c_jde.doc", /// title (Table B.15c: Treatment effects on market participation) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3c_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg psold T3T1 psold_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg psold T3T1 psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3c_jde.doc", append ctitle (ANCOVA) /// drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel D: Income */ /* OLS */ reg inc_pc T3T1 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3d_jde.doc", /// title (Table B.15d: Treatment effects on income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3d_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg inc_pc T3T1 inc_pc_15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg inc_pc T3T1 inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "T1T3d_jde.doc", append ctitle (ANCOVA) /// drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) restore /**************************************************************************** * * * Table B.16: Treatment effect on income * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Rice Income */ /* OLS */ reg lnri treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTri_jde.doc", /// title (Table B.16a: Treatment effects on Rice Income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnri treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTri_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnri treated lnri15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTri_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnri15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnri treated lnri15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTri_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnri15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Other Farm Income */ /* OLS */ reg lnfi treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTfi_jde.doc", /// title (Table B.16b: Treatment effects on Other Farm Income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnfi treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTfi_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnfi treated lnfi15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTfi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnfi15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnfi treated lnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTfi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Non-Farm Income */ /* OLS */ reg lnnfi treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTnfi_jde.doc", /// title (Table B.16c: Treatment effects on Non-Farm Income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnnfi treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTnfi_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnnfi treated lnnfi15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTnfi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnnfi15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnnfi treated lnnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTnfi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) *Mean dependent variable in control group sum lnri lnfi lnnfi if treated == 0 & year == 2016 restore /**************************************************************************** * * * Table B.17: Treatment effect of each contract on income * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Rice Income */ /* OLS */ reg lnri i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPri_jde.doc", /// title (Table B.17a: Treatment effects on Rice Income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnri i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPri_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnri i.typ_contract lnri15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPri_jde.doc", append ctitle (ANCOVA) /// drop( _a* lnri15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnri i.typ_contract lnri15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPri_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnri15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Other Farm Income */ /* OLS */ reg lnfi i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPfi_jde.doc", /// title (Table B.17b: Treatment effects on Other Farm Income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnfi i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPfi_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnfi i.typ_contract lnfi15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPfi_jde.doc", append ctitle (ANCOVA) /// drop( _a* lnfi15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnfi i.typ_contract lnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPfi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Non-Farm Income */ /* OLS */ reg lnnfi i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPnfi_jde.doc", /// title (Table B.17c: Treatment effects on Non-Farm Income) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnnfi i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPnfi_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnnfi i.typ_contract lnnfi15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPnfi_jde.doc", append ctitle (ANCOVA) /// drop( _a* lnnfi15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnnfi i.typ_contract lnnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPnfi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) restore /**************************************************************************** * * * Table B.18: Treatment effect on food security * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Food Insecurity */ /* OLS */ reg hfias treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTe_jde.doc", /// title (Table B.18a: Treatment effects on HFIAS) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg hfias treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTe_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg hfias treated hfias_15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTe_jde.doc", append ctitle (ANCOVA) /// drop(_a* hfias_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg hfias treated hfias_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTe_jde.doc", append ctitle (ANCOVA) /// drop(_a* hfias_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Food Security */ /* OLS */ reg fcs treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTf_jde.doc", /// title (Table B.18b: Treatment effects on FCS) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg fcs treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTf_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg fcs treated fcs_15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTf_jde.doc", append ctitle (ANCOVA) /// drop(_a* fcs_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg fcs treated fcs_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTf_jde.doc", append ctitle (ANCOVA) /// drop(_a* fcs_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) *Mean dependent variable in control group sum hfias fcs if treated == 0 & year == 2016 restore /**************************************************************************** * * * Table B.19: Treatment effect of each contract on food security * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Food Insecurity */ /* OLS */ reg hfias i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPe_jde.doc", /// title (Table B.19a: Treatment effects on HFIAS) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg hfias i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPe_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg hfias i.typ_contract hfias_15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPe_jde.doc", append ctitle (ANCOVA) /// drop( _a* hfias_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg hfias i.typ_contract hfias_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPe_jde.doc", append ctitle (ANCOVA) /// drop(_a* hfias_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Food Security */ /* OLS */ reg fcs i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPf_jde.doc", /// title (Table B.19b: Treatment effects on HFIAS) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg fcs i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPf_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg fcs i.typ_contract fcs_15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPf_jde.doc", append ctitle (ANCOVA) /// drop( _a* fcs_15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg fcs i.typ_contract fcs_15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPf_jde.doc", append ctitle (ANCOVA) /// drop(_a* fcs_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) restore /**************************************************************************** * * * Table B.20: Treatment effect on profit * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Profits - low wages */ /* OLS */ reg lnpil treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTpil_jde.doc", /// title (Table B.20a: Treatment effects on Profits - low wages) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpil treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpil_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnpil treated lnpil15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpil_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpil15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpil treated lnpil15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpil_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpil15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Profits - high wages */ /* OLS */ reg lnpih treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTpih_jde.doc", /// title (Table B.20b: Treatment effects on Profits - high wages) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpih treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpih_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnpih treated lnpih15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpih_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpih15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpih treated lnpih15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpih_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpih15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Profits - self-reported wages */ /* OLS */ reg lnpi treated $arond if year == 2016, vce(cl groupement) outreg2 using "CTpi_jde.doc", /// title (Table B.20c: Treatment effects on Profits - self-reported wages) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpi treated $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpi_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnpi treated lnpi15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpi15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpi treated lnpi15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "CTpi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) *Mean dependent variable in control group sum lnpil lnpih lnpi if treated == 0 & year == 2016 restore /**************************************************************************** * * * Table B.21: Treatment effect of each contract on profit * * * ****************************************************************************/ *Drop attritors preserve drop if typ_contract == 4 /* Panel A: Profits - low wages */ /* OLS */ reg lnpil i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPil_jde.doc", /// title (Table B.21a: Treatment effects on Profits - low wages) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpil i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPil_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnpil i.typ_contract lnpil15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPil_jde.doc", append ctitle (ANCOVA) /// drop( _a* lnpil15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpil i.typ_contract lnpil15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPil_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpil15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel B: Profits - high wages */ /* OLS */ reg lnpih i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPih_jde.doc", /// title (Table B.21b: Treatment effects on Profits - high wages) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpih i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPih_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnpih i.typ_contract lnpih15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPih_jde.doc", append ctitle (ANCOVA) /// drop( _a* lnpih15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpih i.typ_contract lnpih15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPih_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpih15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* Panel C: Profits - self-reported wages */ /* OLS */ reg lnpi i.typ_contract $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPi_jde.doc", /// title (Table B.21c: Treatment effects on Profits - self-reported wages) ctitle (OLS) /// drop(_a*) se dec(3) label replace nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpi i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPi_jde.doc", append ctitle (OLS) /// drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) /* ANCOVA */ reg lnpi i.typ_contract lnpi15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPi_jde.doc", append ctitle (ANCOVA) /// drop( _a* lnpi15) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, No) reg lnpi i.typ_contract lnpi15 $hhvar15 $arond if year == 2016, vce(cl groupement) outreg2 using "iTYPi_jde.doc", append ctitle (ANCOVA) /// drop(_a* lnpi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* /// rice_forma assoc) se dec(3) label nocons /// addtext(Arrondissement FE, Yes, Household Covariates, Yes) restore /* END */