fn f_to_p(f: f64, df1: usize, df2: usize) -> f64
Approximate p-value from F-statistic using the F-distribution.
Uses the regularized incomplete beta function approximation. For F(q, n-p) distribution.