We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent f852803 commit 24de004Copy full SHA for 24de004
src/univariate/continuous/ksdist.jl
@@ -28,7 +28,7 @@ function cdf(d::KSDist,x::Float64)
28
return 0.0
29
elseif b <= 1
30
# accuracy could be improved
31
- return exp(lfact(n)+n*(log(2*b-1)-log(n)))
+ return exp(logfactorial(n)+n*(log(2*b-1)-log(n)))
32
elseif x >= 1
33
return 1.0
34
elseif b >= n-1
@@ -56,7 +56,7 @@ function ccdf(d::KSDist,x::Float64)
56
if b <= 0.5
57
58
59
- return 1-exp(lfact(n)+n*(log(2*b-1)-log(n)))
+ return 1-exp(logfactorial(n)+n*(log(2*b-1)-log(n)))
60
61
62
0 commit comments