sig
  type t = { nbinomial_failures : int; nbinomial_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 : failures:int -> p:float -> Distributions.NegativeBinomial.t
  val mme : int array -> Distributions.NegativeBinomial.t
end