plot_volcano {ideal} | R Documentation |
Volcano plot for log fold changes and log p-values in the ggplot2 framework, with additional support to annotate genes if provided.
plot_volcano( res_obj, FDR = 0.05, ylim_up = NULL, vlines = NULL, title = NULL, intgenes = NULL, intgenes_color = "steelblue", labels_intgenes = TRUE )
res_obj |
A |
FDR |
Numeric value, the significance level for thresholding adjusted p-values |
ylim_up |
Numeric value, Y axis upper limits to restrict the view |
vlines |
The x coordinate (in absolute value) where to draw vertical lines, optional |
title |
A title for the plot, optional |
intgenes |
Vector of genes of interest. Gene symbols if a |
intgenes_color |
The color to use to mark the genes on the main plot. |
labels_intgenes |
Logical, whether to add the gene identifiers/names close to the marked plots |
The genes of interest are to be provided as gene symbols if a symbol
column is provided in res_obj
, or else b< using the identifiers specified
in the row names
An object created by ggplot
library(airway) data(airway) airway dds_airway <- DESeq2::DESeqDataSetFromMatrix(assay(airway), colData = colData(airway), design=~cell+dex) # subsetting for quicker run, ignore the next two commands if regularly using the function gene_subset <- c( "ENSG00000103196", # CRISPLD2 "ENSG00000120129", # DUSP1 "ENSG00000163884", # KLF15 "ENSG00000179094", # PER1 rownames(dds_airway)[rep(c(rep(FALSE,99), TRUE), length.out=nrow(dds_airway))]) # 1% of ids dds_airway <- dds_airway[gene_subset,] dds_airway <- DESeq2::DESeq(dds_airway) res_airway <- DESeq2::results(dds_airway) plot_volcano(res_airway)