Update src/Statistics.jl

Co-authored-by: Milan Bouchet-Valat <[email protected]>

Author: andreasnoack

src/Statistics.jl

@@ -1035,9 +1035,9 @@ end
     @assert n > 0 # this case should never happen here
 
     m = alpha + p * (one(alpha) - alpha - beta)
-    # Using muladd here avoids some rounding errors when aleph is an integer
+    # Using fma here avoids some rounding errors when aleph is an integer
     # The use of oftype supresses the promotion caused by alpha and beta
-    aleph = muladd(n, p, oftype(p, m))
+    aleph = fma(n, p, oftype(p, m))
     j = clamp(trunc(Int, aleph), 1, n - 1)
     γ = clamp(aleph - j, 0, 1)