I am looking at the output of cuffdiff. It has columns:test_id, gene_id, gene, locus, sample_1, sample_2, status, value_1, value_2, log2(fold_change), test_stat, p_value, q_value, significant. Can anybody explain how test_stat, p value, q value calculated ? what is the meaning of p value when i have only one replicate ? Also, in some cases the values of value_1 and value_2 are 0 making fold change infinite. How should I handle these cases ? i tried substituting 0's with small value 0.001 but there are some values that are as less as 0.0000160004. if i replace 0's with 0.0000000001 fold change becomes way high in comparison to other cases. How do I handle these cases ? I don't think I can ignore them because if a gene is not expressed at all in one experiment and expressed in other, it must be important.
thanks in advance for the help!!!
cuffdiff FAQ has moved