|
| 1 | +/************************************************************************************** |
| 2 | +* |
| 3 | +* PROGRAM TO ANALYSE EMPLOYMENT TRANSATIONS |
| 4 | +* |
| 5 | +* Last version: Justin van de Ven, 06 May 2025 |
| 6 | +* First version: Justin van de Ven, 06 May 2025 |
| 7 | +* |
| 8 | +**************************************************************************************/ |
| 9 | + |
| 10 | +clear all |
| 11 | +global moddir = "C:\Justin\dev\CeMPASimPaths\output\20251029155614\csv" |
| 12 | +global outdir = "C:\Justin\analysis\data\" |
| 13 | +cd "$outdir" |
| 14 | + |
| 15 | + |
| 16 | +/************************************************************************************** |
| 17 | +* start |
| 18 | +**************************************************************************************/ |
| 19 | +global year_start = 2015 |
| 20 | +global year_end = 2026 |
| 21 | + |
| 22 | + |
| 23 | +/************************************************************************************** |
| 24 | +* load data |
| 25 | +**************************************************************************************/ |
| 26 | +import delimited using "$moddir/BenefitUnit.csv", clear |
| 27 | +rename *, l |
| 28 | +rename id_benefitunit idbenefitunit |
| 29 | +gsort idbenefitunit time |
| 30 | +save "$outdir/temp0", replace |
| 31 | +import delimited using "$moddir/Person.csv", clear |
| 32 | +rename *, l |
| 33 | +rename id_person idperson |
| 34 | +rename socialcareprovision socialcareprovision_p |
| 35 | +gsort idbenefitunit time idperson |
| 36 | +merge m:1 idbenefitunit time using temp0 |
| 37 | +gsort time idbenefitunit idperson |
| 38 | +gen refbenefitunit = 0 |
| 39 | +replace refbenefitunit = 1 if (idbenefitunit != idbenefitunit[_n-1]) |
| 40 | + |
| 41 | +destring hoursworkedweekly, replace force |
| 42 | +recode hoursworkedweekly (missing=0) |
| 43 | +gen idNotEmployedAdult = (hoursworkedweekly<0.1 & dag>17) |
| 44 | + |
| 45 | +gen led = (deh_c3=="Low") |
| 46 | +gen med = (deh_c3=="Medium") |
| 47 | +gen hed = (deh_c3=="High") |
| 48 | + |
| 49 | +gen male = (dgn=="Male") |
| 50 | + |
| 51 | +gen idbaby = (dag==0) |
| 52 | +bys time idbenefitunit: egen baby = sum(idbaby) |
| 53 | + |
| 54 | +save "$outdir/temp1", replace |
| 55 | + |
| 56 | + |
| 57 | +/************************************************************************************** |
| 58 | +* analysis |
| 59 | +**************************************************************************************/ |
| 60 | +use "$outdir/temp1", clear |
| 61 | +global year_ref = $year_start-1 |
| 62 | + |
| 63 | +// block 1 statistics |
| 64 | +gen emp = (hoursworkedweekly>0.5) |
| 65 | +gen emp_to_nemp = 0 |
| 66 | +gen nemp_to_emp = 0 |
| 67 | +gsort idperson time |
| 68 | +replace emp_to_nemp = 1 if (idperson[_n-1]==idperson & time[_n-1]+1==time & emp[_n-1]==1 & emp==0) |
| 69 | +replace nemp_to_emp = 1 if (idperson[_n-1]==idperson & time[_n-1]+1==time & emp[_n-1]==0 & emp==1) |
| 70 | +order time idperson hoursworkedweekly emp emp_to_nemp |
| 71 | +matrix store1 = J($year_end-$year_ref,1,.) |
| 72 | +forvalues yy = $year_start/$year_end { |
| 73 | + qui{ |
| 74 | + sum baby if (time==`yy' & dag>17 & dag<65), mean |
| 75 | + mat store1[`yy'-${year_ref},1] = r(mean) |
| 76 | + } |
| 77 | +} |
| 78 | +matlist store1 |
0 commit comments