#Load cufflinks example data
cuffDB <- prepareCuffExample()
#Generate SpliceRList from cufflinks data
cuffDB_spliceR <- prepareCuff(cuffDB)
#Reduce dataset size for fast example runtime
cuffDB_spliceR[[1]] <- cuffDB_spliceR[[1]][1:500]
#Run spliceR
mySpliceRList <- spliceR(cuffDB_spliceR, compareTo='preTranscript', filters=c('expressedGenes','geneOK', 'isoOK', 'expressedIso', 'isoClass'))
#Get top dIF transcripts
topIsoShift(mySpliceRList, n=20)
Run the code above in your browser using DataLab