Hello Jianguang,

The tool package CummeRbund would likely be of interest
http://compbio.mit.edu/cummeRbund/

Support can be obtained through the these mailing lists:
http://lists.sourceforge.net/lists/listinfo/bowtie-bio-announce

They would also be best folks to address your new questions about the the underlying mechanism of how splice overload is calculated.

In general, converting to fold starts by converting a data metric into a log value (log2 generally). A google can bring up many resources that describe the calculation and various tools (R, etc.). This is common for expression analysis, but whether this is recommended by the tool authors for Cuffdiff output analysis is another good question to ask them.

Good luck with your project,

Jen
Galaxy team

On 9/6/12 10:35 AM, Du, Jianguang wrote:
Hi Jen,
Thank you for your answer.
However, the output file "transcript differential expression testing" gives the 
ratio (log2 of the fold change) of FPKM of a specific transcript between two conditions, 
which means this fold change in FPKM does not take the overall gene expression into 
consideration (the expression of one gene many be much higher in condition A than 
condidtion B) and therefore can not be used as difference of alternative splicing.
What I am doing is looking for the difference of splicing between two cell types by examining the output file "splicing differential 
expression testing". In this file, column 10 gives the value of "sqrt(JS)" (the explaination of it is "The splice 
overloading of the primary transcript, as measured by the square root of the Jensen-Shannon divergence computed on the relative abundances 
of the splice variants", and the value is never larger than 1). My understanding is that, this value has already taken the overall 
gene expression into consideration. But I do not know how much "sqrt(JS)" equals to 2 fold of change because I want to focus the 
alternative splicing with 2 fold difference between two cell types. Do you know how to convert the value of "sqrt(JS)" into 
"fold"?
In addition, how to understand the sentence "The splice overloading of the primary 
transcript"? If one gene have 3 transcript: A, B and C and the expression of these 
transcripts is: A=60%, B=25%, and C=15%, do you mean the primary transcript is A? Does 
the Cuffdiff take the overall expression (A+B+C=100%) or just the primary transcript 
(A=60%) into consideration when calculates the ratio of transcript B? Actually, it would 
be much easier for us if Cuffdiff calculates the ratio of the expression of the 
alternatively spliced exon to overall gene expression, and then compare between 
conditions.
Thanks in advance,
Jianguang

________________________________________
From: Jennifer Jackson [j...@bx.psu.edu]
Sent: Thursday, September 06, 2012 12:38 PM
To: Du, Jianguang
Cc: galaxy-user@lists.bx.psu.edu; closetic...@galaxyproject.org
Subject: Re: [galaxy-user] Please help to understand the square root of 
Jensen-Shannon divergence

Hello Jianguang,

Fold is included in the Cuffdiff output. Section "Differential
expression tests", first file, column #9.
http://cufflinks.cbcb.umd.edu/manual.html

Hopefully this helps,

Jen
Galaxy team

On 9/4/12 1:16 PM, Du, Jianguang wrote:
Dear All,

I am looking for the differential splicing events between cell types.
However the Cuffdiff gives output using "the square root of
Jensen-shannon divergence" to measure the difference.

Although I tried my best to understand "the definition of the square
root of Jensen-shannon divergence", I still could not understand the
meaning of a specific value of "the square root of Jensen-shannon
divergence". I would appreciate it very much if anyone let me know how
to covert "the square root of Jensen-shannon divergence" into "fold".
For example, how much "the square root of Jensen-shannon divergence" is
"2 fold difference" equal to.

Thanks in advance.

Jianguang



___________________________________________________________
The Galaxy User list should be used for the discussion of
Galaxy analysis and other features on the public server
at usegalaxy.org.  Please keep all replies on the list by
using "reply all" in your mail client.  For discussion of
local Galaxy instances and the Galaxy source code, please
use the Galaxy Development list:

    http://lists.bx.psu.edu/listinfo/galaxy-dev

To manage your subscriptions to this and other Galaxy lists,
please use the interface at:

    http://lists.bx.psu.edu/


--
Jennifer Jackson
http://galaxyproject.org


--
Jennifer Jackson
http://galaxyproject.org
___________________________________________________________
The Galaxy User list should be used for the discussion of
Galaxy analysis and other features on the public server
at usegalaxy.org.  Please keep all replies on the list by
using "reply all" in your mail client.  For discussion of
local Galaxy instances and the Galaxy source code, please
use the Galaxy Development list:

 http://lists.bx.psu.edu/listinfo/galaxy-dev

To manage your subscriptions to this and other Galaxy lists,
please use the interface at:

 http://lists.bx.psu.edu/

Reply via email to