results_all: Extract and annotate all results from a DESeq analysis

View source: R/results_all.R

results_allR Documentation

Extract and annotate all results from a DESeq analysis

Description

Extract all possible contrasts and annotate result tables from a DESeq object. Currently supports simple designs with a single variable.

Usage

results_all(
  object,
  biomart,
  vs = "all",
  subset,
  relevel,
  alpha = 0.05,
  add_columns,
  trt,
  lfcShrink = "normal",
  simplify = TRUE,
  ...
)

Arguments

object

a DESeqDataSet

biomart

annotations from read_biomart with column 1 matching row names in results

vs

either compare all vs. all (default) or all vs specific treatment like control, or see note.

subset

index to subset all pairwise comparisons, try combn(sort(samples$trt),2)

relevel

Levels to compare, if missing then levels(dds$trt)

alpha

the significance cutoff for the adjusted p-value cutoff (FDR)

add_columns

a vector of biomart columns to add to result table, default gene_name, biotype, chromosome, description and human_homolog if present

trt

Compare groups within trt group, default is first term in the design formula

lfcShrink

shrink fold changes using lfcShrink, default normal, ashr or none. Use results_apeglm for apeglm shrinkage. The MLE or unshrunken fold change is saved and added to the results as MLE_log2FC

simplify

return a tibble if only 1 contrast present

...

additional options passed to results

Value

A list of tibbles for each contrast

Note

If you combine factors of interest into a single group following section 3.3 in the DESeq2 vignette, you can set vs = "combined" to limit the comparisons. See check_contrasts for details.

Author(s)

Chris Stubben

Examples

## Not run: 
library(hciRdata)
res <- results_all(pasilla$dds, fly98)
res

## End(Not run)

HuntsmanCancerInstitute/hciR documentation built on Aug. 20, 2024, 12:22 a.m.