cd "${output_path}\tmp\"


eststo clear
use "Linked_data.dta", clear

*Generate variables for age and discontinuities
gen age2=age*age

gen vage2=vage*vage


gen vdisc60=cond((vage+0.001)>=(60+1/12),1,0)
replace vdisc60=cond((vage+0.001)>=(60+2/12),1,0) if vgebjahr==1952 & vgebmonat==1
replace vdisc60=cond((vage+0.001)>=(60+3/12),1,0) if vgebjahr==1952 & vgebmonat==2
replace vdisc60=cond((vage+0.001)>=(60+4/12),1,0) if vgebjahr==1952 & vgebmonat==3
replace vdisc60=cond((vage+0.001)>=(60+5/12),1,0) if vgebjahr==1952 & vgebmonat==4
replace vdisc60=cond((vage+0.001)>=(60+6/12),1,0) if vgebjahr==1952 & vgebmonat==5
replace vdisc60=cond((vage+0.001)>=(60+7/12),1,0) if vgebjahr==1952 & vgebmonat>=6 & vgebmonat<=12
replace vdisc60=cond((vage+0.001)>=(60+8/12),1,0) if vgebjahr==1953 
replace vdisc60=cond((vage+0.001)>=(60+9/12),1,0) if vgebjahr==1954
replace vdisc60=cond((vage+0.001)>=(60+10/12),1,0) if vgebjahr==1955
replace vdisc60=cond((vage+0.001)>=(60+11/12),1,0) if vgebjahr==1956
replace vdisc60=cond((vage+0.001)>=61,1,0) if vgebjahr==1957
 
gen vdisc65=cond((vage+0.001)>=(65+1/12),1,0)
replace vdisc65=cond((vage+0.001)>=(65+2/12),1,0) if vgebjahr==1947
replace vdisc65=cond((vage+0.001)>=(65+3/12),1,0) if vgebjahr==1948
replace vdisc65=cond((vage+0.001)>=(65+4/12),1,0) if vgebjahr==1949
replace vdisc65=cond((vage+0.001)>=(65+5/12),1,0) if vgebjahr==1950
replace vdisc65=cond((vage+0.001)>=(65+6/12),1,0) if vgebjahr==1951
replace vdisc65=cond((vage+0.001)>=(65+7/12),1,0) if vgebjahr==1952
replace vdisc65=cond((vage+0.001)>=(65+8/12),1,0) if vgebjahr==1953
replace vdisc65=cond((vage+0.001)>=(65+9/12),1,0) if vgebjahr==1954
replace vdisc65=cond((vage+0.001)>=(65+10/12),1,0) if vgebjahr==1955
replace vdisc65=cond((vage+0.001)>=(65+11/12),1,0) if vgebjahr==1956
replace vdisc65=cond((vage+0.001)>=66,1,0) if vgebjahr==1957

sort pid syear
xtset pid syear

*Generate outcome variable
gen fbirth=f.birth



eststo: xi: xtivreg2 fbirth age age2 vage vage2 i.month i.syear (vretired = vdisc60) 				if vage>=50 & vage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear) first savefirst savefprefix(_all_)
eststo: xi: xtivreg2 fbirth age age2 vage vage2 i.month i.syear (vretired = vdisc60) 				if sex==2 & vage>=50 & vage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear) first savefirst savefprefix(_daught_)
eststo: xi: xtivreg2 fbirth age age2 vage vage2 i.month i.syear (vretired = vdisc60) 				if sex==1 & vage>=50 & vage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear) first savefirst savefprefix(_son_)

eststo: xi: xtreg fbirth age age2 vage vage2 i.month i.syear vdisc60				if vage>=50 & vage<=70 & age>=25 & age<=40, fe cluster(pid) 
eststo: xi: xtreg fbirth age age2 vage vage2 i.month i.syear vdisc60 				if sex==2 & vage>=50 & vage<=70 & age>=25 & age<=40, fe cluster(pid) 
eststo: xi: xtreg fbirth age age2 vage vage2 i.month i.syear vdisc60 				if sex==1 & vage>=50 & vage<=70 & age>=25 & age<=40, fe cluster(pid) 


estout est1 est2 est3 using "${output_path}\Results\Table2.csv" , keep(vretired) cells(b(star fmt(3)) se(par(`"="("' `")""'))) mlabels(,depvars) style(fixed)  starlevels(* 0.10 ** 0.050 *** 0.010)  replace 
estout _all_vretired _daught_vretired _son_vretired using "${output_path}\Results\Table2.csv" , keep(vdisc60) cells(b(star fmt(3)) se(par(`"="("' `")""'))) mlabels(none) collabels(none)  style(fixed) starlevels(* 0.10 ** 0.050 *** 0.010)  append
estout est1 est2 est3 using "${output_path}\Results\Table2.csv" , keep() cells(none) mlabels(none) collabels(none)  style(fixed)  stats(widstat) append
estout est4 est5 est6 using "${output_path}\Results\Table2.csv" , keep(vdisc60) cells(b(star fmt(3)) se(par(`"="("' `")""'))) mlabels(,depvars) style(fixed)  starlevels(* 0.10 ** 0.050 *** 0.010)  append 
estout est1 est2 est3 using "${output_path}\Results\Table2.csv" , keep() cells(none) mlabels(none) collabels(none)  style(fixed)  stats(N) append

eststo clear


*Estimation for mothers

gen mage2=mage*mage



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






eststo: xi: xtivreg2 fbirth age age2 mage mage2 i.month i.syear (mretired = mdisc60) 				if mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear) first savefirst savefprefix(_all_)
eststo: xi: xtivreg2 fbirth age age2 mage mage2 i.month i.syear (mretired = mdisc60) 				if sex==2 & mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear) first savefirst savefprefix(_daught_)
eststo: xi: xtivreg2 fbirth age age2 mage mage2 i.month i.syear (mretired = mdisc60) 				if sex==1 & mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear) first savefirst savefprefix(_son_)

eststo: xi: xtreg fbirth age age2 mage mage2 i.month i.syear mdisc60 				if mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) 
eststo: xi: xtreg fbirth age age2 mage mage2 i.month i.syear mdisc60 				if sex==2 & mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) 
eststo: xi: xtreg fbirth age age2 mage mage2 i.month i.syear mdisc60				if sex==1 & mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) 


estout est1 est2 est3 using "${output_path}\Results\Table2.csv" , keep(mretired) cells(b(star fmt(3)) se(par(`"="("' `")""'))) mlabels(,depvars) style(fixed)  starlevels(* 0.10 ** 0.050 *** 0.010)  append 
estout _all_mretired _daught_mretired _son_mretired using "${output_path}\Results\Table2.csv" , keep(mdisc60) cells(b(star fmt(3)) se(par(`"="("' `")""'))) mlabels(none) collabels(none)  style(fixed) starlevels(* 0.10 ** 0.050 *** 0.010)  append
estout est1 est2 est3 using "${output_path}\Results\Table2.csv" , keep() cells(none) mlabels(none) collabels(none)  style(fixed)  stats(widstat) append
estout est4 est5 est6 using "${output_path}\Results\Table2.csv" , keep(mdisc60) cells(b(star fmt(3)) se(par(`"="("' `")""'))) mlabels(,depvars) style(fixed)  starlevels(* 0.10 ** 0.050 *** 0.010)  append 
estout est1 est2 est3 using "${output_path}\Results\Table2.csv" , keep() cells(none) mlabels(none) collabels(none)  style(fixed)  stats(N) append

eststo clear

