i am a newbie in RNA seq data analysis. I want to check a gene as DEG in a dataset of RNA-chip seq experiment. The data provided is in the form of a single column for each treatment type and lists the expression level of each gene normalized to transcripts per million (TPM).i need help on using what type of analysis using R to find out the DEGs.
The snapshot of data is like this.
Gene EntrezID Normal_TPM Diabetes_TPM
Apoa1 11806 14668.15 2875.06
Mup3 17842 9992.58 1697.63
Serpina3k 20714 8031.3 2849.67
Alb 11657 6801.26 6912.08
Cyp2e1 13106 6580.8 7816.79
Thanks in advance
edgeR works with raw counts, so maybe EBseq could be a better advice for TPMs?
http://deweylab.biostat.wisc.edu/rsem/README.html
Hi Govardhan, I think you're mixing up CPM (counts per million) with TPM (transcripts per million). See here how it's computed.
yeah, so you can get TPM formula here then.