The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
BioUtils is an end-to-end R toolkit for analyzing gene expression data from GEO datasets. It provides a unified workflow for differential expression, statistical testing, visualization, and biological interpretation.
install.packages("remotes")
remotes::install_github("spencertreadway/BioUtils")# Load data
eset <- load.geo.soft("GDS507.soft", log.transform = TRUE)
geo <- extract.expression(eset)
# PCA visualization
pca.plot(geo$expression, geo$phenotype, color.by = "disease.state")
# Differential expression
de.results <- run.limma.de(geo)
# Volcano plot
volcano.plot(de.results, fc.threshold = 0.3)
# Select top genes
top.genes <- head(rownames(de.results[order(de.results$adj.P.Val), ]), 5)
probe.ids <- find.probe.by.gene(geo$gene, top.genes)
# Single gene analysis
expr <- get.gene.expression(geo$expression, probe.ids[1])
df <- build.analysis.df(expr, geo$phenotype, geo$gene)
gene.analysis.plot(df)
# LASSO model
phenotype.binary <- ifelse(geo$phenotype$disease.state == "disease", 1, 0)
lasso.fit <- fit.lasso(geo$expression, phenotype.binary)BioUtils integrates multiple layers of analysis:
Full documentation is available at https://spencertreadway.github.io/BioUtils/ or via:
help(package = "BioUtils")MIT
Spencer Treadway
These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.