Normalization using the quantile function
1
0
Entering edit mode
3.9 years ago
Leite ★ 1.2k

Hello everyone,

I normalized a data using the quantile function in the Bioconductor package limma, boxplots for each group before (A) and after (B) the process this does not look like what I usually see.

Could anyone, please, explain to me why by that black part passed to down of the boxplot in B?

Is that correct?

Rplot
image uploader

Best Regards,

Leite

R normalization quantile • 1.8k views
ADD COMMENT
3
Entering edit mode
3.9 years ago

Through quantile normalisation, which I assume you used, it's virtually impossible for the highly-expressed outliers in A to be the same as the lowly-expressed outliers in B. I'm not sure if that's what you were implying?

It may be due to the long, lower tail that's present in one of your samples (9 samples from the right in A), which, through quantile normalisation, would affect the distribution in all other samples.

I neither believe that this is a problem in any way. I have seen elsewhere how, after quantile normalisation, lowly-expressed outliers will exist. Before deciding if this is a problem, you should take a look at:

  • PCA bi-plot
  • density plots of normalised counts
  • pairwise scatterplot matrices
  • quantile-quantile plots after differential expression analysis
ADD COMMENT
1
Entering edit mode

Hey Kevin,

Thank you so much! I

boxplot(log(as.matrix(x)),las=2,ylab="Log2(Intensity)")

I checked my code and found an error, I used this function to create the boxplot, but my sample x was already in the log function, so it was as if I was log twice.

Now I gotta fix the error.

Best regards,

Leite

ADD REPLY

Login before adding your answer.

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