DESeq2 color coding
1
0
Entering edit mode
ChIP-Tease • 0
@chip-tease-8339
Last seen 7.4 years ago
Germany

Hello everybody,

i was wondering, if it is possible to change the color coding in the DESeq2 plotMA function?

Since i cannot see how the function works, i don't know how to modify anything.

Additionally i would like to see how the plotMA works, which points to the same problem: I cannot see any code here. When i type the following, i can see a bit more code, but no possibility to change it.

> getMethod("plotMA","DESeqDataSet")
Method Definition:

function (object, ...)
{
    .local <- function (object, alpha = 0.1, main = "", ylim,
        ...)
    {
        res <- results(object, ...)
        plotMA.DESeqResults(res, alpha = alpha, main = main,
            ylim = ylim)
    }
    .local(object, ...)
}
<environment: namespace:DESeq2>

Thanks a lot, Alex

DESeq2 color R package • 4.1k views
ADD COMMENT
4
Entering edit mode
Mike Smith ★ 6.5k
@mike-smith
Last seen 22 hours ago
EMBL Heidelberg

The manual page for this function states that this is essentially a wrapper for the plotMA() function in the geneplotter package, and suggests you create a data.frame yourself and call that function directly if you want to change plotting parameters.

You can see from the code you've printed above, that you can create the data.frame with something like:

df <- results(dds)

And then you can plot it with different colour points like so:

plotMA(df, colNonSig = "blue")

If you want to read the help file for the geneplotter function to look at the option you can do ?geneplotter::plotMA

ADD COMMENT
0
Entering edit mode

Thanks a lot, that helped a lot!

ADD REPLY

Login before adding your answer.

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