/* create status employment */ /* creating a dummy variable to indicate worksers who are self-employed, employers, and family workers */ recode k6r26 ( 1/3 = 0 ) ( 6 = 0 ) gen status = k6r26 == 0 replace work = work* status keep k1* rhead mnoart age sex mwert agr_income nagr_income wk_wage work /* generating variables indicating the number of self-employed workers in each age group by household */ egen p1014 = count(age) if age >= 10 & age <= 14 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p1519 = count(age) if age >= 15 & age <= 19 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p2024 = count(age) if age >= 20 & age <= 24 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p2529 = count(age) if age >= 25 & age <= 29 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p3034 = count(age) if age >= 30 & age <= 34 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p3539 = count(age) if age >= 35 & age <= 39 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p4044 = count(age) if age >= 40 & age <= 44 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p4549 = count(age) if age >= 45 & age <= 49 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p5054 = count(age) if age >= 50 & age <= 54 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p5559 = count(age) if age >= 55 & age <= 59 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen p6099 = count(age) if age >= 60 & age <= 99 & work == 1 , by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) /* generating total household self-employment income */ egen nl_income = rsum(agr_income nagr_income) recode p* ( . = 0 ) /* labor component of household self-employment income */ egen nl_income23 = (2/3) * nl_income egen mp1014 = max(p1014), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp1519 = max(p1519), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp2024 = max(p2024), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp2529 = max(p2529), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp3034 = max(p3034), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp3539 = max(p3539), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp4044 = max(p4044), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp4549 = max(p4549), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp5054 = max(p5054), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp5559 = max(p5559), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) egen mp6099 = max(p6099), by(k1r1 k1r2 k1r3 k1r4 k1r5 k1r9 k1r10) /* creating the fractions of self-employed workers in the household */ /* total self-employed in household */ egen nwork = rsum(mp1014 mp1519 mp2024 mp2529 mp3034 mp3539 mp4044 mp4549 mp5054 mp5559 mp6099) /* fraction of household self-employed in each age group */ replace mp1014 = mp1014/nwork replace mp1519 = mp1519/nwork replace mp2024 = mp2024/nwork replace mp2529 = mp2529/nwork replace mp3034 = mp3034/nwork replace mp3539 = mp3539/nwork replace mp4044 = mp4044/nwork replace mp4549 = mp4549/nwork replace mp5054 = mp5054/nwork replace mp5559 = mp5559/nwork replace mp6099 = mp6099/nwork drop p1014 p1519 p2024 p2529 p3034 p3539 p4044 p4549 p5054 p5559 p6099 /* regressing, without an intercept, labor component of household self-employment income on the number of self-employed workers in the household and the fraction of self-employed in each age group by household */ reg nl_income23 nwork mp1014 mp1519 mp2024 mp2529 mp3034 mp3539 mp4044 mp4549 mp5054 mp5559 mp6099 if rhead == 1 [ w = mwert ], robust noconstant /* saving coefficients and assigning them to household members based on age */ gen bp1014 = _b[mp1014] if age >= 10 & age <= 14 gen bp1519 = _b[mp1519] if age >= 15 & age <= 19 gen bp2024 = _b[mp2024] if age >= 20 & age <= 24 gen bp2529 = _b[mp2529] if age >= 25 & age <= 29 gen bp3034 = _b[mp3034] if age >= 30 & age <= 34 gen bp3539 = _b[mp3539] if age >= 35 & age <= 39 gen bp4044 = _b[mp4044] if age >= 40 & age <= 44 gen bp4549 = _b[mp4549] if age >= 45 & age <= 49 gen bp5054 = _b[mp5054] if age >= 50 & age <= 54 gen bp5559 = _b[mp5559] if age >= 55 & age <= 59 gen bp6099 = _b[mp6099] if age >= 60 & age <= 99 egen t_share = rsum(bp1014 bp1519 bp2024 bp2529 bp3034 bp3539 bp4044 bp4549 bp5054 bp5559 bp6099) /* allocating labor component of household non-labor income to those self-employed */ gen sbp1014 = work * nl_income * bp1014/t_share if age >= 10 & age <= 14 gen sbp1519 = work * nl_income * bp1519/t_share if age >= 15 & age <= 19 gen sbp2024 = work * nl_income * bp2024/t_share if age >= 20 & age <= 24 gen sbp2529 = work * nl_income * bp2529/t_share if age >= 25 & age <= 29 gen sbp3034 = work * nl_income * bp3034/t_share if age >= 30 & age <= 34 gen sbp3539 = work * nl_income * bp3539/t_share if age >= 35 & age <= 39 gen sbp4044 = work * nl_income * bp4044/t_share if age >= 40 & age <= 44 gen sbp4549 = work * nl_income * bp4549/t_share if age >= 45 & age <= 49 gen sbp5054 = work * nl_income * bp5054/t_share if age >= 50 & age <= 54 gen sbp5559 = work * nl_income * bp5559/t_share if age >= 55 & age <= 59 gen sbp6099 = work * nl_income * bp6099/t_share if age >= 60 & age <= 99 /* estimated non-labor (entrepreneurial) income */ drop bp* recode sbp* ( . = 0 ) egen nl_income_e = rsum(sbp1014 sbp1519 sbp2024 sbp2529 sbp3034 sbp3539 sbp4044 sbp4549 sbp5054 sbp5559 sbp6099) label var nl_income_e "estimated non-labor income" drop sbp* /* employee monthly earnings */ gen m_wage = wk_wage * 4 label var m_wage "monthly wage estimated from sakernas" /* individual labor income */ egen laborincome = rsum(nl_income_e m_wage) /* age profiles of labor income */ table age [ w = mwert ], content (mean laborincome)