Code in this repository aims to expand the functionality of Preseq and make it available in the R statistical computing enviroment. There are five ways this is supposed to work:
-
The basic functionality of the
preseqprogram, initially focusing only on library complexity, is available. These functions contain the string "rfa" as part of their names. -
The mathematical routines for doing rational function approximation via continued fractions were implemented as a wrapper for existing functionality in C++. This changed since version 3.1.1, and now all functionality is in R.
-
Fitting a zero-truncated negative binomial distribution to the sample is available. These functions include the string "ztnb" as part of the names.
-
The simulation module is used to generate samples based on mixture of Poisson.
-
Extra functions are provided to estimate the number of species represented at least r times in a random sample.
See https://cran.r-project.org/package=preseqR for details.
- Improve the user interface for core functions
- Add functions to optimize the depth of single-cell whole-genome sequencing experiments and whole-exome sequencing experiments
- Add functions to predict the sample coverage, which is the probability of sampling an observed species from a population
- Add functions to predict the fraction of k-mers represented at least r times in a sequencing experiment
- Fixed a bug for removing defects
- Substituted embedded C++ code with R code
- Removed the dependencies on the software preseq
- Fixed a bug in Chao's estimator
- Fixed issues for a Solaris C++ compiler.
- We have changed the return types of many functions in the package. These functions no longer generate estimated accumulative curves. Instead, they return function types, which are estimators for the number of species represented by at least r indivdiduals in a random sample.
- We added several estimators for predicting the number of species represented by at least r individuals in a random sample
We have changed the interfaces for most of our exported functions. We add new estimators for the number of species represented by at least r individuals in a random sample.
-
We recommand everyone to install the package
preseqRfrom CRAN. It can be easily done by in an R session by typing:>install.packages("preseqR")
-
The following instructions are for installing the package from the source. Assume the source code of
preseqRhas been pulled from the git repo and it is under the current directory. Start a session in an R interpreter and type:>install.packages("polynom") >install.packages("preseqR", repos=NULL, type="source")
Note that the package
polynomis required bypreseqR.