*Do file für Blogeintrag 21.03.2020
version 15.1 // Stata/MP 15.1 for Windows

use pequiv_short.dta, clear
merge 1:1 pid syear using "$Ordner\pl.dta", keepusing(pli0080 plh0252 plj0078)
replace pli0080=. if pli0080<0
replace plh0252=. if plh0252<0

recode pli0080 (1=5 "täglich") (2=4 "wöchentlich") (3=3 "monatlich") (4=2 "seltener") (5=1 "nie"), gen(freunde_bek_treffen)

recode plh0252 (1=1 "stimme überhaupt nicht zu") (2=2) (3=3) (4=4) (5=5) (6=6) (7=7 "7 stimme voll zu"), gen(wenig_kontrolle_über_leben)
label var wenig_kontrolle_über_leben "Wenig Kontrolle über Leben"

replace plj0078=. if plj0078<0

recode plj0078  (1=5 "voll und ganz") (2=4 "überwiegend") (3=3 "in mancher Beziehung") (4=2 "kaum") (5=1 "gar nicht"), gen(deutsch_fuehlen)

sort pid syear
save pequiv_temp.dta, replace
use pequiv_temp, clear

* Freunde / Bekannte treffen
eststo clear
eststo fe: xtreg life_sat i.pgmode i.age_cat5 i.duration5 ib5.freunde_bek_treffen , fe vce(robust)
coefplot , keep(*.freunde_bek_treffen) legend(region(color(%0)) col(1) ring(0) position(3)) ylabel() nooffset xtitle(Freunde und Bekannte treffen) sort(,descending) ytitle(Verringerung Lebenszufriedenheit 0-100) $coefplotopts

*Kontrolle über Leben
eststo clear
eststo fe_plh0252: xtreg life_sat i.pgmode i.duration5 i.age_cat5 i.wenig_kontrolle_über_leben, fe vce(robust)
coefplot , keep(*.wenig_kontrolle_über_leben) legend(region(color(%0)) col(1) ring(0) position(3)) ylabel() nooffset xtitle("Zustimmung: Ich habe wenig Kontrolle über die Dinge, die in meinem Leben passieren") sort(,descending) ytitle(Verringerung Lebenszufriedenheit 0-100) xlabel(,alternate) scale(.9) $coefplotopts

* Gemeinschaftsgefühl
eststo clear
eststo FE: xtreg life_sat i.pgmode i.duration5 i.age_cat5 i.syear i.deutsch_fuehlen, fe vce(robust) 
coefplot FE, keep(*.deutsch_fuehlen)  xtitle("Zustimmung: Sich als Deutscher fühlen", margin(0 0 0 3)) nooffset xlabel(,alternate) ylabel(0(1)6) ytitle("Zufriedenheitszugewinn 0-100") $coefplotopts