************************************************************
*
*	Project: "Copayment and NH use"
*
*	Dofile "D.2.CP_NHuse_CP_Graphs.do"
*
*	Purpose:	Graphs based on schedule of co-payments
*				Use the schedule coded in D.1.CP_NHuse_Copayments_ForGraphs.do
*
*	Author:		Marianne Tenand
*
*	Created:	March 2020
*	Last modif:	October 24 2021
*
************************************************************

global CBS_RA = 0 
	/* Global indicating whether this script is run within the CBS remote access environment (=1), or outside (=0) */

if $CBS_RA == 1 {
	* Globals for directories
	do	""

	global file ""
	global out  ""
	}

if $CBS_RA == 0 {
    global folder		""
	global file 		""
	global out  		""
	global base_perso  	""
	}


global N=100000
global income_max = 50000
*global wealth_max = 1000000
	
*****************************************
*** A. 2012 vs 2013 on the same graph ***
*****************************************

*** Dataset (fictitious) for the poor ***
*****************************************

* No box_3 wealth
clear
set obs $N
set seed 12345

gen income				=0+_n*($income_max/$N)
gen wea_ass				=0

gen year = 2012

gen type = 1
save "${base_perso}\data_temp_1.dta", replace


* median box_3 wealth
clear
set obs $N
set seed 12345

gen income				=0+_n*($income_max/$N)
gen wea_ass				=27932    
	
gen year = 2012

gen type = 2
save "${base_perso}\data_temp_2.dta", replace


* p75 box_3 wealth
clear
set obs $N
set seed 12345

gen income				=0+_n*($income_max/$N)
gen wea_ass				=92089 

gen year = 2012

gen type = 3
save "${base_perso}\data_temp_3.dta", replace


* p95 box_3 wealth
clear
set obs $N
set seed 12345

gen income				=0+_n*($income_max/$N)
gen wea_ass				=455159

gen year = 2012

gen type = 4

forvalues i=1/3{
	append using "${base_perso}\data_temp_`i'.dta"
	}

sort type

gen inc_box123		 = income+0.04*wea_ass
gen inc_box12		 = income
gen income_disposable = inc_box123

save "${base_perso}\data_temp_4.dta", replace

replace year=2013

append using "${base_perso}\data_temp_4.dta"

	
* Simulate CP *

qui do "${file}"


*** Graphs: high_rate and low-rate co-payments on NH care ***
*************************************************************

foreach var in low high {
#delimit ;
	graph twoway (scatter CP_NH_`var'_pre income if type==1 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_NH_`var'_post income if type==1 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("`var'-rate co-payment""(in euros/month)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(500)2500)
				legend(label(1 "No financial wealth""Pre-reform (2012)") 
				label(2 "No financial""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
					;
	graph export "${out}CP_`var'_schedule_group1.png", as(png) replace
		;
#delimit cr   

#delimit ;
	graph twoway (scatter CP_NH_`var'_pre income if type==2 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_NH_`var'_post income if type==2 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("`var'-rate co-payment""(in euros/month)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(500)2500)
				legend(label(1 "Median financial wealth""Pre-reform (2012)") 
				label(2 "Median financial wealth""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
					;
	graph export "${out}CP_`var'_schedule_group2.png", as(png) replace
		;
#delimit cr   		

#delimit ;
	graph twoway (scatter CP_NH_`var'_pre income if type==3 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_NH_`var'_post income if type==3 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("`var'-rate co-payment""(in euros/month)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(500)2500)
				legend(label(1 "p75 financial wealth""Pre-reform (2012)") 
				label(2 "p75 financial wealth""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
				;
	graph export "${out}CP_`var'_schedule_group3.png", as(png) replace
		;
#delimit cr   

#delimit ;
	graph twoway (scatter CP_NH_`var'_pre income if type==4 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_NH_`var'_post income if type==4 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("`var'-rate co-payment""(in euros/month)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(500)2500)
				legend(label(1 "p95 financial wealth""Pre-reform (2012)") 
				label(2 "p95 financial wealth""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
					;
	graph export "${out}CP_`var'_schedule_group4.png", as(png) replace
		;
#delimit cr   	
}



*** Graphs: maximum co-payment on home care ***
***********************************************

#delimit ;
	graph twoway (scatter CP_HC_max_pre income if type==1 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_HC_max_post income if type==1 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("Maximum co-payment""(in euros/4-week period)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(250)1250)
				legend(label(1 "No financial wealth""Pre-reform (2012)") 
				label(2 "No financial""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
					;
	graph export "${out}CP_max_schedule_group1.png", as(png) replace
		;
#delimit cr   

#delimit ;
	graph twoway (scatter CP_HC_max_pre income if type==2 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_HC_max_post income if type==2 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("Maximum co-payment""(in euros/4-week period)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(250)1250)
				legend(label(1 "Median financial wealth""Pre-reform (2012)") 
				label(2 "Median financial wealth""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
					;
	graph export "${out}CP_max_schedule_group2.png", as(png) replace
		;
#delimit cr   		

#delimit ;
	graph twoway (scatter CP_HC_max_pre income if type==3 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_HC_max_post income if type==3 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("Maximum co-payment""(in euros/4-week period)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(250)1250)
				legend(label(1 "p75 financial wealth""Pre-reform (2012)") 
				label(2 "p75 financial wealth""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
				;
	graph export "${out}CP_max_schedule_group3.png", as(png) replace
		;
#delimit cr   

#delimit ;
	graph twoway (scatter CP_HC_max_pre income if type==4 & year == 2012,
					msize(tiny) mcolor(gs8))
				(scatter CP_HC_max_post income if type==4 & year == 2013,
					msize(tiny) mcolor(gs4)),
				xtitle("Income (in euros/year)") 
				ytitle("Maximum co-payment""(in euros/4-week period)")
				xline(14739 17220 21515, lpattern(dash) lcolor(black))
				ylabel(,angle(0)) ysca(titlegap(2)) 
				xsca(titlegap(2)) ylabel(0(250)1250)
				legend(label(1 "p95 financial wealth""Pre-reform (2012)") 
				label(2 "p95 financial wealth""Post-reform (2013)"))
				plotregion(margin(zero) lcolor(none)) bgcolor(white) graphregion(color(white))
					;
	graph export "${out}CP_max_schedule_group4.png", as(png) replace
		;
#delimit cr   	
  
*** Cleaning ***
                     
forvalues i=1/4{
	erase "${base_perso}\data_temp_`i'.dta"
	}