comparisson of RNA-seq FPKMs across samples
Entering edit mode
5.6 years ago

Dear BioStars community,

I am trying to "recycle" some RNA-seq data published as supplementary material of the publication "High-throughput screening using patient-derived tumor xenografts to predict clinical trial drug response" but I am not sure if I can really use it for my purpose.

I have a matrix of fpkm values per gene for 376 samples. For each gene in each sample, I would like to know whether the gene in this sample is highly expressed compared to the whole population (e.g obtain a Z-score for the expression of the gene in the sample given the expression of this gene in the whole population).

I have read that fpkm values depend on the sequencing depth of each sample and, therefore, cannot be directly compared across samples. Is there a way to somehow transform fpkm values into something that could be compared across samples?

I have never worked with RNA-seq data so any contribution will be more than welcome, even if you think it is a simple and basic comment or suggestion. I am comfortable programming with Python but I can also use R, perl or command line tools.

Many thanks in advance,


RNA-Seq fpkm • 2.3k views
Entering edit mode
4.0 years ago
FeiZhao • 0

You might as well try htseq-count and deseq2 to compare gene expression in different samples


Login before adding your answer.

Traffic: 1772 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6