updog

Build Status Appveyor Build status codecov License: GPL v3 CRAN_Status_Badge

Updog provides a suite of methods for genotyping polyploids from next-generation sequencing (NGS) data. It does this while accounting for many common features of NGS data: allele bias, overdispersion, sequencing error, and (possibly) outlying observations. It is named updog for “Using Parental Data for Offspring Genotyping” because we originally developed the method for full-sib populations, but it works now for more general populations. The method is described in detail Gerard et. al. (2018) <doi:10.1534/genetics.118.301468>. Additional details concerning prior specification are described in Gerard and Ferrão (2019) <doi:10.1093/bioinformatics/btz852>.

The main function is flexdog(), which provides many options for the distribution of the genotypes in your sample. Novel genotype distributions include the class of proportional normal distributions (model = "norm") and the class of discrete unimodal distributions (model = "ash"). The default is model = "norm" because it is the most robust to varying genotype distributions, but feel free to use more specialized priors if you have more information on the data.

Also provided are:

The original updog package is now named updogAlpha and may be found here.

See also ebg, fitPoly, and TET, and polyRAD. Our best “competitor” is probably fitPoly, though polyRAD has some nice ideas for utilizing population structure and linkage disequilibrium.

See NEWS for the latest updates on the package.

Vignettes

I’ve included many vignettes in updog, which you can access online here.

Bug Reports

If you find a bug or want an enhancement, please submit an issue here.

Installation

You can install updog from CRAN in the usual way:

install.packages("updog")

You can install the current (unstable) version of updog from GitHub with:

# install.packages("devtools")
devtools::install_github("dcgerard/updog")

CVXR

If you want to use the use_cvxr = TRUE option in flexdog (not generally recommended), you will need to install the CVXR package. Before I could install CVXR in Ubuntu, I had to run in the terminal

sudo apt-get install libmpfr-dev

and then run in R

install.packages("Rmpfr")

How to Cite

Please cite

Gerard, D., Ferrão, L. F. V., Garcia, A. A. F., & Stephens, M. (2018). Genotyping Polyploids from Messy Sequencing Data. Genetics, 210(3), 789-807. doi: 10.1534/genetics.118.301468.

Or, using BibTex:

@article {gerard2018genotyping,
    author = {Gerard, David and Ferr{\~a}o, Lu{\'i}s Felipe Ventorim and Garcia, Antonio Augusto Franco and Stephens, Matthew},
    title = {Genotyping Polyploids from Messy Sequencing Data},
    volume = {210},
    number = {3},
    pages = {789--807},
    year = {2018},
    doi = {10.1534/genetics.118.301468},
    publisher = {Genetics},
    issn = {0016-6731},
    URL = {https://doi.org/10.1534/genetics.118.301468},
    journal = {Genetics}
}

If you are using the proportional normal prior class (model = "norm") or the unimodal prior class (model = "ash"), then please also cite

Gerard, D. & Ferrão L. F. V. (2019). “Priors for Genotyping Polyploids.” Bioinformatics (in press). doi: 10.1093/bioinformatics/btz852

Or, using BibTex:

@article{gerard2019priors,
    author = {Gerard, David and Ferr{\~a}o, Lu{\'i}s Felipe Ventorim},
    title = {Priors for Genotyping Polyploids},
    journal = {Bioinformatics},
    year = {2019},
    month = {11},
    issn = {1367-4803},
    doi = {10.1093/bioinformatics/btz852},
    note = {btz852},
}

Code of Conduct

Please note that this project is released with a Contributor Code of Conduct. By participating in this project you agree to abide by its terms.