sig
  type bandwidth = Silverman | Scott
  type kernel = Gaussian
  val estimate_pdf :
    ?kernel:Sample.KDE.kernel ->
    ?bandwidth:Sample.KDE.bandwidth ->
    ?n_points:int -> float array -> float array * float array
end