pub fn hannan_rissanen(x: &[f64], config: HrConfig) -> Result<HrResult>
Hannan–Rissanen stages 1–2: select (p, q) and return the (consistent, not efficient) stage-2 coefficient estimates.
(p, q)