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?