cd "${output_path}\tmp\"



use "Linked_data.dta", clear






gen age2=age*age


gen mage2=mage*mage


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 vage2=vage*vage


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 mvretired=cond(mretired==1 & vretired==1,1,cond(mretired!=. & vretired!=.,0,.))
gen mvdisc60=cond(vdisc60==1 & mdisc60==1,1,cond(vdisc60!=. & mdisc60!=.,0,.))

sort pid syear
gen fbirth=f.birth

eststo: xi: xtivreg2 fbirth age age2 vage vage2 mage mage2 i.month i.syear (vretired mretired mvretired = vdisc60 mdisc60 mvdisc60) if vage>=50 & vage<=70 & mage>=50 & mage<=70 & age>=25 & age<=40, fe first cluster(pid) partial(i.month i.syear)
eststo: xi: xtivreg2 fbirth age age2 vage vage2 mage mage2 i.month i.syear (vretired mretired mvretired = vdisc60 mdisc60 mvdisc60) if sex==2 & vage>=50 & vage<=70 & mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear)
eststo: xi: xtivreg2 fbirth age age2 vage vage2 mage mage2 i.month i.syear (vretired mretired mvretired = vdisc60 mdisc60 mvdisc60) if sex==1 & vage>=50 & vage<=70 & mage>=50 & mage<=70 & age>=25 & age<=40, fe cluster(pid) partial(i.month i.syear)

estout _all using "${output_path}\Results\TableA8.csv" , keep(vretired mretired mvretired) cells(b(star fmt(3)) se(par(`"="("' `")""'))) mlabels(,depvars) style(fixed)  stats(widstat N) ///
starlevels(* 0.10 ** 0.050 *** 0.010)  replace
eststo clear





