Hello,
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 exon usage for a given exon of the particular gene across samples.
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, just like I would plot expression across samples to visualize a couple of top DEG hits. How do I compute this relative usage matrix across all my samples to construct this heatmap? 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!
-
The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...11-06-2024, 07:24 PM -
Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...10-18-2024, 07:11 AM