get_phhs_single function

Get Pseudo-Households (PHH) for a single region