set graphics off cd "P:\EVALUATION\STATA\kh1u2\D_kh1u2" log using "logfiles\D2_kh1u2.smcl", name(D2log) replace /////////////////////////////////////////////////////////////////////////////////// // schätz-ablauf puhani und ai-norton /////////////////////////////////////// // greift zurück auf die ado-files: inteff_puhani.ado mit dem befehl inteff_puhani // und auf das ado-file probit_puhani_dd.ado für den fall der interaktion zweier dummie-variablen // die x-e müssen 01 kodiert sein (i.sex#i.m geht nicht) // do inteff_puhani.ado // do probit_puhani_dd.ado /* probit depvar d1 d2 d1d2 x, robust etc estimates store name_dd *puhani* inteff_puhani "gleiches modell s. o." if _est_name_dd==1, savedata(name, replace) savegraph1(name, replace) savegraph2(name_z, replace) *ai-norton* inteff "gleiches modell s. o." if _est_name_dd==1, savedata(name, replace) savegraph1(name, replace) savegraph2(name_z, replace) graph combine a.gph b.gph etc., col(2) saving("name",replace) graph export "name", as(eps) replace */ // Fallzahlen ermitteln: //////////////////////// /* für die Noten: ============== foreach v of varlist dm*{ gen `v'_nomiss=. replace `v'_nomiss=0 if `v'==. replace `v'_nomiss=1 if `v'~=. } contract dm72_nomiss-dm92_nomiss, freq(count) */ /* Korrektur für n=68 Fälle: dm81 (statt dm72) mit dm92 bzw. dm72 mit dm91 (statt dm92) gekennzeichnet durch dm_korr{Dummy: 0,1} */ /////////////////////////////////////////////////////////////////////////////////// global indep_var_41 /// i.sex i.m i.sex#i.m /// i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss /// i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 /// L1_mean_sozkog_nomiss /// diff_LSelbBP_kat i.diff_LSelbBP_miss /// i.BTVja i.BerufstätigkeitV_miss /// i.SchVRSAbi i.SchulabschlussV_miss /// i.e113_kateg /// i.e113_miss A_05_pc global indep_var_41_hand /// sex m sex_m kh /// PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss /// WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 /// L1_mean_sozkog_nomiss /// diff_LSelbBP_kat diff_LSelbBP_miss /// BTVja BerufstätigkeitV_miss /// SchVRSAbi SchulabschlussV_miss /// e113_kateg /// e113_miss A_05_pc global NotenD_M /// d72 m72 global Noten_meanDM /// dm72 global Indiv_f /// i.sex##i.m i.kh global Eltern /// i.BTVja i.BerufstätigkeitV_miss /// i.SchVRSAbi i.SchulabschlussV_miss /// i.e113_kateg i.e113_miss /// global A05 /// A_05_pc global BO /// PraktAnz82_kat i.PraktAnz82_miss /// PraktAnz92_kat i.PraktAnz92_miss /// i.WB_jn_91_kat i.WB_jn_91_miss /// i.Stand_BP_912 i.Stand_BP_913 global Lehrer_abs /// L1_mean_sozkog_nomiss /// i.L1_mean_sozkog_miss i.L1_einzan_d global Lehrer_diff /// L1_mean_sozkog_nomiss /// diff_LSelbBP_kat i.diff_LSelbBP_miss global leer /// global NotenSpez /// ${leer} ${Noten_meanDM} ${NotenD_M} /////////////////////////////////////////////////////////////////////////////////// // LISTE der GLOBALS: // ${Indiv_f} ${Eltern} ${A05} ${BO} ${Lehrer_abs} ${Lehrer_diff} // ${Noten_meanDM} ${NotenD_M} // ${leer} // ${NotenSpez} /* ================================================================================== Sept-2012 NEU: mit A_05*100 und Interatkionen f#m alle alte Version vom Juli-2012 in Ablage 3_2a-d.) Regression auf Noten 3_3a-d.) Notenverbesserung ================================= Aufbau in Abschnitten ================================= 3_2a. 3_3a. nur Individualdaten des Schülers 3_2b. 3_3b. mit Elternvariablen 3_2c. 3_3c. mit Einschätzung des Lehrers (Klasse 8) 3_2d. 3_3d. mit Schülervariablen (verändern sich im Zeitverlauf) und mit Einschätzung des Lehrers (Klasse 8 u 9) ? Veränderung b) Kommentar: Indiv und: Variablen der Eltern: ========================= BerufstätigkeitV3 Vater anderer Beruf BerufstätigkeitV4 Vater erlernter Beruf Ref.: arbeitslos/Hausmann ==> neu: BTVja dummy für V3uV4 SchulabschlussV1 Vater ohne SA SchulabschlussV2 Vater HS = Ref SchulabschlussV3 Vater RS SchulabschlussV4 Vater HS-Reife Abitur ==> neu: SchVRSAbi für dummy V3uV4 Ref.: <=HS e113 Alleinerziehend FamStatus12 leben zusammen/verheiratet Ref.: FamStatus34 geschieden/getrennt/verwitwet A_05 Ausländeranteil am Wohnort c) Kommentar: Indiv und Eltern sowie: Lehrervariablen: ================ L1_mean_sozkog Mittelwert über 5 Einschätzungsfragen (ist offen, zu Kompromissen bereit, findet Lsg., ist integrativ, fördert positiv das Gruppenklima - jeweils: -2 -1 +1 +2) L1_einzan_d gezielte pädagogische Einzelansprache ist wünschenswert (janein) d) Kommentar: Indiv und Eltern sowie: Lehrervariablen: ================ L1_mean_sozkog Mittelwert über 5 Einschätzungsfragen (ist offen, zu Kompromissen bereit, findet Lsg., ist integrativ, fördert positiv das Gruppenklima - jeweils: -2 -1 +1 +2) zusätzlich: =========== Anzahl Praktika in Klasse 8 und 9 Wunschberuf vorhanden janein in Klasse 9 (Zwischenzeugnis) Veränderung in der Selbständigkeit die Berufsplanung zu verantworten diff_LSelbBP == "L2_MotKomp - L1_Einzelansprache" +=+ -=- als i.diff_LSelbBPpos als Dummyvariable für positive Veränderung ca. 59% oder diff_LSelbBP_kat als -2,-1,0-1,2 kontinuierliche Variable ================================================================================== */ forvalues i=0/1 { display as text _dup(30) "=" _newline(1) as text "Untersuchung für " as result "treatment= `i'" _newline(1) as text _dup(30) "=" _newline(3) as result "nur Individualdaten" _newline(1) as text _dup(30) "=" #delimit ; regress dm92 i.sex i.m i.sex#i.m i.kh if treatment==`i' , robust ; estimates store ols_3_2a_`i'; display as text "estimates store: "as result "ols_3_2a_`i'" ; regress dm92 i.sex i.m i.sex#i.m i.kh dm72 if treatment==`i' , robust ; #delimit cr estimates store ols_3_3a_`i' display as text "estimates store: "as result "ols_3_3a_`i'" } forvalues i=0/1 { display _newline(3) display as text _dup(30) "=" _newline(1) as text "Untersuchung für " as result "treatment= `i'" _newline(1) as text _dup(30) "=" _newline(1) as result "Individualdaten und Elternvariablen" _newline(1) as text _dup(30) "=" #delimit ; regress dm92 i.sex i.m i.sex#i.m i.kh i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc if treatment==`i' , robust ; #delimit cr estimates store ols_3_2b_`i' display as text "estimates store: "as result "ols_3_2b_`i'" #delimit ; regress dm92 i.sex i.m i.sex#i.m i.kh i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc dm72 if treatment==`i' , robust coefl; #delimit cr estimates store ols_3_3b_`i' display as text "estimates store: "as result "ols_3_3b_`i'" } forvalues i=0/1 { display as text _dup(30) "=" _newline(1) as text "Untersuchung für " as result "treatment= `i'" _newline(1) as text _dup(30) "=" _newline(3) as result "Individualdaten, Elternvariablen und Einschätzung des Lehrers (Klasse 8)" _newline(1) as text _dup(30) "=" #delimit ; regress dm92 i.sex i.m i.sex#i.m i.kh i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss L1_mean_sozkog_nomiss i.L1_mean_sozkog_miss i.L1_einzan_d A_05_pc if treatment==`i' , robust ; #delimit cr estimates store ols_3_2c_`i' display as text "estimates store: "as result "ols_3_2c_`i'" #delimit ; regress dm92 i.sex i.m i.sex#i.m i.kh i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc L1_mean_sozkog_nomiss i.L1_mean_sozkog_miss i.L1_einzan_d dm72 if treatment==`i' , robust ; #delimit cr estimates store ols_3_3c_`i' display as text "estimates store: "as result "ols_3_3c_`i'" } forvalues i=0/1 { display as text _dup(30) "=" _newline(1) as text "Untersuchung für " as result "treatment= `i'" _newline(1) as text _dup(30) "=" display _newline(3) as result "Individualdaten, Elternvariablen, Einschätzung des Lehrers (Klasse 8)" _newline(1) as text _dup(30) "=" display _newline(3) as result "und Variablen zur BO in Klasse 8 und 9" _newline(1) as text _dup(30) "=" #delimit ; regress dm92 i.sex i.m i.sex#i.m i.kh i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc L1_mean_sozkog_nomiss PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 diff_LSelbBP_kat i.diff_LSelbBP_miss if treatment==`i' , robust ; #delimit cr estimates store ols_3_2d_`i' display as text "estimates store: "as result "ols_3_2d_`i'" #delimit ; regress dm92 i.sex i.m i.sex#i.m i.kh i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc L1_mean_sozkog_nomiss PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 diff_LSelbBP_kat i.diff_LSelbBP_miss dm72 if treatment==`i' , robust ; #delimit cr estimates store ols_3_3d_`i' display as text "estimates store: "as result "ols_3_3d_`i'" } ** ================================================================================== display "Verbleib nach Hauptschule (Info aus Lehrerbefragung in 92 (inkl. Schüler, falls L=.)," _newline " aufgefüllt mit Info aus nS1" _newline(2) tab verbleib3_92_pool /* ================================================================================== Sept-2012 Probit und ME - TRANSITION (3) NEU: mit A_05*100 probit Marginaler Effekt des Interaktionsterms nach: inteff_puhani inteff (Ai&Norton) alte Version vom Juli-2012 in Ablage 41 bzw. 41_b mit dm72 Verlauf: mit ANZAHL Praktikum wrszu ist raus probit_41_ÜBERGANG ==> neu: SchVRSAbi für dummy V3uV4 Ref.: ohne/HS Ref.: arbeitslos/Hausmann ==> neu: BTVja dummy für V3uV4 Berufstätigkeit NEU mit und ohne Kontrolle für das Notenviveau in (72) anstatt mit diffdm79 ==> in Ablage unter 6.) verschoben 1) probit (mit i. c. Kennzeichnung) if treatment==1, vce(cluster nr_schule) estimates store probit_41_* margins, dydx(*) post estimates store margin_41_* 2) probit ohne i. etc. if treatment ==1, vce(cluster nr_schule) inteff if _est_probit_41_*==1, savedata(41b_inteff_sex_m_`v', replace) savegraph1(41b_inteff_sex_m_`v', replace) savegraph2(41b_z_sex_m_`v', replace) 3) probit ohne i. etc. if treatment ==1, vce(cluster nr_schule) inteff_puhani if _est_probit_41_*==1, savedata(41b_inteff_pu_sex_m_`v', replace) savegraph1(41b_inteff_pu_sex_m_`v', replace) savegraph2(41b_z_pu_sex_m_`v', replace) */ // 41 * * foreach v of varlist V92pool_BV V92pool_WFS V92pool_AB{ display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_41_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; qui probit `v' i.sex i.m i.sex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_41_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "margin_41_`v'" _newline(1) /// as text _dup(50) "=" margins sex#m qui margins, dydx(*) post estimates store margin_41_`v' #delimit ; qui probit `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if treatment==1 , vce(cluster nr_schule); inteff `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if _est_probit_41_`v'==1, savedata(dta\41_inteff_sex_m_`v', replace) savegraph1(graph\41_inteff_sex_m_`v', replace) savegraph2(graph\41_z_sex_m_`v', replace) ; qui probit `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if treatment==1 , vce(cluster nr_schule); inteff_puhani `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if _est_probit_41_`v'==1, savedata(dta\41_inteff_pu_sex_m_`v', replace) savegraph1(graph\41_inteff_pu_sex_m_`v', replace) savegraph2(graph\41_z_pu_sex_m_`v', replace) ; #delimit cr display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_41b_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; qui probit `v' i.sex i.m i.sex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc dm72 if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_41b_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "margin_41b_`v'" _newline(1) /// as text _dup(50) "=" margins sex#m qui margins, dydx(*) post estimates store margin_41b_`v' #delimit ; qui probit `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if treatment==1 , vce(cluster nr_schule); inteff `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if _est_probit_41b_`v'==1, savedata(dta\41b_inteff_sex_m_`v', replace) savegraph1(graph\41b_inteff_sex_m_`v', replace) savegraph2(graph\41b_z_sex_m_`v', replace) ; qui probit `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if treatment==1 , vce(cluster nr_schule); inteff_puhani `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if _est_probit_41b_`v'==1, savedata(dta\41b_inteff_pu_sex_m_`v', replace) savegraph1(graph\41b_inteff_pu_sex_m_`v', replace) savegraph2(graph\41b_z_pu_sex_m_`v', replace) ; #delimit cr graph combine graph\41_inteff_sex_m_`v'.gph graph\41b_inteff_sex_m_`v'.gph graph\41_z_sex_m_`v'.gph graph\41b_z_sex_m_`v'.gph, col(2) title("Ai & Norton (right side: with dm72") saving("graph\41ub_inteff_sex_m_`v'", replace) graph export "graph\41ub_inteff_sex_m_`v'.eps", as(eps) replace graph combine graph\41_inteff_pu_sex_m_`v'.gph graph\41b_inteff_pu_sex_m_`v'.gph graph\41_z_pu_sex_m_`v'.gph graph\41b_z_pu_sex_m_`v'.gph, col(2) title("Puhani(right side: with dm72") saving("graph\41ub_inteff_pu_sex_m_`v'", replace) graph export "graph\41ub_inteff_pu_sex_m_`v'.eps", as(eps) replace } * * * * // DEUTSCH MATHE 51 // mit d72 und m72 in der Spezifikation zur Kontrolle des Ausgangsniveaus zu // Beginn der Maßnahme * * foreach v of varlist V92pool_BV V92pool_WFS V92pool_AB{ #delimit ; probit `v' $indep_var_41 $NotenD_M if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_51_`v' margins sex#m margins, dydx(*) post estimates store margin_51_`v' #delimit ; qui probit `v' $indep_var_41_hand $NotenD_M if treatment==1 , vce(cluster nr_schule); inteff `v' $indep_var_41_hand $NotenD_M if _est_probit_51_`v'==1, savedata(dta\51_inteff_sex_m_`v', replace) savegraph1(graph\51_inteff_sex_m_`v', replace) savegraph2(graph\51_z_sex_m_`v', replace) ; qui probit `v' $indep_var_41_hand $NotenD_M if treatment==1 , vce(cluster nr_schule); inteff_puhani `v' $indep_var_41_hand $NotenD_M if _est_probit_51_`v'==1, savedata(dta\51_inteff_pu_sex_m_`v', replace) savegraph1(graph\51_inteff_pu_sex_m_`v', replace) savegraph2(graph\51_z_pu_sex_m_`v', replace) ; #delimit cr graph combine "graph\51_inteff_sex_m_`v'.gph" "graph\51_z_sex_m_`v'.gph" , col(2) title("Ai & Norton (right side: with d+m") saving("graph\51_inteff_sex_m_`v'", replace) graph export "graph\51_inteff_sex_m_`v'.eps", as(eps) replace graph combine "graph\51_inteff_pu_sex_m_`v'.gph" "graph\51_z_pu_sex_m_`v'.gph" , col(2) title("Puhani(right side: with d+m") saving("graph\51_inteff_pu_sex_m_`v'", replace) graph export "graph\51_inteff_pu_sex_m_`v'.eps", as(eps) replace } * * // Interaktionseffekte mit margins # * * * foreach v of varlist V92pool_BV V92pool_WFS V92pool_AB{ probit `v' $indep_var_41 if treatment==1 ,vce(cluster nr_schule) estimates store probit_41_`v' margins sex#m margins, dydx(*) post estimates store margin_41_`v' qui probit `v' $indep_var_41 if treatment==1 ,vce(cluster nr_schule) display _newline(2) as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as text "Transition into:" as result "`v'" _newline(2) /// as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as result "male as female vs. male as male " _newline(2) /// as text _dup(10) " " as result "female as male vs. female as female" _newline(1) /// as text _dup(50) "=" margins, over(sex) at(sex=(0 1)) contrast(atcontrast(r._at) wald) margins, over(sex) at(sex=(0 1)) post margins, coeflegend test _b[2._at#0bn.sex] = _b[1bn._at#0bn.sex] // male as female vs. male as male test _b[2._at#1.sex] = _b[1bn._at#1.sex] // female as male vs. female as female qui probit `v' $indep_var_41 if treatment==1 ,vce(cluster nr_schule) display _newline(2) as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as text "Transition into:" as result "`v'" _newline(2) /// as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as result "ger as migr. vs. ger as ger " _newline(2) /// as text _dup(10) " " as result "migr. as ger vs. migr. as migr. " _newline(1) /// as text _dup(50) "=" margins, over(m) at(m=(0 1)) contrast(atcontrast(r._at) wald) margins, over(m) at(m=(0 1)) post margins, coeflegend test _b[2._at#0bn.m] = _b[1bn._at#0bn.m] // ger as migr. vs. ger as ger test _b[2._at#1.m] = _b[1bn._at#1.m] // migr. as ger vs. migr. as migr. qui probit `v' $indep_var_41 if treatment==1 ,vce(cluster nr_schule) display _newline(2) as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as text "Transition into:" as result "`v'" _newline(2) /// as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as result "migr as female vs migr as male" _newline(2) /// as text _dup(10) " " as result "ger as female vs ger as male " _newline(1) /// as text _dup(50) "=" margins, over(m) at(sex=(0 1)) contrast(atcontrast(r._at) wald) margins, over(m) at(sex=(0 1)) post margins, coeflegend test _b[2._at#1.m] = _b[1bn._at#1.m] // migr as female vs migr as male test _b[2._at#0bn.m] = _b[1bn._at#0bn.m] // ger as female vs ger as male qui probit `v' $indep_var_41 if treatment==1 ,vce(cluster nr_schule) display _newline(2) as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as text "Transition into:" as result "`v'" _newline(2) /// as text _dup(50) "=" _newline(1) /// as text _dup(10) " " as result "male as ger vs. male as migr" _newline(2) /// as text _dup(10) " " as result "female as ger vs. female as migr" _newline(1) /// as text _dup(50) "=" margins, over(sex) at(m=(0 1)) contrast(atcontrast(r._at) wald) margins, over(sex) at(m=(0 1)) post margins, coeflegend test _b[2._at#0bn.sex] = _b[1bn._at#0bn.sex] // male as ger vs. male as migr test _b[2._at#1.sex] = _b[1bn._at#1.sex] // female as ger vs. female as migr } ************************************************************************************************* *** zusätzliche Spezifikationen *** // SEX=MÄNNLICH // malesex 61 * * foreach v of varlist V92pool_BV V92pool_WFS V92pool_AB{ display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_61_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; qui probit `v' i.malesex i.m i.malesex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_61_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "margin_61_`v'" _newline(1) /// as text _dup(50) "=" margins malesex#m qui margins, dydx(*) post estimates store margin_61_`v' #delimit ; qui probit `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if treatment==1 , vce(cluster nr_schule); inteff `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if _est_probit_61_`v'==1, savedata(dta\61_inteff_malesex_m_`v', replace) savegraph1(graph\61_inteff_malesex_m_`v', replace) savegraph2(graph\61_z_malesex_m_`v', replace) ; qui probit `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if treatment==1 , vce(cluster nr_schule); inteff_puhani `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc if _est_probit_61_`v'==1, savedata(dta\61_inteff_pu_malesex_m_`v', replace) savegraph1(graph\61_inteff_pu_malesex_m_`v', replace) savegraph2(graph\61_z_pu_malesex_m_`v', replace) ; #delimit cr display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_61b_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; qui probit `v' i.malesex i.m i.malesex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05_pc dm72 if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_61b_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "margin_61b_`v'" _newline(1) /// as text _dup(50) "=" margins malesex#m qui margins, dydx(*) post estimates store margin_61b_`v' #delimit ; qui probit `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if treatment==1 , vce(cluster nr_schule); inteff `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if _est_probit_61b_`v'==1, savedata(dta\61b_inteff_malesex_m_`v', replace) savegraph1(graph\61b_inteff_malesex_m_`v', replace) savegraph2(graph\61b_z_malesex_m_`v', replace) ; qui probit `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if treatment==1 , vce(cluster nr_schule); inteff_puhani `v' malesex m malesex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05_pc dm72 if _est_probit_61b_`v'==1, savedata(dta\61b_inteff_pu_malesex_m_`v', replace) savegraph1(graph\61b_inteff_pu_malesex_m_`v', replace) savegraph2(graph\61b_z_pu_malesex_m_`v', replace) ; #delimit cr graph combine "graph\61_inteff_malesex_m_`v'.gph" "graph\61b_inteff_malesex_m_`v'.gph" "graph\61_z_malesex_m_`v'.gph" "graph\61b_z_malesex_m_`v'.gph", col(2) title("Ai & Norton (right side: with dm72") saving("graph\61ub_inteff_malesex_m_`v'", replace) graph export "graph\61ub_inteff_malesex_m_`v'.eps", as(eps) replace graph combine "graph\61_inteff_pu_malesex_m_`v'.gph" "graph\61b_inteff_pu_malesex_m_`v'.gph" "graph\61_z_pu_malesex_m_`v'.gph" "graph\61b_z_pu_malesex_m_`v'.gph", col(2) title("Puhani(right side: with dm72") saving("graph\61ub_inteff_pu_malesex_m_`v'", replace) graph export "graph\61ub_inteff_pu_malesex_m_`v'.eps", as(eps) replace } ********************************************************************** * * * * 44_ und 44b_ (mit dm72) * * mit Interaktion zwischen A_05 und Migrationshintergrund m: A_05#i.m * ========================================================== * // INTERAKTIONSEFFEKT VON M IST FALSCH!!!!!! * inteff einmal mit A_05 m A_05_m vorne * und einmal mit sex m sex_m vorne * * * foreach v of varlist V92pool_BV V92pool_WFS V92pool_AB{ display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_44_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; qui probit `v' i.sex i.m i.sex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss c.A_05#i.m A_05 if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_44_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "margin_44_`v'" _newline(1) /// as text _dup(50) "=" margins sex#m qui margins, dydx(*) post estimates store margin_44_`v' #delimit ; qui probit `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05 A_05_m if treatment==1 , vce(cluster nr_schule); inteff `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05 A_05_m if _est_probit_44_`v'==1, savedata(dta\44_inteff_sex_m_`v', replace) savegraph1(graph\44_inteff_sex_m_`v', replace) savegraph2(graph\44_z_sex_m_`v', replace) ; inteff `v' A_05 m A_05_m sex sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss if _est_probit_44_`v'==1, savedata(dta\44_inteff_A_05_m_`v', replace) savegraph1(graph\44_inteff_A_05_m_`v', replace) savegraph2(graph\44_z_A_05_m_`v', replace) ; #delimit cr display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_44b_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; qui probit `v' i.sex i.m i.sex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss c.A_05#i.m A_05 dm72 if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_44b_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "margin_44b_`v'" _newline(1) /// as text _dup(50) "=" margins sex#m qui margins, dydx(*) post estimates store margin_44b_`v' #delimit ; qui probit `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05 A_05_m dm72 if treatment==1 , vce(cluster nr_schule); inteff `v' sex m sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss A_05 A_05_m dm72 if _est_probit_41b_`v'==1, savedata(dta\44b_inteff_sex_m_`v', replace) savegraph1(graph\44b_inteff_sex_m_`v', replace) savegraph2(graph\44b_z_sex_m_`v', replace) ; inteff `v' A_05 m A_05_m sex sex_m kh PraktAnz82_kat PraktAnz82_miss PraktAnz92_kat PraktAnz92_miss WB_jn_91_kat WB_jn_91_miss Stand_BP_912 Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat diff_LSelbBP_miss BTVja BerufstätigkeitV_miss SchVRSAbi SchulabschlussV_miss e113_kateg e113_miss dm72 if _est_probit_44_`v'==1, savedata(44b_inteff_A_05_m_`v', replace) savegraph1(44b_inteff_A_05_m_`v', replace) savegraph2(44b_z_A_05_m_`v', replace) ; #delimit cr graph combine 44_inteff_sex_m_`v'.gph 44b_inteff_sex_m_`v'.gph 44_z_sex_m_`v'.gph 44b_z_sex_m_`v'.gph, col(2) saving("graph\44ub_inteff_sex_m_`v'", replace) graph export "graph\44ub_inteff_sex_m_`v'.eps", as(eps) replace graph combine 44_inteff_A_05_m_`v'.gph 44b_inteff_A_05_m_`v'.gph 44_z_A_05_m_`v'.gph 44b_z_A_05_m_`v'.gph, col(2) saving("graph\44ub_inteff_A_05_m_`v'", replace) graph export "graph\44ub_inteff_A_05_m_`v'.eps", as(eps) replace } * * * * * * * /*************************************************************** * 4_2 mit Noten zu unterschiedlichen Zeitpunkten => als pooled * OLS ausarbeiten (mehr Verlauf) *=============================================================== * mit Veränderung in der Selbständigkeit der Berufsplanung... * keine beseren Ergebnisse (drop) foreach v of varlist V92pool_BV V92pool_WFS V92pool_AB{ foreach note of numlist 72 81 82 91 92{ display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_42_`v'_`note'" _newline(1) /// as text _dup(50) "=" #delimit ; qui probit `v' i.sex i.m i.sex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss i.L1_Steigerung1 i.L1_Steigerung3 i.L1_Steigerung4 i.L1_Steigerung_miss i.L2_MotKomp3 i.L2_MotKomp4u5 i.L2_MotKomp_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss dm`note' if treatment==1 , robust ; #delimit cr estimates store probit_42_`v'_`note' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "probit_marg_42_`v'_`note'" _newline(1) /// as text _dup(50) "=" margins sex#m qui margins, dydx(*) post estimates store margin_42_`v'_`note' } } esttab *42_*BV_* , /// /*using "tex\results_probit.tex", append */ /// /*title(Regressionen auf Noten\lab{reg_32})*/ /// /*booktabs onecell */ /// nogaps nolines /// drop(*0b.*) se(%9.4f) brackets /// pr2 /*r2 ar2 */ /// alignment(D{.}{.}{-1}) /// star(* 0.10 ** 0.05 *** 0.01) /// /*label*/ /*oder coefflabel()*/ /// /*coeflabel(1.sex#1.m "sex*m" 1.sex#1.m#1.kh "sex*m*kh")*/ /// addnote("") esttab *42_*WFS_* , /// /*using "tex\results_probit.tex", append */ /// /*title(Regressionen auf Noten\lab{reg_32})*/ /// /*booktabs onecell */ /// nogaps nolines /// drop(*0b.*) se(%9.4f) brackets /// pr2 /*r2 ar2 */ /// alignment(D{.}{.}{-1}) /// star(* 0.10 ** 0.05 *** 0.01) /// /*label*/ /*oder coefflabel()*/ /// /*coeflabel(1.sex#1.m "sex*m" 1.sex#1.m#1.kh "sex*m*kh")*/ /// addnote("") esttab *42_*AB_* , /// /*using "tex\results_probit.tex", append */ /// /*title(Regressionen auf Noten\lab{reg_32})*/ /// /*booktabs onecell */ /// nogaps nolines /// drop(*0b.*) se(%9.4f) brackets /// pr2 /*r2 ar2 */ /// alignment(D{.}{.}{-1}) /// star(* 0.10 ** 0.05 *** 0.01) /// /*label*/ /*oder coefflabel()*/ /// /*coeflabel(1.sex#1.m "sex*m" 1.sex#1.m#1.kh "sex*m*kh")*/ /// addnote("") *****************************************************************/ /***************************************************************** * Übergänge in "further schooling" differenziert: nach WFS und BFS * * ohne ME nach Ai und Norton 2003 *================================================================= foreach v of varlist V92pool_WR { display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_43_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; probit `v' i.sex i.m i.sex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg A_05 diffdm79 if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_43_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "probit_marg_43_`v'" _newline(1) /// as text _dup(50) "=" margins sex#m qui margins, dydx(*) post estimates store margin_43_`v' } foreach v of varlist V92pool_BFS { display _newline(2) as text _dup(50) "=" _newline(1) /// as text "estimates store: " as result "probit_43_`v'" _newline(1) /// as text _dup(50) "=" #delimit ; probit `v' i.sex i.m i.sex#i.m i.kh PraktAnz82_kat i.PraktAnz82_miss PraktAnz92_kat i.PraktAnz92_miss i.WB_jn_91_kat i.WB_jn_91_miss i.Stand_BP_912 i.Stand_BP_913 L1_mean_sozkog_nomiss diff_LSelbBP_kat i.diff_LSelbBP_miss i.BTVja i.BerufstätigkeitV_miss i.SchVRSAbi i.SchulabschlussV_miss i.e113_kateg i.e113_miss A_05 diffdm79 if treatment==1 , vce(cluster nr_schule); #delimit cr estimates store probit_43_`v' display _newline(2) as text _dup(50) "=" _newline(1) /// as text "margins store: " as result "probit_marg_43_`v'" _newline(1) /// as text _dup(50) "=" margins sex#m qui margins, dydx(*) post estimates store margin_43_`v' } esttab probit_43_* /// using "results_probit43.csv", replace /// /// title(Probit estimation\lab{prob_41}) /// /// booktabs onecell /// nogaps nolines /// drop(*0b.*) se(%9.4f) brackets /// pr2 /*r2 ar2 */ /// alignment(D{.}{.}{-1}) /// star(* 0.10 ** 0.05 *** 0.01) /// /*label*/ /*oder coefflabel()*/ /// coeflabel(1.sex#1.m "female*migration background" 1.sex "female" 1.m "migration background" 1.kh "cohort 2" /// PraktAnz82_kat "number of internships (8)" /// 1.PraktAnz82_miss "number of internships (8) missing" /// PraktAnz92_kat "number of internships (9)" /// 1.PraktAnz92_miss "number of internships (9) missing" /// 1.WB_jn_91_kat "knows desired occupation (91)" /// 1.WB_jn_91_miss "knows desired occupation (91) missing" /// 1.Stand_BP_912 "state of career planning: moderate improvement" /// 1.Stand_BP_913 "state of career planning: substantial improvement" /// L1_mean_sozkog_nomiss "social and kognitive skills (82)" /// diff_LSelbBP_kat "career planning: improvement in autonomy (82-92)" /// 1.diff_LSelbBP_miss "teacher response missing" /// 1.BTVja "father employed" 1.BerufstätigkeitV_miss "father employed missing" /// 1.SchVRSAbi "father's degree of schooling: medium-high" /// 1.SchulabschlussV_miss "father's degree of schooling missing" /// 1.e113_kateg "single-parent family" /// 1.e113_miss "single-parent family missing" /// A_05 "share of foreigners in neighborhood" /// diffdm79 "improvement in grades (72-92)" /// _cons "constant") /// addnote("") *****************************************************************/ log close D2log * ======================================================================