The locus metadata supplied by DArT has FreqHomSnp included, but the frequency of the homozygous alternate will change when some individuals are removed from the dataset.
utils.recalc.freqhomsnp(x, verbose = NULL)
The modified genlight object.
Name of the genlight object containing the SNP data [required].
Verbosity: 0, silent or fatal errors; 1, begin and end; 2, progress log; 3, progress and results summary; 5, full report [default 2].
Custodian: Luis Mijangos (Post to https://groups.google.com/d/forum/dartr)
This script recalculates the FreqHomSnp and places these recalculated values in the appropriate place in the genlight object.
Note that the frequency of the homozygote alternate SNPS is calculated from the individuals that could be scored.
utils.recalc.metrics
for recalculating all metrics,
utils.recalc.callrate
for recalculating CallRate,
utils.recalc.freqhomref
for recalculating frequency of homozygous
reference, utils.recalc.avgpic
for recalculating AvgPIC,
utils.recalc.freqhet
for recalculating frequency of heterozygotes,
gl.recalc.maf
for recalculating minor allele frequency,
gl.recalc.rdepth
for recalculating average read depth
#out <- utils.recalc.freqhomsnp(testset.gl)
Run the code above in your browser using DataLab