Package: updog 2.1.7

David Gerard

updog: Flexible Genotyping for Polyploids

Implements empirical Bayes approaches to genotype polyploids from next generation sequencing data while accounting for allele bias, overdispersion, and sequencing error. The main functions are flexdog() and multidog(), which allow the specification of many different genotype distributions. Also provided are functions to simulate genotypes, rgeno(), and read-counts, rflexdog(), as well as functions to calculate oracle genotyping error rates, oracle_mis(), and correlation with the true genotypes, oracle_cor(). These latter two functions are useful for read depth calculations. Run browseVignettes(package = "updog") in R for example usage. See Gerard et al. (2018) <doi:10.1534/genetics.118.301468> and Gerard and Ferrao (2020) <doi:10.1093/bioinformatics/btz852> for details on the implemented methods.

Authors:David Gerard [aut, cre]

updog_2.1.7.tar.gz
updog_2.1.7.tar.gz(r-4.7-arm64)updog_2.1.7.tar.gz(r-4.7-x86_64)updog_2.1.7.tar.gz(r-4.6-arm64)updog_2.1.7.tar.gz(r-4.6-x86_64)
updog_2.1.7.tgz(r-4.6-emscripten)
manual.pdf |manual.html✨
card.svg |card.png
updog/json (API)
NEWS

# Install 'updog' in R:
install.packages('updog', repos = c('https://cran.r-universe.dev', 'https://cloud.r-project.org'))

Bug tracker:https://github.com/dcgerard/updog/issues

Pkgdown/docs site:https://dcgerard.github.io

Uses libs:
  • openblas– Optimized BLAS
  • c++– GNU Standard C++ Library v3
  • openmp– GCC OpenMP (GOMP) support library
Datasets:
  • snpdat - GBS data from Shirasawa et al
  • uitdewilligen - Subset of individuals and SNPs from Uitdewilligen et al (2013).

On CRAN:

Conda:

openblascppopenmp

6.61 score 3 packages 114 scripts 763 downloads 10 mentions 26 exports 37 dependencies

Last updated from:4aaeea5912. Checks:6 OK. Indexed: no.

TargetResultTimeFilesSyslog
linux-devel-arm64OK220
linux-devel-x86_64OK268
source / vignettesOK339
linux-release-arm64OK220
linux-release-x86_64OK213
wasm-releaseOK183

Exports:dbetabinomfilter_snpflexdogflexdog_fullformat_multidogget_q_arrayis.flexdogis.multidoglog_sum_explog_sum_exp_2multidogoracle_cororacle_cor_from_jointoracle_jointoracle_misoracle_mis_from_jointoracle_mis_vecoracle_mis_vec_from_jointoracle_plotpbetabinomplot_genoqbetabinomrbetabinomrflexdogrgenowem

Dependencies:assertthatclicodetoolscpp11digestdoFuturedoRNGfarverforeachfuturefuture.applyggplot2globalsgluegtableisobanditeratorslabelinglifecyclelistenvmagrittrparallellyplyrR6RColorBrewerRcppRcppArmadilloreshape2rlangrngtoolsS7scalesstringistringrvctrsviridisLitewithr

Example Use of Updog

David Gerard

Rendered fromsmells_like_updog.Rmdusingknitr::rmarkdownon Jun 15 2026.

Last update: 2025-09-24
Started: 2018-07-10

Genotyping Many SNPs with multidog()

David Gerard

Rendered frommultidog.Rmdusingknitr::rmarkdownon Jun 15 2026.

Last update: 2026-03-17
Started: 2020-01-28

Oracle Calculations

David Gerard

Rendered fromoracle_calculations.Rmdusingknitr::rmarkdownon Jun 15 2026.

Last update: 2025-09-24
Started: 2018-07-10

Simulate Next-Generation Sequencing Data

David Gerard

Rendered fromsimulate_ngs.Rmdusingknitr::rmarkdownon Jun 15 2026.

Last update: 2025-09-24
Started: 2018-07-10

Readme and manuals

Help Manual

Help pageTopics
The Beta-Binomial Distributionbetabinom dbetabinom pbetabinom qbetabinom rbetabinom
Filter SNPs based on the output of 'multidog()'.filter_snp
Flexible genotyping for polyploids from next-generation sequencing data.flexdog
Flexible genotyping for polyploids from next-generation sequencing data.flexdog_full
Return arrayicized elements from the output of 'multidog'.format_multidog
Return the probabilities of an offspring's genotype given its parental genotypes for all possible combinations of parental and offspring genotypes. This is for species with polysomal inheritance and bivalent, non-preferential pairing.get_q_array
Tests if an argument is a 'flexdog' object.is.flexdog
Tests if an argument is a 'multidog' object.is.multidog
Log-sum-exponential trick.log_sum_exp
Log-sum-exponential trick using just two doubles.log_sum_exp_2
Fit 'flexdog' to multiple SNPs.multidog
Calculates the correlation between the true genotype and an oracle estimator.oracle_cor
Calculate the correlation of the oracle estimator with the true genotype from the joint distribution matrix.oracle_cor_from_joint
The joint probability of the genotype and the genotype estimate of an oracle estimator.oracle_joint
Calculate oracle misclassification error rate.oracle_mis
Get the oracle misclassification error rate directly from the joint distribution of the genotype and the oracle estimator.oracle_mis_from_joint
Returns the oracle misclassification rates for each genotype.oracle_mis_vec
Get the oracle misclassification error rates (conditional on true genotype) directly from the joint distribution of the genotype and the oracle estimator.oracle_mis_vec_from_joint
Construct an oracle plot from the output of 'oracle_joint'.oracle_plot
Make a genotype plot.plot_geno
Draw a genotype plot from the output of 'flexdog'.plot.flexdog
Plot the output of 'multidog'.plot.multidog
Simulate GBS data from the 'flexdog' likelihood.rflexdog
Simulate individual genotypes from one of the supported 'flexdog' models.rgeno
GBS data from Shirasawa et al (2017)snpdat
Subset of individuals and SNPs from Uitdewilligen et al (2013).uitdewilligen
EM algorithm to fit weighted ash objective.wem