Package: updog 2.1.7

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.zip(r-4.7)updog_2.1.7.zip(r-4.6)updog_2.1.7.zip(r-4.5)
updog_2.1.7.tgz(r-4.6-x86_64)updog_2.1.7.tgz(r-4.6-arm64)updog_2.1.7.tgz(r-4.5-x86_64)updog_2.1.7.tgz(r-4.5-arm64)
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✨
DESCRIPTION |NEWS
card.svg |card.png
updog/json (API)

# Install 'updog' in R:
install.packages('updog', repos = c('https://dcgerard.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

9.43 score 33 stars 3 packages 114 scripts 790 downloads 10 mentions 26 exports 37 dependencies

Last updated from:04e8c4243f. Checks:13 OK. Indexed: yes.

TargetResultTimeFilesSyslog
linux-devel-arm64OK236
linux-devel-x86_64OK231
source / vignettesOK378
linux-release-arm64OK224
linux-release-x86_64OK237
macos-release-arm64OK207
macos-release-x86_64OK231
macos-oldrel-arm64OK128
macos-oldrel-x86_64OK388
windows-develOK229
windows-releaseOK206
windows-oldrelOK201
wasm-releaseOK212

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

Genotyping Many SNPs with multidog()
Abstract | Fit multidog() | multidog() Output | References
Rendered frommultidog.Rmdusingknitr::rmarkdown

Last update: 2026-02-03
Started: 2019-12-05

Example Use of Updog
What's Updog? | Example from an S1 Population | Fit updog | Analyze Output | Filtering SNPs | References
Rendered fromsmells_like_updog.Rmdusingknitr::rmarkdown

Last update: 2020-02-27
Started: 2018-05-04

Oracle Calculations
Abstract | Controlling Misclassification Error | Visualizing the Joint Distribution | References
Rendered fromoracle_calculations.Rmdusingknitr::rmarkdown

Last update: 2020-02-27
Started: 2018-05-11

Simulate Next-Generation Sequencing Data
Abstract | Analysis | F1 Population | HWE Population | References
Rendered fromsimulate_ngs.Rmdusingknitr::rmarkdown

Last update: 2020-02-27
Started: 2018-05-04

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