Nothing
## ----message=FALSE------------------------------------------------------------
library(psichomics)
library(plyr)
## ---- include=FALSE-----------------------------------------------------------
suppaOutput <- system.file("extdata/eventsAnnotSample/suppa_output/suppaEvents",
package="psichomics")
suppaFile <- tempfile(fileext=".RDS")
## ---- results="hide"----------------------------------------------------------
# suppaOutput <- "path/to/SUPPA/output"
# Replace `genome` for the string with the identifier before the first
# underscore in the filenames of that directory (for instance, if one of your
# filenames of interest is "hg19_A3.ioe", the string would be "hg19")
suppa <- parseSuppaAnnotation(suppaOutput, genome="hg19")
annot <- prepareAnnotationFromEvents(suppa)
# suppaFile <- "suppa_hg19_annotation.RDS"
saveRDS(annot, file=suppaFile)
## ---- include=FALSE-----------------------------------------------------------
matsOutput <- system.file("extdata/eventsAnnotSample/mats_output/ASEvents/",
package="psichomics")
matsFile <- tempfile("mats", fileext=".RDS")
## ---- results="hide"----------------------------------------------------------
# matsOutput <- "path/to/rMATS/output"
mats <- parseMatsAnnotation(
matsOutput, # Output directory from rMATS
genome = "fromGTF", # Identifier of the filenames
novelEvents=TRUE) # Parse novel events?
annot <- prepareAnnotationFromEvents(mats)
# matsFile <- "mats_hg19_annotation.RDS"
saveRDS(annot, file=matsFile)
## ---- include=FALSE-----------------------------------------------------------
misoAnnotation <- system.file("extdata/eventsAnnotSample/miso_annotation",
package="psichomics")
misoFile <- tempfile("miso", fileext=".RDS")
## ---- results="hide"----------------------------------------------------------
# misoAnnotation <- "path/to/MISO/annotation"
miso <- parseMisoAnnotation(misoAnnotation)
annot <- prepareAnnotationFromEvents(miso)
# misoFile <- "miso_AS_annotation_hg19.RDS"
saveRDS(annot, file=misoFile)
## ---- include=FALSE-----------------------------------------------------------
vastAnnotation <- system.file("extdata/eventsAnnotSample/VASTDB/Hsa/TEMPLATES/",
package="psichomics")
vastFile <- tempfile("vast", fileext=".RDS")
## ---- results="hide"----------------------------------------------------------
# vastAnnotation <- "path/to/VASTDB/libs/TEMPLATES"
vast <- parseVastToolsAnnotation(vastAnnotation, genome="Hsa")
annot <- prepareAnnotationFromEvents(vast)
# vastFile <- "vast_AS_annotation_hg19.RDS"
saveRDS(annot, file=vastFile)
## ---- include=FALSE-----------------------------------------------------------
annotFile <- tempfile(fileext=".RDS")
## ---- results="hide"----------------------------------------------------------
# Combine the annotation from SUPPA, MISO, rMATS and VAST-TOOLS
annot <- prepareAnnotationFromEvents(suppa, vast, mats, miso)
# annotFile <- "AS_annotation_hg19.RDS"
saveRDS(annot, file=annotFile)
## ---- results="hide"----------------------------------------------------------
annot <- readRDS(annotFile) # "annotFile" is the path to the annotation file
junctionQuant <- readFile("ex_junctionQuant.RDS") # example set
psi <- quantifySplicing(annot, junctionQuant)
## -----------------------------------------------------------------------------
psi # may have 0 rows because of the small junction quantification set
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.