* Autor: Holger Langkabel (Grundlage: gen_p50p - 5.2 - 1R02 - diff-in-diff niea) 09-05-2012 global settings "R:\C-Projekte\C07-Perspektive50plus\DO - 0 - Master & Settings\gen_p50p - 0.0 - settings.do" include "$settings" global settings_analyse "R:\C-Projekte\C07-Perspektive50plus\DO - 6.x - Analyse\gen_p50p - 6.0 - settings analyse.do" include "$settings_analyse" global days 30(30)540 // <-- bestimmt Tagesintervall f𲠓ch嵺ungen capture log close log using "$logdir\gen_p50p - 6.2 - 1R05 - diff-in-diff iea.log", replace cd "$analysedir" use "$proddir\p50plus_diff_in_diff.dta", clear * Integrationswahrscheinlichkeiten ----------------------------------------------------------- foreach dep in iea lb { preserve bysort t1 d: summarize `dep'_* foreach days of numlist $days { summarize `dep'_`days' generate `dep'_mean_`days' = `r(mean)' summarize `dep'_`days' if d == 0 & t1 == 0 generate `dep'_mean_control_t0_`days' = `r(mean)' summarize `dep'_`days' if d == 1 & t1 == 0 generate `dep'_mean_treat_t0_`days' = `r(mean)' summarize `dep'_`days' if d == 0 & t1 == 1 generate `dep'_mean_control_t1_`days' = `r(mean)' summarize `dep'_`days' if d == 1 & t1 == 1 generate `dep'_mean_treat_t1_`days' = `r(mean)' } keep *_mean* keep if _n == 1 generate aux=1 reshape long `dep'_mean_ `dep'_mean_control_t0_ `dep'_mean_treat_t0_ `dep'_mean_control_t1_ `dep'_mean_treat_t1_ , i( aux ) j( days ) drop aux twoway connected `dep'_mean_control_t0_ `dep'_mean_treat_t0 days graph export "DiD-Signifikanz `dep'\mean t0.png", replace twoway connected `dep'_mean_control_t1_ `dep'_mean_treat_t1 days graph export "DiD-Signifikanz `dep'\mean t1.png", replace save "DiD-Signifikanz `dep'\mean.dta", replace restore * Conditional DiD ------------------------------------------------------------------------------ preserve drop if P50 == 0 levelsof wz_ieb, local(wz_levels) levelsof berufsgruppe, local(beruf_levels) foreach days of numlist $days { display _newline(3) as result "`dep'_`days'" probit `dep'_`days' t1##d /*alol1*/ /*alol24*/ /*alol57*/ bel1 bel24 bel57 prgl1 /*prgl24 prgl57*/ /*olfl1 olfl24 olfl57*/ ib6.ausb_k /*i.iaw_ausb_ieb*/ /*i.gro_id_lhg*/ /*azw_id_lhg*/ lhg_bg_anz_lhg /*lhg_ehb_anz_lhg*/ /*lhg_m_ehb_anz_lhg*/ /*lhg_nef_anz_lhg*/ /*lhg_u15_anz_lhg*/ /*lhg_a65_anz_lhg*/ /*ib1.b_abz_id_such_asu*/ /*i.b_jan_id_cberrueck_asu*/ i.b_jan_id_ges_ein_asu ib1.bsb_k /*ib4.mig_gruppe_ieb*/ i.ieb_regio_typ09_int_lhg ib3.sbs_k i.sex_ieb i.iaw_sgb_ii_tr_art /*i.staat_k*/ i.wz_ieb /*i.ieb_sgb_ii_tr_dst_lhg*/ i.berufsgruppe /*i.bkz_ausb_ba_schl_k*/ i.berufstellg_ieb /*des_bip_et*/ des_bev_a_ant /*des_wkg_fort des_wkg_zu*/ /*des_wkg_saldo des_bevdichte*/ /*i.regionaldirektion*/ /*des_an_cp_ant*/ /*des_bev_50p*/ /*alter*/ i.altersgruppe i.iaw_mig_gruppe /*des_alo*/, cluster(IDTDS) summarize `e(depvar)' generate intprob_`days' = `r(mean)' summarize p50plus if e(sample) == 1 generate p50plus_`days' = `r(mean)' generate n_`days' = e(N) test _cons generate cons_p_`days' = `r(p)' generate cons_c_`days' = _b[_cons] generate cons_u_`days' = _b[_cons] + invnormal( 1 - $siglev / 2 ) * _se[_cons] generate cons_l_`days' = _b[_cons] - invnormal( 1 - $siglev / 2 ) * _se[_cons] testparm 1.d generate d_p_`days' = `r(p)' generate d_c_`days' = _b[1.d] generate d_u_`days' = _b[1.d] + invnormal( 1 - $siglev / 2 ) * _se[1.d] generate d_l_`days' = _b[1.d] - invnormal( 1 - $siglev / 2 ) * _se[1.d] testparm 1.t1 generate t1_p_`days' = `r(p)' generate t1_c_`days' = _b[1.t1] generate t1_u_`days' = _b[1.t1] + invnormal( 1 - $siglev / 2 ) * _se[1.t1] generate t1_l_`days' = _b[1.t1] - invnormal( 1 - $siglev / 2 ) * _se[1.t1] testparm t1#d generate t1d_p_`days' = `r(p)' generate t1d_c_`days' = _b[1.t1#1.d] generate t1d_u_`days' = _b[1.t1#1.d] + invnormal( 1 - $siglev / 2 ) * _se[1.t1#1.d] generate t1d_l_`days' = _b[1.t1#1.d] - invnormal( 1 - $siglev / 2 ) * _se[1.t1#1.d] } * table d t1 if e(sample) == 1 keep *_p_* *_c_* *_u_* *_l_* p50plus_* intprob_* n_* keep if _n == 1 generate aux=1 rename *_p_30 *_p_ describe *_p_, varlist local pvars `r(varlist)' display "`pvars'" rename *_p_ *_p_30 rename *_c_30 *_c_ rename *_l_30 *_l_ rename *_u_30 *_u_ describe *_c_ *_u_ *_l_, varlist local coeffvars `r(varlist)' display "`coeffvars'" rename *_c_ *_c_30 rename *_l_ *_l_30 rename *_u_ *_u_30 reshape long `pvars' `coeffvars' p50plus_ intprob_ n_ , i( aux ) j( days ) drop aux save "DiD-Signifikanz `dep'\conditinal did probit `dep'.dta", replace restore } * profiler off profiler report log close