Description Usage Arguments Details Value Author(s) References Examples
Applies a number of filters on a spliceR object to reduce data set size before running downstream analyses.
1 | preSpliceRFilter(spliceRobject, filters, expressionCutoff=0)
|
spliceRobject |
a |
filters |
vector, giving the filters that should be applied - any combinations of 'geneOK', 'expressedGenes', 'sigGenes', 'isoOK', 'expressedIso', 'isoClass' and/or 'sigIso'. Works only for data from cufflinks, as a manually generated |
expressionCutoff |
Numeric, giving the expression threshold (often in FPKM) used for the 'expressedGenes' and 'expressedIso' filter. Default value is 0. |
Often, many genes and isoforms are flagged as not "OK" or "LOWDATA" by Cufflinks, indicating low confidence in these. This function is handy for reducing the data size of a Cufflinks data set to reduce running times for downstream analyses.
Note, that preSpliceRFilter
removes trancsripts from the dataset permanently, reducing size, while the filter options of spliceR
and annotatePTC
only selects transcripts for analysis, but does not remove any data.
A SpliceRList
with transcripts after filtering.
Kristoffer Vitting-Seerup, Johannes Waage
Vitting-Seerup K , Porse BT, Sandelin A, Waage J. (2014) spliceR: an R package for classification of alternative splicing and prediction of coding potential from RNA-seq data. BMC Bioinformatics 15:81.
1 2 3 4 5 6 7 8 | #Load cufflinks example data
cuffDB <- prepareCuffExample()
#Generate SpliceRList from cufflinks data
cuffDB_spliceR <- prepareCuff(cuffDB)
#Filter
cuffDB_spliceR_filtered <- preSpliceRFilter(cuffDB_spliceR, filters=c("expressedIso", "isoOK", "expressedGenes", "geneOK"))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.