normalization RNASeq Data for QTL analysis
0
1
Entering edit mode
8 months ago
Linda ▴ 10

Hi all,

I am looking for the best way to normalize RNA-Seq data generated using the 3’prime technology of Lexogen. Due to this technology, I do not need to include gene length.

However, the data should follow a normal distribution, as I will proceed with linear modeling in matrixeQTL.

I've read quite a lot about RNASeq normalization, but I feel like there are many options and opinions and at the same time not so much info about my concrete situation. For example, here https://pubmed.ncbi.nlm.nih.gov/31861985/ TPM seems to perform best (but I do not need to account for gene length), while log10 was evaluated as not good regarding biological signal.

Starting from raw counts, I see the options now to start either with

  • DESeq / TMM / quantile

to normalize between samples.

But then I think I need further transformation to make the data "normal". For example

  • log, vst or inverse rank transformation

I would be thankful for any advice from you!

RNA-Seq R edgeR DESeq matrixQTL • 362 views
ADD COMMENT
0
Entering edit mode

I found some advice on matrixeQTL faq:

If you go to the "Outliers in expression. Quantile normalization and Kruskal-Wallis test." paragraph, they also suggest some R code for the normalization.

ADD REPLY

Login before adding your answer.

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