cd "${output_path}\tmp\"

eststo clear
use "Linked_data.dta", clear

keep if age>=25 & age<=40
gen age2=age*age
gen age3=age2*age

gen vage2=vage*vage
gen vage3=vage2*vage
gen vage4=vage3*vage
gen vage5=vage4*vage

*We create variables with the cut-off point for each observation to calculate new running variables centered around the cut-off
gen vdisc60=60+1/12
replace vdisc60=60+2/12 if vgebjahr==1952 & vgebmonat==1
replace vdisc60=60+3/12 if vgebjahr==1952 & vgebmonat==2
replace vdisc60=60+4/12 if vgebjahr==1952 & vgebmonat==3
replace vdisc60=60+5/12 if vgebjahr==1952 & vgebmonat==4
replace vdisc60=60+6/12 if vgebjahr==1952 & vgebmonat==5
replace vdisc60=60+7/12 if vgebjahr==1952 & vgebmonat>=6 & vgebmonat<=12
replace vdisc60=60+8/12 if vgebjahr==1953 
replace vdisc60=60+9/12 if vgebjahr==1954
replace vdisc60=60+10/12 if vgebjahr==1955
replace vdisc60=60+11/12 if vgebjahr==1956
replace vdisc60=61 if vgebjahr==1957
 
gen vdisc65=65+1/12
replace vdisc65=65+2/12 if vgebjahr==1947
replace vdisc65=65+3/12 if vgebjahr==1948
replace vdisc65=65+4/12 if vgebjahr==1949
replace vdisc65=65+5/12 if vgebjahr==1950
replace vdisc65=65+6/12 if vgebjahr==1951
replace vdisc65=65+7/12 if vgebjahr==1952
replace vdisc65=65+8/12 if vgebjahr==1953
replace vdisc65=65+9/12 if vgebjahr==1954
replace vdisc65=65+10/12 if vgebjahr==1955
replace vdisc65=65+11/12 if vgebjahr==1956
replace vdisc65=66 if vgebjahr==1957

sort pid syear
xtset pid syear
gen fbirth=f.birth





gen mage2=mage*mage
gen mage3=mage2*mage
gen mage4=mage3*mage


gen mdisc60=60+1/12
replace mdisc60=60+2/12 if mgebjahr==1952 & mgebmonat==1
replace mdisc60=60+3/12 if mgebjahr==1952 & mgebmonat==2
replace mdisc60=60+4/12 if mgebjahr==1952 & mgebmonat==3
replace mdisc60=60+5/12 if mgebjahr==1952 & mgebmonat==4
replace mdisc60=60+6/12 if mgebjahr==1952 & mgebmonat==5
replace mdisc60=60+7/12 if mgebjahr==1952 & mgebmonat>=6 & mgebmonat<=12
replace mdisc60=60+8/12 if mgebjahr==1953 
replace mdisc60=60+9/12 if mgebjahr==1954
replace mdisc60=60+10/12 if mgebjahr==1955
replace mdisc60=60+11/12 if mgebjahr==1956
replace mdisc60=61 if mgebjahr==1957
 
gen mdisc65=65+1/12
replace mdisc65=65+2/12 if mgebjahr==1947
replace mdisc65=65+3/12 if mgebjahr==1948
replace mdisc65=65+4/12 if mgebjahr==1949
replace mdisc65=65+5/12 if mgebjahr==1950
replace mdisc65=65+6/12 if mgebjahr==1951
replace mdisc65=65+7/12 if mgebjahr==1952
replace mdisc65=65+8/12 if mgebjahr==1953
replace mdisc65=65+9/12 if mgebjahr==1954
replace mdisc65=65+10/12 if mgebjahr==1955
replace mdisc65=65+11/12 if mgebjahr==1956
replace mdisc65=66 if mgebjahr==1957


gen mage_disc60=mdisc60*mage
gen mage2_disc60=mdisc60*mage2
gen mage3_disc60=mdisc60*mage3

gen mage_disc65=mdisc65*mage
gen mage2_disc65=mdisc65*mage2
gen mage3_disc65=mdisc65*mage3


gen vage_run=vage-vdisc60
gen mage_run=mage-mdisc60

cap log close
set logtype text
log using "${output_path}\Results\TableA3.txt", replace
preserve
keep if vage>=50 & vage<=70

rddensity vage_run, c(0) plot graph_options(title("Density continuity test") xtitle("Father's age (centered)"))
graph export "${output_path}\Results\FigA9_A1.tif", replace width(1200)
rddensity vage, c(60) plot graph_options(title("Density continuity test") xtitle("Father's age"))
graph export "${output_path}\Results\FigA9_A2.tif", replace width(1200)


restore
keep if mage>=50 & mage<=70


rddensity mage_run, c(0) plot graph_options(title("Density continuity test") xtitle("Mother's age (centered)"))
graph export "${output_path}\Results\FigA9_B1.tif", replace width(1200)

rddensity mage, c(60) plot graph_options(title("Density continuity test") xtitle("Mother's age"))
graph export "${output_path}\Results\FigA9_B2.tif", replace width(1200)


log close