seqgendiff-package: seqgendiff: Sequence Generation/Modification for Differential...

Description seqgendiff Functions Author(s) References

Description

This package is designed to take real RNA-seq data and alter it by adding a known amount of signal. You can then use this modified dataset in simulation studies for differential expression analysis, factor analysis, confounder adjustment, or library size adjustment. The advantage of this way of simulating data is that you can see how your method behaves when the simulated data exhibit common (and annoying) features of real data. For example, in the real world data are not normally (or negative binomially) distributed and unobserved confounding is a major issue. This package will simulate data that exhibit these characteristics. The methods used in this package are described in detail in Gerard (2019).

seqgendiff Functions

select_counts

Subsample the columns and rows of a real RNA-seq count matrix. You would then feed this sub-matrix into one of the thinning functions below.

thin_diff

The function most users should be using for general-purpose binomial thinning. For the special applications of the two-group model or library/gene thinning, see the functions listed below.

thin_2group

The specific application of thinning in the two-group model.

thin_lib

The specific application of library size thinning.

thin_gene

The specific application of total gene expression thinning.

thin_all

The specific application of thinning all counts.

effective_cor

Returns an estimate of the actual correlation between surrogate variables and a user-specified design matrix.

ThinDataToSummarizedExperiment

Converts a ThinData object to a SummarizedExperiment object.

ThinDataToDESeqDataSet

Converts a ThinData object to a DESeqDataSet object.

Author(s)

David Gerard

References


seqgendiff documentation built on July 8, 2020, 5:55 p.m.