joinby ${wave}hidp ${wave}pno using ${inpath}\${wave}indall${dta}, unmatched(using)
if ${cohort}== 1 {
	cap drop ${wave}nchild_dv
	cap gen ${wave}aivfio=.
	cap gen ${wave}ivfio=.	
	}
recode ${wave}relationship_dv (.=0)
bys ${wave}hidp ${wave}pno:  ///
    egen ${wave}nchild_dv = sum((${wave}relationship_dv==9 | ${wave}relationship_dv==10 | ${wave}relationship_dv==12) ///
                            & ((${wave}aage_dv<16) | (${wave}aage_dv<0 & (${wave}ivfio>=21 & ${wave}ivfio<=24)))) 
lab var ${wave}nchild_dv "Number of own children in household"
bysort ${wave}hidp ${wave}pno: keep if _n==1
keep ${wave}hidp ${wave}pno ${wave}nchild_dv
lab var ${wave}pno "person number"
order ${wave}hidp ${wave}pno ${wave}nchild_dv
sort ${wave}hidp ${wave}pno
lab dat "Number of own children, person-level DV"
save ${outpath}/${wave}nchild_dv.dta , replace