Question: How to use TPM in plotMA in DESeq2
gravatar for georgewwp
26 days ago by
georgewwp0 wrote:

Hi all,

I'm doing RNAseq analysis with DESeq2. When making MA plot, I want to position each gene on the x-axis using its mean abundance from one condition instead of rowmeans from counts(dds) across samples. This way, x-axis would carry information I'm interested in. 

I quantified the transcripts using salmon and summarized at the gene level using tximport. 

I'm thinking I could do txi <- tximport(files, type = "salmon", tx2gene = tx2gene, countsFromAbundance = "lengthScaledTPM") and use txi$counts. Is this a good option?

Is there a way I can access normalized TPM from within dds so that I can use the plotMA function?

I discovered counts(dds) do not change after dds <- estimateSizeFactors(dds). How do I access the counts across samples after normalization? 

Thank you!


ADD COMMENTlink modified 22 days ago by Michael Love20k • written 26 days ago by georgewwp0
gravatar for Michael Love
22 days ago by
Michael Love20k
United States
Michael Love20k wrote:

You can use txi$abundance for the TPM measure. plotMA is a very simple function, you can probably most easily just write your own function that plots mean TPM on the x and log2FoldChange on the y, coloring points a different color if padj < threshold.

ADD COMMENTlink written 22 days ago by Michael Love20k

Thank you, Michael!

ADD REPLYlink written 20 days ago by georgewwp0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 390 users visited in the last hour