Change comparison operator to avoid problematic p=0 edge case

JuliaStats · Nov 12, 2023 · f935edc · f935edc
1 parent 7e232ca
commit f935edc
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/univariate/discrete/bernoulli.jl b/src/univariate/discrete/bernoulli.jl
@@ -114,7 +114,7 @@ cf(d::Bernoulli, t::Real) = failprob(d) + succprob(d) * cis(t)
 
 #### Sampling
 
-rand(rng::AbstractRNG, d::Bernoulli) = rand(rng) <= succprob(d)
+rand(rng::AbstractRNG, d::Bernoulli) = rand(rng) < succprob(d)
 
 #### MLE fitting