No description, website, or topics provided.
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Failed to load latest commit information.
R
data Commit made by the Bioconductor Git-SVN bridge. Oct 13, 2014
inst
man
src Fixed the issue #2 Mar 4, 2016
vignettes Merged github and bioconductor versions Apr 29, 2018
visual_test setting up git-svn bridge Sep 30, 2014
.Rbuildignore
.gitignore
DESCRIPTION
GOsummaries.Rproj
NAMESPACE Fixed #4 #5 Mar 27, 2018
NEWS Commit made by the Bioconductor Git-SVN bridge. Oct 13, 2014
README.md

README.md

GOsummaries

An R package that visualizes the GO enrichment results as word clouds and arranges them together with figures of experimental data. This allows us to draw informative summary plots for analyses such as differential expression or clustering, where for each gene list we display its behaviour in the experiment alongside with its GO annotations. The approach is especially interesting for Principal Component Analysis (PCA), where we can annotate the principal axes functionally based on the weights of the genes.

Citation

Raivo Kolde, Jaak Vilo. 2015. “GOsummaries: an R Package for Visual Functional Annotation of Experimental Data.” F1000Research 4.

Installation

The package is available ar Bioconductor, so the installation can be done as usual.

source("https://bioconductor.org/biocLite.R")
biocLite("GOsummaries")

More comprehensive user guide can be found in the vignette.

Examples

PCA

# Perform PCA on samples
pcr = prcomp(t(tissue_example$exp))

# Create gosummaries object
gs_pca = gosummaries(pcr, annotation = tissue_example$annot)

# Plot
plot(gs_pca, classes = "Tissue")

pca

K-means clustering

library(GOsummaries)

data(tissue_example)

# Filter genes and perform k-means
sd = apply(tissue_example$exp, 1, sd)
exp2 = tissue_example$exp[sd > 0.75,]
exp2 = exp2 - apply(exp2, 1, mean)
kmr = kmeans(exp2, centers = 6, iter.max = 100)

# Create gosummaries object
gs_kmeans = gosummaries(kmr, exp = exp2, annotation = tissue_example$annot)
plot(gs_kmeans, components = 1:2)

kmeans