calculate.residuals: Calculate residuals

In OutSeekR: Statistical Approach to Outlier Detection in RNA-Seq and Related Data

Calculate residuals

Description

Calculate residuals between quantiles of the input and quantiles of one of four distributions: normal, log-normal, exponential, or gamma.

Usage

calculate.residuals(x, distribution)

Arguments

x	A numeric vector.
distribution	A number corresponding to the optimal distribution of x as returned by, e.g., identify.bic.optimal.data.distribution(). One of 1 = normal, 2 = log-normal, 3 = exponential, and 4 = gamma.

Value

A numeric vector of the same length as x. Names are not retained.

Examples

# Generate fake data.
set.seed(1234);
x <- rgamma(
    n = 20,
    shape = 2,
    scale = 2
    );
names(x) <- paste(
    'Sample',
    seq_along(x),
    sep = '.'
    );
calculate.residuals(
    x = x,
    distribution = 4
    );

OutSeekR documentation built on April 11, 2025, 5:39 p.m.