** The impact of extreme rainfall and flooding on sanitation coverage in rural Malawi pre- and post- tropical Cyclone Freddy
* codebook


insheet using "cyclone_freddy_sanitation_Malawi_data.csv", comma clear

clear


*------------------------------------------------------------------------------------------------------------------------------
** JMP sanitation pre- and post-cyclone -- Table 1 & Figure 2
*------------------------------------------------------------------------------------------------------------------------------

* Percentage of households at each level of the JMP sanitation service pre-cyclone
tab jmp_sanitation_precyclone

	* 1 = basic
	* 2 = limited
	* 3 = unimproved

* Percentage of households at each level of the JMP sanitation service post-cyclone
tab jmp_sanitation_postcyclone

	* 1 = basic
	* 2 = limited
	* 3 = unimproved
	* 4 = no sanitation

*------------------------------------------------------------------------------------------------------------------------------
** Sanitation facility collapse -- Table 2
*------------------------------------------------------------------------------------------------------------------------------

* whether sanitation facility collapsed 
tab latrine_collapsed_cyclone

	* 0 = no
	* 1 = yes

* types of sanitation facilities that collapsed/were damaged
tab jmp_sanitation_precyclone if latrine_collapsed_cyclone == 1

	* 1 = basic
	* 2 = limited
	* 3 = unimproved

* reasons why sanitation facility did not collapse
tab latrine_strength

	* 0 = not sure
	* 1 = the materials used were durable
	* 2 = the cyclone was not strong in our community
	* 3 = other
	* .

	
*------------------------------------------------------------------------------------------------------------------------------
** Sanitation facility rebuilding
*------------------------------------------------------------------------------------------------------------------------------

* percentage of collapsed sanitation facilities that were rebuilt
tab latrine_rebuilt


* types of sanitation facilities that were rebuilt
tab jmp_sanitation_postcyclone if latrine_rebuilt == 1

	* 1 = basic
	* 2 = limited
	* 3 = unimproved

* reasons why the type of sanitation facility was selected
tab latrine_designchoice if latrine_rebuilt == 1 & jmp_sanitation_postcyclone == 3

	* 1 = it's durable
	* 2 = it's cheap
	* 3 = these are the types common in the community
	* 4 = had no option
	* 999 = don't know
	* other

*------------------------------------------------------------------------------------------------------------------------------
** Analysis - latrine collapse
*------------------------------------------------------------------------------------------------------------------------------

* combine basic and limited --> improved sanitation
replace jmp_sanitation_precyclone = 1 if jmp_sanitation_precyclone == 1 | jmp_sanitation_precyclone == 2
replace jmp_sanitation_precyclone = 2 if jmp_sanitation_precyclone == 3
tab jmp_sanitation_precyclone

	* 1 = basic or limited
	* 2 = unimproved

* risk ratio
glm latrine_collapsed_cyclone jmp_sanitation_precyclone, fam(bin) link(log) nolog eform


*------------------------------------------------------------------------------------------------------------------------------
** Analysis - latrine rebuilding
*------------------------------------------------------------------------------------------------------------------------------	

* odds ratio
logit latrine_rebuilt under_five if !missing(household_size), or
logit latrine_rebuilt household_size if !missing(household_size), or


*------------------------------------------------------------------------------------------------------------------------------
** Sankey diagram -- Figure 2
*------------------------------------------------------------------------------------------------------------------------------

rename jmp_sanitation_precyclone baseline
rename jmp_sanitation_postcyclone endline
rename id count
rename latrine_collapsed_cyclone cyclone

* relabel baseline
label define baseline_labels 1 "Basic" 2 "Limited" 3 "Unimproved"
label list baseline_labels
label values baseline baseline_labels

* relabel endline
label define endline_labels 1 "Basic" 2 "Limited" 3 "Unimproved" 4 "No sanitation"
label list endline_labels
label values endline endline_labels

* relabel cyclone
label define cyclone_labels 0 "Latrine did not collapse" 1 "Latrine collapsed"
label list cyclone_labels
label values cyclone cyclone_labels

* collapse from baseline to cyclone
collapse (count) count, by(baseline cyclone)

* collapse from cyclone to endline
collapse (count) count, by(cyclone endline)

* sankey diagram
sankey count, from(source) to(destination) by(layer) palette(CET C6) sort1(name, reverse) nolabel
graph export sankey_diagram_v2.svg, as(svg)