Hope for help about quantro and qsmooth-library size normalization
0
0
Entering edit mode
Apr • 0
@apr-21255
Last seen 5.4 years ago
PKU

I was hoping to follow quantro and qsmooth methods in my RNA-seq data normalization. My problem is that for my data, the sequence depth (total read counts) of the different groups were not consistent (e.g. clean base number of two of these groups were respect about 5G and 7G) and the quantro results showed there were global differences between these groups. And after qsmooth, the global differences between groups were preserved and it contained a reflection of their total read counts differences. So my question is that in this case, whether an extra total read counts normalization among groups is needed to be done before using the quantro or qsmooth in the RNA-seq data analysis. Thanks a lot!

normalization quantro q'smooth library size normalization • 847 views
ADD COMMENT
0
Entering edit mode

Do you think differences in the sequencing depth are related to biological or technical variation?

If it's technical, quantro actually performs an ANOVA to test if the medians of the distributions are different across groups. There is a parameter in quantro (useMedianNormalized = TRUE) to make the medians the same if the medians are different (https://github.com/stephaniehicks/quantro/blob/3fe77cae1c7f8ede0b1a1dda52e17981fd79fbcf/R/AllClasses.R#L196). The default is TRUE, but you can set this to false if you suspect the difference is due to biological variation.

If it is biological, you can also try using control genes (e.g. https://support.bioconductor.org/p/113059/#113061).

ADD REPLY

Login before adding your answer.

Traffic: 606 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