sig
  type t = { bernoulli_p : float; }
  type elt = int
  val sample : ?rng:Gsl.Rng.t -> size:int -> t -> elt array
  val cumulative_probability : t -> n:elt -> float
  val probability : t -> n:elt -> float
  val mean : t -> float
  val variance : t -> float
  val skewness : t -> float
  val kurtosis : t -> float
  val create : p:float -> Distributions.Bernoulli.t
  val mle : int array -> Distributions.Bernoulli.t
end