Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DEXSeq error in estimateDispersions: match.arg(start.method, c("log(y)", "mean")) fpadilla Bioinformatics 14 07-03-2013 03:11 PM
DEXSeq gene level counts Julien Roux Bioinformatics 3 11-28-2012 01:31 AM
Using Cufflinks to calculate cassette exon inclusion rate chaodai Bioinformatics 1 07-30-2012 05:40 AM
Getting Gene-Level "Expression" Value srp33 RNA Sequencing 1 10-08-2010 02:58 PM
"Systems biology and administration" & "Genome generation: no engineering allowed" seb567 Bioinformatics 0 05-25-2010 01:19 PM

Thread Tools
Old 02-12-2013, 02:42 AM   #1
Junior Member
Location: Russia

Join Date: Feb 2012
Posts: 7
Default [DEXSeq] exon counts to "PSI" (exon inclusion level)

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.
yerbol is offline   Reply With Quote
Old 10-08-2015, 11:18 PM   #2
Location: Singapore

Join Date: Oct 2009
Posts: 44

Hi yerbol,

we have done just that:

There is also a script available that will calculate the PSI scores as explained in the protocol for you. You can either use count files or read files as input.

All the best,
DerSeb is offline   Reply With Quote
Old 11-23-2015, 07:00 AM   #3
Location: Goettingen, Germany

Join Date: Aug 2011
Posts: 28

Hi Sebastian, any chance you could provide fulltext to the article you linked? Can't get it though my university....


Neuromancer is offline   Reply With Quote
Old 11-23-2015, 06:32 PM   #4
Location: Singapore

Join Date: Oct 2009
Posts: 44

Dear Roman,

feel free to contact me at:

sebastian.schaefer and at the domain

Thanks for your interest in our work,
DerSeb is offline   Reply With Quote

dexseq, exon counts, inclusion level, psi

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:17 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO