pub fn chi_squared_test(observed: &[f64], expected: &[f64]) -> f64
Pearson’s χ² goodness-of-fit statistic.
Returns Σ (O_i − E_i)² / E_i.
Panics if the slices have different lengths. Bins with expected count zero are skipped.