Tweet
Hello all,
I`m using DEXSeq to find differential splicing. It operates with raw counts of reads aligned to the exonic segments. But frequently it is worth to have inclusion levels ("PSI") assigned to exons, which is proportion of transcripts including given exonic segment among all transcripts of given gene.
Calculation of absolute expression level of exon is seems to be not very difficult (something like FPKM): exon count * size factor (or divided to the size of library) / length of exon.
But for inclusion level we need to find ratio of expression level of exon (FPKM/counts normalized to length) to the expression level of whole gene:
something like "PSI" = "exonFpkm"/"geneFpkm".
Is there any smart way to calculate something like PSI from exon count objects used in DEXSeq?
Thank you.
I`m using DEXSeq to find differential splicing. It operates with raw counts of reads aligned to the exonic segments. But frequently it is worth to have inclusion levels ("PSI") assigned to exons, which is proportion of transcripts including given exonic segment among all transcripts of given gene.
Calculation of absolute expression level of exon is seems to be not very difficult (something like FPKM): exon count * size factor (or divided to the size of library) / length of exon.
But for inclusion level we need to find ratio of expression level of exon (FPKM/counts normalized to length) to the expression level of whole gene:
something like "PSI" = "exonFpkm"/"geneFpkm".
Is there any smart way to calculate something like PSI from exon count objects used in DEXSeq?
Thank you.
Comment