Description Usage Arguments Value Author(s) Examples
View source: R/quickAnalysis.R
Compare open reading frames for whippet differentially spliced events
1 2 3 | leafcutterTranscriptChangeSummary(significantEvents,
combineGeneEvents = FALSE, exons, BSgenome, NMD = FALSE,
showProgressBar = TRUE, exportGTF = NULL)
|
significantEvents |
data.frame containing information from the per_intron_results.tab file output from leafcutter. |
combineGeneEvents |
combine clusters occuring in the same gene? Currently not reccomended. |
exons |
GRanges gtf annotation of exons |
BSgenome |
BSGenome object containing the genome for the species analysed |
NMD |
Use NMD predictions? (Note: notNMD must be installed to use this feature) |
showProgressBar |
show a progress bar of alternative isoform generation? |
exportGTF |
file name to export alternative isoform GTFs (default=NULL) |
data.frame containing signficant whippet diff data and ORF change summaries
Beth Signal
1 2 3 4 5 6 7 8 9 10 | leafcutterFiles <- list.files(system.file("extdata","leafcutter/",
package = "GeneStructureTools"), full.names = TRUE)
leafcutterIntrons <- read.delim(leafcutterFiles[
grep("intron_results", leafcutterFiles)],stringsAsFactors=FALSE)
gtf <- rtracklayer::import(system.file("extdata","example_gtf.gtf",
package = "GeneStructureTools"))
exons <- gtf[gtf$type=="exon"]
g <- BSgenome.Mmusculus.UCSC.mm10::BSgenome.Mmusculus.UCSC.mm10
leafcutterTranscriptChangeSummary(significantEvents = leafcutterIntrons,
exons=exons,BSgenome = g,NMD=FALSE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.