let quantile_of_freedom_degree =
      fun ?(p=`p0_95) ->
        let pi = indice_of_probability p in
        fun n ->
          if n < 1 then
            raise (Invalid_degree (Printf.sprintf "Invalid freedom degree %d" n))
          else if n > Array.length quantiles then
            raise (Invalid_degree
                     (Printf.sprintf "No predefined quantile for probability %s and freedom degree %n"
                        (string_of_probability p) n)
                  )
          else
            quantiles.(n-1).(pi)