Skip to contents

This script is meant to be used prior to gl.nhybrids to maximise the information content of the SNPs used to identify hybrids (currently newhybrids does allow only 200 SNPs). The idea is to use first all loci that have fixed alleles between the potential source populations and then 'fill up' to 200 loci using loci that have private alleles between those. The functions filters for those loci (if invers is set to TRUE, the opposite is returned (all loci that are not fixed and have no private alleles - not sure why yet, but maybe useful.)

Usage, pop1, pop2, invers = FALSE, verbose = NULL)



Name of the genlight object containing the SNP data [required].


Name of the first parental population (in quotes) [required].


Name of the second parental population (in quotes) [required].


Switch to filter for all loci that have no private alleles and are not fixed [FALSE].


Verbosity: 0, silent or fatal errors; 1, begin and end; 2, progress log; 3, progress and results summary; 5, full report [default 2, unless specified using gl.set.verbosity].


The reduced genlight dataset, containing now only fixed and private alleles.


Authors: Bernd Gruber & Ella Kelly (University of Melbourne); Custodian: Luis Mijangos -- Post to


result <-, pop1=pop([1], 
#> Starting 
#>   Processing genlight object with SNP data
#> Completed: