I am using voom/limma to detect differential splicing via diffSplice, with visualization via plotSplice. I am trying to figure out how to make a heatmap of relative exon usage for my top hits, analagous to a heatmap of expression which illustrates the top hits of ordinary differential expression. This heatmap should clearly show the change in relative logFC for a given exon of the particular gene which is called differentially spliced.
So, for example suppose I have (for some expression cutoff
is_expr) an exon count matrix
ecounts over some samples, then:
dge_exon <- DGEList( counts = ecounts, genes = e_annotation)
dge_exon <- dge_exon[is_expr,]
dge_exon <- calcNormFactors(dge_exon)
v_exon <- voom( dge_exon, design)
fx <- eBayes(lmFit(v_exon, design))
ex <- diffSplice(fx)
and I want to plot relative exon usage as a heatmap across all my samples using the top few hits in
ex. How do I compute this relative usage matrix across all my samples? Would it be something like
eusage <- aggregate( v_exon$E, by = list(as.matrix(v_exon$genes$GeneID)), function(x) x/sum(x) )
Thank you in advance for your help!