SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
RIN values for RNAseq pogaora RNA Sequencing 14 07-15-2014 10:18 AM
RNASeq Differential Gene Expression with two-sample Mann-Whitney test on FPKM? ndovu9 Bioinformatics 4 01-27-2014 03:58 AM
How is calculated the gene expression values of microarrays from TCGA? dsmarcoantonio Bioinformatics 0 12-18-2012 04:24 AM
[NGS - analysis of gene expression data] Machine Learning + RNAseq data Chuckytah Bioinformatics 7 03-05-2012 04:16 AM

Reply
 
Thread Tools
Old 01-29-2014, 06:33 PM   #1
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default How to present a subset of gene expression values from RNAseq in a publication?

Hi!
Im curious, my supervisor want to include 5 genes from the Cuffdiff results in a table for publication.

He just wants fold changes and standard deviation (and star if significant etc), but first of all I don't know how to calculate SD for FC.. Any suggestions?

If FC+SD is not the solution, what should I send him?


Thank you very much.
sindrle is offline   Reply With Quote
Old 01-29-2014, 06:44 PM   #2
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

My suggestion so far is:
1. Take FPKM values and test values

Here is an example:

q1_FPKM q1_conf_lo q1_conf_hi
1,20921 0,404779 2,01361

q2_FPKM q2_conf_lo q2_conf_hi
0,833905 0,24662 1,4278

Gives

FC loFC hiFC
0,689627939 0,609270738 0,709074746

I present this as: 0,69 +/-0.05*

*just marks if significant.. (based on FDR)
sindrle is offline   Reply With Quote
Old 01-30-2014, 02:24 AM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Note that the confidence intervals aren't symemetric, so "0.69, 95% CI [0.61, 0.71]" would seem to make more sense.

Coming up with reasonable confidence intervals or SDs by hand when dispersion is being squeezed with an empirical Bayes method is going to be non-trivial. I remember you mentioning edgeR and such before, so you might ask on the bioconductor email list if there's a convenient function to do this built-in. The authors are pretty responsive on that list.
dpryan is offline   Reply With Quote
Old 01-30-2014, 03:14 AM   #4
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

Thank you, its my concerns as well.. But Im on a time limit and must present something very soon.

Ill try the mailing list!


Problem with edgeR is the calculation of FPKM and my problem reproducing fold changes etc from CPM/logCPM.

With Cuffdiff it is possible, FPKM calculations may be reproduced and it has a lot more easy accessible results in the output files.

Last edited by sindrle; 01-30-2014 at 03:22 AM.
sindrle is offline   Reply With Quote
Old 01-30-2014, 03:25 AM   #5
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Ah, yeah, the fold change is also shrunken (I think, I'd have to double check), so again you won't normally get the same answer by hand. As an aside, your supervisor is being a bit unreasonable to request these metrics when you actually require a pretty high level of statistical understanding to know how they were created (and your supervisor lacks that understanding). Best of luck dealing with that!
dpryan is offline   Reply With Quote
Old 01-30-2014, 03:29 AM   #6
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

I know :S

But, that is why Im asking here, to get your expert advise how to best present a subfraction of Cuffdiff results in a publication, without doing lots of stupid **** (like Im doing now) :P
sindrle is offline   Reply With Quote
Reply

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 11:00 PM.


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