docs: def insertion_rate(Nref, Ninsertions, eps=numpy.spacing(1)):
      """Insertion rate
  
      Parameters
      ----------
      Nref : int >=0
          Number of entries in the reference.
  
      Ninsertions : int >=0
          Number of insertions.
  
      eps : float
          eps.
          Default value numpy.spacing(1)
  
      Returns
      -------
      insertion_rate: float
          Insertion rate
  
      """
  
      return float(Ninsertions / (Nref + eps))