/****************************************************************************
* *
* 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 */