Determines the domain mode (overall, binary, or factor) from `yz_data` and calls `dispatch_estimator_one_domain()` for each domain, returning a unified result list.
dispatch_estimator(build, yz_data)A list with components:
- `type`: either `"single"` or `"multi"`. - `labels`: character vector of domain labels. - `estimates`: a list of per-domain results, each with `mean` and `variance`.
A `pw_fit` object returned by the build step.
A list returned by `process_na_yz()`, containing `Y`, `w`, `X`, and `domain`.