This package and included functions were developed for the American chestnut landscape genomics manuscript. The aim was to develop a method for estimating the number of individuals to sample from a population in order to capture a predefined percentage of diversity within a breeding population.
A VCF file with GT information or a genotype matrix can be used as input
# castgen is available from CRAN and can be installed using the following command:
install.packages("castgen")
The output will provide the following information:
Number of individuals to sample = 22.0
95% Confidence Intervals = (16.378127142052662, 27.621872857947338)
Iterations performed = 5
If you use this package, please cite:
A.M. Sandercock, J.W. Westbrook, Q. Zhang, & J.A. Holliday, A genome-guided strategy for climate resilience in American chestnut restoration populations, Proc. Natl. Acad. Sci. U.S.A. 121 (30) e2403505121, https://doi.org/10.1073/pnas.2403505121 (2024).
The capture_diversity functions were converted from a Python function https://github.com/alex-sandercock/Capture_genomic_diversity. These R versions of the function are much (much) faster to run and have the option to run in parallel depending on cores available.
