code: def count_hom_ref(self, axis=None):
          
          b = self.is_hom_ref()
          return np.sum(b, axis=axis)
docs: Count homozygous reference genotypes.
  
          Parameters
          ----------
          axis : int, optional
              Axis over which to count, or None to perform overall count.