Biostar Beta. Not for public use.
boxplot in edgeR
0
Entering edit mode
2.3 years ago
Björn • 40
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 ?

ADD COMMENTlink
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 REPLYlink
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink
1
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

Urrgh... maybe in Irish

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode
2.2 years ago
theobroma22 ♦ 1.1k

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

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3