Description Usage Arguments Details Value Author(s) References See Also Examples

Decompose the gene-specific variance into biological and technical components for single-cell RNA-seq data.

1 2 3 4 5 | ```
## S4 method for signature 'matrix,list'
decomposeVar(x, fit, design=NA, subset.row=NULL, ...)
## S4 method for signature 'SCESet,list'
decomposeVar(x, fit, subset.row=NULL, ..., assay="exprs", get.spikes=FALSE)
``` |

`x` |
A numeric matrix of normalized log-expression values, where each column corresponds to a cell and each row corresponds to an endogenous gene. Alternatively, a SCESet object containing such a matrix. |

`fit` |
A list containing the output of |

`design` |
A numeric matrix describing the uninteresting factors contributing to expression in each cell. |

`subset.row` |
A logical, integer or character scalar indicating the rows of |

`...` |
For |

`assay` |
A string specifying which assay values to use, e.g., |

`get.spikes` |
A logical scalar specifying whether decomposition should be performed for spike-ins. |

This function computes the variance of the normalized log-counts for each endogenous gene.
The technical component of the variance for each gene is determined by interpolating the fitted trend in `fit`

at the mean log-count for that gene.
This represents variance due to sequencing noise, variability in capture efficiency, etc.
The biological component is determined by subtracting the technical component from the total variance.

Highly variable genes (HVGs) can be identified as those with large biological components. Unlike other methods for decomposition, this approach estimates the variance of the log-counts rather than of the counts themselves. The log-transformation blunts the impact of large positive outliers and ensures that the HVG list is not dominated by outliers. Interpretation is not compromised – HVGs will still be so, regardless of whether counts or log-counts are considered.

The design matrix can be set if there are factors that should be blocked, e.g., batch effects, known (and uninteresting) clusters.
If `NULL`

, it will be set to an all-ones matrix, i.e., all cells are replicates.
If `NA`

, it will be extracted from `fit$design`

, assuming that the same cells were used to fit the trend.

Users can also directly specify which rows to use with `subset.row`

.
This is equivalent to running `decomposeVar`

on `x[subset.row,]`

, but is more efficient as it avoids the construction of large temporary matrices.

If `assay="exprs"`

and the size factors are not centred at unity, a warning will be raised - see `?trendVar`

for details.

A data frame is returned where each row corresponds to and is named after a row of `x`

(if `subset.row=NULL`

; otherwise, each row corresponds to an element of `subset.row`

).
This contains the numeric fields:

`mean`

:Mean normalized log-count per gene.

`total`

:Variance of the normalized log-counts per gene.

`bio`

:Biological component of the variance.

`tech`

:Technical component of the variance.

`p.value, FDR`

:Raw and adjusted p-values for the test against the null hypothesis that

`bio=0`

.

Rows corresponding to spike-in transcripts have their p-value and FDR fields set to `NA`

unless `get.spikes=TRUE`

.

Aaron Lun

Lun ATL, McCarthy DJ and Marioni JC (2016).
A step-by-step workflow for low-level analysis of single-cell RNA-seq data with Bioconductor.
*F1000Res.* 5:2122

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | ```
example(computeSpikeFactors) # Using the mocked-up data 'y' from this example.
y <- computeSumFactors(y) # Size factors for the the endogenous genes.
y <- computeSpikeFactors(y, general.use=FALSE) # Size factors for spike-ins.
y <- normalize(y) # Normalizing the counts by the size factors.
# Decomposing technical and biological noise.
fit <- trendVar(y)
results <- decomposeVar(y, fit)
head(results)
plot(results$mean, results$total)
o <- order(results$mean)
lines(results$mean[o], results$tech[o], col="red", lwd=2)
plot(results$mean, results$bio)
# A trend fitted to endogenous genes can also be used, pending assumptions.
fit.g <- trendVar(y, use.spikes=FALSE)
results.g <- decomposeVar(y, fit.g)
head(results.g)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.