boxplot in edgeR
1
0
Entering edit mode
6.8 years ago
Björn ▴ 110
par(mfrow=c(1,2))      
lcpm <- cpm(y2, log=TRUE)             
boxplot(lcpm, las=2, col=group$Sample, main="")       
title(main="A. Example: Unnormalised data",ylab="Log-cpm")

The boxplot I got was

https://ibb.co/bZ8RQn

Why don't I get normal boxplot after using above command ?

edgeR boxplot lcpm • 3.4k views
ADD COMMENT
1
Entering edit mode

Maybe you have too many outliers (the bars seem indeed to be squeezed at the bottom of the plot). What happens if you use:

boxplot(lcpm, las = 2, col = group$Sample, main = "", outline = FALSE)

?

ADD REPLY
0
Entering edit mode

Try a violin or joy plot rather than a box plot.

ADD REPLY
0
Entering edit mode

@kevin , I tried with your code but did not work, unfortunately.

ADD REPLY
0
Entering edit mode

Me? ...or do you mean Devon (sounds similar), or my colleague in Ireland, SelenoCysteine?

ADD REPLY
1
Entering edit mode

I'd like to hear the accent where Kevin/Devon sound like selenocysteine :)

ADD REPLY
0
Entering edit mode

Urrgh... maybe in Irish

ADD REPLY
0
Entering edit mode

Ancient gaelic indeed. Anyway, what do you mean it doesn't work? Which error message do you get?

ADD REPLY
0
Entering edit mode
6.8 years ago
theobroma22 ★ 1.2k

Try mean-centering your data with the log scale transformation.

ADD COMMENT

Login before adding your answer.

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