data(splicingFactors)
fs <- firmaAnalysis(useToyData=TRUE)
gs <- readGeneSummaries(useToyData=TRUE)
tra <- aberrantExonUsage(1.0, fs)
# The aberrantExonUsage function also creates the 'quantiles' object with
# upper and lower threshold values for accepting aberrant exon usage, and
# the list object 'aberrantExons' with the sample-wise number of exons
# outside the threshold values.
aberrantExonsPerms <- probesetPermutations(fs, quantiles)
# Create scatter plot with the samples
scatterPlot("scatter.png", TRUE, aberrantExons, aberrantExonsPerms)
Run the code above in your browser using DataLab