pub fn bic(n: usize, rss: f64, p: usize) -> f64
Compute BIC from residual sum of squares.
BIC = n * ln(RSS/n) + ln(n) * p