transform_RNAseq: Transform RNA-seq data using common transformations

Description Usage Arguments Value Examples

Description

Application of common transformations for RNA-seq data prior to fitting a normal mixture model

Usage

1
2
transform_RNAseq(y, norm = "TMM", transformation = "arcsin",
  geneLength = NA, meanFilterCutoff = NULL, verbose = TRUE)

Arguments

y

(n x q) matrix of observed counts for n observations and q variables

norm

The type of estimator to be used to normalize for differences in library size: “TC” for total count, “DESeq” for the normalization method in the DESeq package, and “TMM” for the TMM normalization method (Robinson and Oshlack, 2010). Can also be a vector (of length q) containing pre-estimated library size estimates for each sample.

transformation

Transformation type to be used: “arcsin”, “logit”, “logMedianRef”, “profile”, “voom”, “logRPKM” (if geneLength is provided by user), “none

geneLength

Vector of length equal to the number of rows in “y” providing the gene length (bp) for RPKM calculation

meanFilterCutoff

Value used to filter low mean normalized counts

verbose

If TRUE, include verbose output

Value

tcounts

Transformed counts

normCounts

Normalized counts

snorm

Per-sample normalization factors divided by mean normalization factor

ellnorm

Per-sample normalization factors

Examples

1
2
3
4
5
6
7
8
9
set.seed(12345)
countmat <- matrix(runif(300*4, min=0, max=500), nrow=300, ncol=4)
countmat <- countmat[which(rowSums(countmat) > 0),]
conds <- rep(c("A","B","C","D"), each=2)

## Arcsin transformation, TMM normalization
arcsin <- transform_RNAseq(countmat, norm="TMM", transformation="arcsin")$tcounts
## Logit transformation, TMM normalization
logit <- transform_RNAseq(countmat, norm="TMM", transformation="logit")$tcounts

coseq documentation built on May 2, 2019, 4:55 p.m.