diffexpnb: Function for differential expression analysis
In FateID: Quantification of Fate Bias in Multipotent Progenitors

diffexpnb

R Documentation

Function for differential expression analysis

Description

This function performs differential expression analysis between two sets of single cell transcriptomes. The inference is based on a noise model or relies on the DESeq2 approach.

Usage

diffexpnb(
  x,
  A,
  B,
  DESeq = FALSE,
  method = "pooled",
  norm = FALSE,
  vfit = NULL,
  locreg = FALSE,
  ...
)

Arguments

`x`	expression data frame with genes as rows and cells as columns. Gene IDs should be given as row names and cell IDs should be given as column names. This can be a reduced expression table only including the features (genes) to be used in the analysis. This input has to be provided if `g` (see below) is given and corresponds to a valid gene ID, i. e. one of the rownames of `x`. The default value is `NULL`. In this case, cluster identities are highlighted in the plot.
`A`	vector of cell IDs corresponding column names of `x`. Differential expression in set `A` versus set `B` will be evaluated.
`B`	vector of cell IDs corresponding column names of `x`. Differential expression in set `A` versus set `B` will be evaluated.
`DESeq`	logical value. If `TRUE`, then DESeq2 is used for the inference of differentially expressed genes. In this case, it is recommended to provide non-normalized input data `x`. Default value is `FALSE`
`method`	either "per-condition" or "pooled". If DESeq is not used, this parameter determines, if the noise model is fitted for each set separately ("per-condition") or for the pooled set comprising all cells in `A` and `B`. Default value is "pooled".
`norm`	logical value. If `TRUE` then the total transcript count in each cell is normalized to the minimum number of transcripts across all cells in set `A` and `B`. Default value is `FALSE`.
`vfit`	function describing the background noise model. Inference of differentially expressed genes can be performed with a user-specified noise model describing the expression variance as a function of the mean expression. Default value is `NULL`.
`locreg`	logical value. If `FALSE` then regression of a second order polynomial is perfomed to determine the relation of variance and mean. If `TRUE` a local regression is performed instead. Default value is `FALSE`.
`...`	additional arguments to be passed to the low level function `DESeqDataSetFromMatrix`.

Value

If DESeq equals TRUE, the function returns the output of DESeq2. In this case list of the following two components is returned:

`cds`	object returned by the DESeq2 function `DESeqDataSetFromMatrix`.
`res`	data frame containing the results of the DESeq2 analysis.

Otherwise, a list of three components is returned:

`vf1`	a data frame of three columns, indicating the mean `m`, the variance `v` and the fitted variance `vm` for set `A`.
`vf2`	a data frame of three columns, indicating the mean `m`, the variance `v` and the fitted variance `vm` for set `B`.
`res`	a data frame with the results of the differential gene expression analysis with the structure of the `DESeq` output, displaying mean expression of the two sets, fold change and log2 fold change between the two sets, the p-value for differential expression (`pval`) and the Benjamini-Hochberg corrected false discovery rate (`padj`).

Examples


x <- intestine$x
y <- intestine$y
v <- intestine$v

tar <- c(6,9,13)
fb <- fateBias(x,y,tar,z=NULL,minnr=5,minnrh=10,nbfactor=5,use.dist=FALSE,seed=NULL,nbtree=NULL)

thr <- .3

A <- rownames(fb$probs)[fb$probs[,"t6"]  > .3]
B <- rownames(fb$probs)[fb$probs[,"t13"] > .3]
de <- diffexpnb(v,A=A,B=B)

FateID documentation built on June 14, 2022, 5:06 p.m.