get.numalts: Helper function for getParams

In leekgroup/derfinder: Differential expression analysis of RNA-seq data at base-pair resolution

Description

Find estimated number of alternative statistics above a set percentile of the estimated null distribution

Usage

get.numalts(pctil, null.mean, null.sd, null.prop, vals, up = TRUE)

Arguments

pctil	percentile, in (0,1), of the null distribution for which the number of alternative statistics above (if pctil is greater than 0.5) or below (if pctil is less than 0.5) is desired.
null.mean	estimated mean of null distribution (usually found with locfdrFit)
null.sd	estimated standard deviation of null distribution (usually found with locfdrFit)
null.prop	estimated proportion of statistics that came from the null distribution
vals	vector of all the observed values from the mixture distribution
up	if TRUE, get the number of overexpressed statistics above the 100pctil-th percentile of the null distribution, else get the number of underexpressed statistics below the 100pctil-th percentile of the null distribution

Details

This function is for experienced users or debugging only - all other users should use getParams, which calls this function.

Value

A list with elements

num	the estimated number of alternative values above/below val (see val below)
val	the 100pctil-th percentile of the null distribution

Author(s)

Alyssa Frazee

See Also

getParams

leekgroup/derfinder documentation built on May 20, 2019, 11:30 p.m.