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