sig
  type estimate = {
    point : float;
    lower_bound : float;
    upper_bound : float;
    confidence_level : float;
  }
  val bca :
    ?rng:Gsl.Rng.t ->
    ?confidence_level:float ->
    estimator:(float array -> float) ->
    n_iter:int -> float array -> Resampling.Bootstrap.estimate
end