fn sse_ols(x: &[f64], y: &[f64], n_pts: usize, x_cols: usize) -> f64
Sum of squared errors for OLS regression: SSE = ||y - X beta||^2.