// Creates w_nonepar_dv - Number of lone parents in hh
// INPUT FILE - w_indall.dta 
// OUTPUT FILE - w_nonepar_dv.dta

// WRITE FILEPATH WHERE YOU HAVE DOWNLOADED THE DATA, e.g. "D:/stata/stata13_se"
global in ""
// WRITE FILEPATH WHERE YOU WOULD LIKE TO SAVE THE OUTPUT FILE, e.g. "D:/nonepar_dv"
global outpath ""


// 1 is for BHPS, 2 is for UKHLS
global cohort "2" 
// INCLUDE WAVE LETTERS FOR THE WAVES NEEDED
global waves a b c d e f g h i j k l

 
**** Don’t make changes after this point ****
foreach w of global waves {
global wave "`w'_"
global num=strpos("abcdefghijklmnopqrstuvwxyz","`w'")
global dta ""
global waveno=$num
global inpath "${in}/bhps_w${num}"
if $cohort==2 {
	global waveno=$num+18	
	global inpath "${in}/ukhls"
}


use "${inpath}/${wave}indall${dta}", clear



/*-------------ORIGINAL CODE--------------*/

gen lonepar_=(${wave}single_dv==1 & (${wave}respf16_dv==1 | ${wave}respm16_dv==1))
bys ${wave}hidp: egen ${wave}nonepar_dv=sum(lonepar_)
tab ${wave}nonepar_dv
lab var ${wave}nonepar_dv "Number of lone parents in hh"
bysort ${wave}hidp: keep if _n==1
keep ${wave}hidp ${wave}nonepar_dv
save "${outpath}/${wave}nonepar_dv", replace
}