Filter loci that contain private (and fixed alleles) between two populations.
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.)
gl.filter.pa(x, pop1, pop2, invers = FALSE, verbose = NULL)
x |
– name of the genlight object containing the SNP data [required] |
pop1 |
– name of the first parental population (in quotes) [required] |
pop2 |
– name of the second parental population (in quotes) [required] |
invers |
– switch to filter for all loci that have no private alleles and are not fixed [FALSE] |
verbose |
– 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
Bernd Gruber & Ella Kelly (University of Melbourne) (Post to https://groups.google.com/d/forum/dartr)
result <- gl.filter.pa(testset.gl, pop1=pop(testset.gl)[1], pop2=pop(testset.gl)[2],verbose=3)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.