escoEstDE: Estimate esco expression DE parameters for discrete cell...

Description Usage Arguments Details Value

View source: R/esco-estimate.R

Description

Parameters are estimated by comparing means of individual genes to the median mean expression level.

Usage

1
escoEstDE(nonDE.normcounts, DE.normcounts, degenes, cellinfo, params)

Arguments

nonDE.normcounts

library size normalised counts matrix for nonDE genes.

DE.normcounts

library size normalised counts matrix for DE genes.

degenes

a dataframe of two columns: "genes" and "clust", where "genes" contains the position of a degenes, and "clust" contains the corresponding cell group the degenes marks.

cellinfo

a vector of length n, where n is the number of cells. Each entries is the group identity of a cell.

params

escoParams object to store estimated values in.

Details

Expression DE factor are detected using the Median Absolute Deviation (MAD) from median method. If the log2 mean expression of a gene is greater than two MADs above the median log2 mean expression it is designated as an outlier. The proportion of outlier genes is used to estimate the outlier probability. Factors for each outlier gene are calculated by dividing mean expression by the median mean expression. A log-normal distribution is then fitted to these factors in order to estimate the outlier factor location and scale parameters using fitdist.

Value

escoParams object with estimated values.


JINJINT/ESCO documentation built on May 13, 2021, 7:25 p.m.