SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Cuffdiff-confidence interval (http://seqanswers.com/forums/showthread.php?t=23109)

found 09-06-2012 06:32 AM

Cuffdiff-confidence interval
 
Hi, everybody, I have a problem about cuffdiff. I used cuffdiff to test differential expression. However, I checked the "genes.fpkm_tracking". All the FPKM estimate have q1_conf_lo=0, and q1_conf_hi a big number, which means the 95% confidence interval are very large. It seems not reasonable. I am not sure if the test statistic is calculated by this intervals. If it is, then the number of the deferentially expressed genes would be smaller.

I also try to use cufflinks to estimate the expression level. and "genes.fpkm_tracking" produced by cufflinks have reasonable intervals. I remember cuffdiff used cufflinks quantification engine. So they should have similar result, right? Anyone have any ideas for this problems?

Thanks a lot

IBseq 10-17-2012 01:02 PM

did you get any answer?

ib

IBseq 10-17-2012 01:03 PM

did you find a solution? I have the same problem...

ib

found 10-23-2012 05:31 AM

No, I did not get reply.
However, I think the variance difference is just caused by numerical reasons. Since I do not have experimental replicate, using cufflinks to estimate the expression only consider the variance from the model, such as the variation resulting from the assignment a read to which isoform or possion distribution. this variance do not have any biological meaning. The true variance should include both the model variance and the biological replicate variance.

In cuffdiff, which compare expression difference from two samples/conditions. I believe it estimate the variance from two different samples/conditions. but in this case, the hypothesis that most genes have similar expression level is not necessarily hold. So I think that is the reason why the variance is larger in cuffdiff result.

What's your opinion?


All times are GMT -8. The time now is 06:04 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.