plot_gsea: Plot enriched Gene Sets

In squirrelandr/DEP: Differential Enrichment analysis of Proteomics data

Description

plot_gsea plots enriched gene sets from Gene Set Enrichment Analysis.

Usage

plot_gsea(gsea_results, number = 10, alpha = 0.05, contrasts = NULL,
  databases = NULL, nrow = 1, term_size = 8)

Arguments

gsea_results	Data.frame, Gene Set Enrichment Analysis results object. (output from test_gsea()).
number	Numeric(1), Sets the number of enriched terms per contrast to be plotted.
alpha	Numeric(1), Sets the threshold for the adjusted P value.
contrasts	Character, Specifies the contrast(s) to plot. If 'NULL' all contrasts will be plotted.
databases	Character, Specifies the database(s) to plot. If 'NULL' all databases will be plotted.
nrow	Numeric(1), Sets the number of rows for the plot.
term_size	Numeric(1), Sets the text size of the terms.

Value

A barplot of the enriched terms (generated by ggplot).

Examples

# Load example
data <- UbiLength
data <- data[data$Reverse != "+" & data$Potential.contaminant != "+",]
data_unique <- make_unique(data, "Gene.names", "Protein.IDs", delim = ";")

# Make SummarizedExperiment
columns <- grep("LFQ.", colnames(data_unique))
exp_design <- UbiLength_ExpDesign
se <- make_se(data_unique, columns, exp_design)

# Filter, normalize and impute missing values
filt <- filter_missval(se, thr = 0)
norm <- normalize_vsn(filt)
imputed <- impute(norm, fun = "MinProb", q = 0.01)

# Test for differentially expressed proteins
diff <- diff <- test_diff(imputed, "control", "Ctrl")
dep <- add_rejections(diff, alpha = 0.05, lfc = 1)

## Not run: 

# Test enrichments
gsea_results <- test_gsea(dep)
plot_gsea(gsea_results)


## End(Not run)

squirrelandr/DEP documentation built on May 7, 2019, 9:31 a.m.