A two part geoadditive model for geographical domain estimation