Entering edit mode
Dear,
May I know the difference between the txi$counts and txi$abundance when the countFromAbundance is set to true(scaledTPM or lengthScaledTPM or dtuScaledTPM)? If set to xxxTPM, the txi$counts is still the raw counts value, right? As it's so close to the value when countFromAbundance is set to no.
Thanks a lot!
Thank you very much!
After you use one of the xxxTPM, if you check the txi$counts and txi$abundance, you can find relatively big difference. Then would you choose txi$counts or txi$abundance to export as xxxTPM value/matrix?
Best,
I cannot find any suggestion (like anywhere) that you should use txi$abundance 'to export as xxxTPM value/matrix'. Can you point to any documentation that says you should use anything but the counts list item?
In other words, the different countsFrom methods are intended to provide counts that are adjusted using the abundance measurements, in order to account for length biases and whatnot. There is no reason for you to do anything directly with the abundance values.
Make it easy on yourself, use the counts assay together with the length assay and calculate your TPMs as suggested: DESeq2: Is it possible to convert read counts to expression values via TPM and return these values?