/* www.felix-bittmann.de Version 1.0 November 2017 */ clear all //Speicher leeren, Daten loeschen *version 12 //Eigene Stata Version angeben *set more off *Datei oeffnen* use "Allbus2014.dta", clear //Pfad anpassen *V275 koerpergroesse *V81 geschlecht *V849 attraktivitaet fremdeinschaetzung *V84 alter *V834 nummer interviewer *V836 alter interviewer *V417 einkommen offene abfrage tab1 V275 V81 V849 V84 V834 ***Einkommen Log*** gen log_einkommen = log(V417) label variable log_einkommen "Einkommen log" ***Bildung rekodieren recode V86 (1 2 = 0 "HS") ( 3=1 "RS") (4 5 = 2 "ABI") (else = .), gen(bildung) label variable bildung "Hoechster Bildungsabschluss" tab bildung *Personen mit fehlenden Werten entfernen* foreach var of varlist V275 V81 V849 V84 V417 bildung { drop if missing(`var') } *describe *********************************** *Deskriptive Statistiken *********************************** histogram V849, discrete histogram V275, by(V81) histogram V417 histogram log_einkommen histogram V84 ttest V275, by(V7) //Ost VS West *********************************** *Normale Regression //Ohne Multilevel reg V849 i.V81 c.V275 c.V84 c.log_einkommen i.bildung *********************************** *******Multilevel Regression******* *********************************** *1. Mann befragt Mann xtmixed V849 c.V275 c.V84 c.log_einkommen i.bildung || V834: /// if V835== 1 & V81 == 1 *2. Mann befragt Frau xtmixed V849 c.V275 c.V84 c.log_einkommen i.bildung || V834: /// if V835== 1 & V81 == 2 *3. Frau befragt Frau xtmixed V849 c.V275 c.V84 c.log_einkommen i.bildung || V834: /// if V835== 2 & V81 == 2 *4. Frau befragt Mann xtmixed V849 c.V275 c.V84 c.log_einkommen i.bildung || V834: /// if V835== 2 & V81 == 1 estat icc //Intraclass Korrelation ***Diagnostik*** *Direkt nach Modell ausfuehren predict resid, residuals predict resid_std, rstandard qnorm resid_std //Normalverteilung der Residuen *********************************************************** *BONUS: Likelihood-Ratio (LR) Test: *Modell ohne zusaetzliche erklarende Variable berechnen und Ergebnisse speichern: xtmixed V849 c.V84 c.log_einkommen i.bildung || V834: /// if V835== 2 & V81 == 1 est store m1 *Modell mit zusaetzlicher Variable berechnen: xtmixed V849 c.V275 c.V84 c.log_einkommen i.bildung || V834: /// if V835== 2 & V81 == 1 est store m2 lrtest m1 m2, stats //LR Test *Ist das Ergebnis hochsignifikant, verbessert die zuaetzliche Variable das Modell. *Niedrigere Werte bei AIC und BIC deuten ebenfalls auf einen besseren Modellfit hin.