code: def getIndexOfHeteroMen(genotypes, menIndex):
      
      toRemove = set()
      for i in menIndex[0]:
          for genotype in [set(j.split(" ")) for j in genotypes[:, i]]:
              if len(genotype) != 1:
                  # We have an heterozygous
                  toRemove.add(i)
  
      toRemove = list(toRemove)
      toRemove.sort()
  
      return (np.array(toRemove, dtype=int),)
docs: Get the indexes of heterozygous men.
  
      :param genotypes: the genotypes of everybody.
      :param menIndex: the indexes of the men (for the genotypes).
  
      :type genotypes: numpy.array
      :type menIndex: numpy.array
  
      :returns: a :py:class:`numpy.array` containing the indexes of the genotypes
                to remove.
  
      Finds the mean that have a heterozygous genotype for this current marker.
      Usually used on sexual chromosomes.