R arm of morie.fn.ksamp: the Scholz-Stephens (1987) k-sample
Anderson-Darling statistic in its midrank (ties) form \(A^2_{akN}\),
standardised \(T = (A^2_{akN} - (k - 1))/\sigma_N\) with the exact
finite-sample variance, and the p-value interpolated as in scipy: a
least-squares quadratic of the log significance levels in the Scholz-
Stephens critical values, clipped to the range 0.001 to 0.25.
