hclimente / spada

⚔️Find splicing-led, functional changes of the proteome.
MIT License
10 stars 1 forks source link

Switch detection procedure #5

Closed hildemann closed 2 years ago

hildemann commented 3 years ago

Hi @hclimente , working with spada i do like the results. But since you said in a previous issue that the current version of spada is a 'complete re-implementation of the software used in the paper', I wondered whether the switch detection conditions are still the same ones, as describes in the paper. For example in the paper it is stated that : 'We only kept those cases for which the tumor isoform PSI was higher than the normal isoform in the tumor sample and the normal isoform PSI in the normal sample was higher than the value for the tumor isoform' However in my analysis, spada reported the situation seen below as a switch, even though in both samples the normal isoform had a bigger psi (y axis = psi).

So are the switch detection criteria now (in the current spada version) still the same as described in the paper ? Your advice would be greatly appreciated. Thanks!

hclimente commented 3 years ago

Thanks for your feedback.

Indeed, I did not implement that check in spada. For reference, these are the criteria that a switch must meet:

  1. "Big enough" delta PSI change in both normal and tumor isoforms.
  2. Gene not-differentially expressed.
  3. Gene expression higher than the specified threshold in both cases and controls.