Comparing gene expression of specific genes between samples in DESeq2 / R
0
0
Entering edit mode
2.3 years ago
masicioglu • 0

Hi folks,

Using DESeq2, I've obtained log2 fold change (used to compare between two conditions), plot, heatmap ext. In addition to, I have two questions:

First, I've learnt TPM normalisation is suitable for compare gene expression level in the same sample but DESeq2 doesn't provide TPM normalisation code. How can I do TPM normalization in one sample for all interested genes from un-normalizate counts in R?

Secondly, I am wondering how can I compare one gene (for example p53) between a couple of samples(or two sample for example liver vs spleen)? Actually, heatmap gives this info but I want to graph comparing samples for one gene. Is it possible?

Could you help me, please?

Thanks

RNA-Seq DESeq2 tpm • 1.8k views
ADD COMMENT
0
Entering edit mode

I've learnt TPM normalisation is suitable for compare gene expression level in the same sample

Can you quote your source for that?

Secondly, I am wondering how can I compare one gene (for example p53) between a couple of samples(or two sample for example liver vs spleen)? Actually, heatmap gives this info but I want to graph comparing samples for one gene. Is it possible?

You could take the normalised counts ( counts(dds, normalized=TRUE) ) and plot them as barplots (as in A: Network/Pathway Analysis from Mass Spec data ): f

You could also take the regularised log counts and plot these as boxplots with scatter-plot overlays (as in A: Boxplot in ggplot2 ):

boxscatter

ADD REPLY

Login before adding your answer.

Traffic: 1275 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6