Derivation
variable: informant (table[s]: all_marital_reports [derived table])
variable: spousenum (table[s]: all_marital_reports [derived table])
replace spousenum=. if informant==1 //only want husband number of spouses to know if union is polygamous
bys wife husband: egen heverpolyn=max(spousenum)