whoa: Evaluation of Genotyping Error in Genotype-by-Sequencing Data

This is a small, lightweight package that lets users investigate the distribution of genotypes in genotype-by-sequencing (GBS) data where they expect (by and large) Hardy-Weinberg equilibrium, in order to assess rates of genotyping errors and the dependence of those rates on read depth. It implements a Markov chain Monte Carlo (MCMC) sampler using 'Rcpp' to compute a Bayesian estimate of what we call the heterozygote miscall rate for restriction-associated digest (RAD) sequencing data and other types of reduced representation GBS data. It also provides functions to generate plots of expected and observed genotype frequencies. Some background on these topics can be found in a recent paper "Recent advances in conservation and population genomics data analysis" by Hendricks et al. (2018) <doi:10.1111/eva.12659>, and another paper describing the MCMC approach is in preparation with Gordon Luikart and Thierry Gosselin.

Version: 0.0.2
Depends: R (≥ 3.3.0)
Imports: dplyr, magrittr, tibble, tidyr, Rcpp (≥ 0.12.16), vcfR, viridis, ggplot2
LinkingTo: Rcpp
Suggests: knitr, rmarkdown
Published: 2021-08-11
Author: Eric C. Anderson [aut, cre]
Maintainer: Eric C. Anderson <eric.anderson at noaa.gov>
License: CC0
NeedsCompilation: yes
Materials: NEWS
CRAN checks: whoa results


Reference manual: whoa.pdf
Vignettes: whoa Tutorial
Package source: whoa_0.0.2.tar.gz
Windows binaries: r-devel: whoa_0.0.2.zip, r-devel-UCRT: whoa_0.0.2.zip, r-release: whoa_0.0.2.zip, r-oldrel: whoa_0.0.2.zip
macOS binaries: r-release (arm64): whoa_0.0.2.tgz, r-release (x86_64): whoa_0.0.2.tgz, r-oldrel: whoa_0.0.2.tgz
Old sources: whoa archive


Please use the canonical form https://CRAN.R-project.org/package=whoa to link to this page.