Hi all,
I have two RNA-Seq samples. Calculating p-value for two samples is not possible. So, I want to calculate fold change for Sample1 vs Sample2.
Based on that I want to detect differential expressed genes.
Each sample has the following information.
tracking_id, class_code, nearest_ref_id, gene_id, gene_short_name, tss_id, locus length, coverage, FPKM, FPKM_conf_lo, FPKM_conf_hi, FPKM_status
Can anyone tell me how can I do that for that samples?
Is there any R package to calculate FC for two smaples?
1 answer
Look at Cuffdiff . The information you gave suggests that cufflinks was used. There are other tools to read in FPKM values but Cuffdiff would be a good starting point
http://cole-trapnell-lab.github.io/cufflinks/cuffdiff/#cuffdiff-input-files
Log in to answer this question.