Question: Looking for code showing how estimateDispersions works.
gravatar for jarrod.ellingson
5 weeks ago by
jarrod.ellingson0 wrote:

I'm trying to understand better how DESeq2 gets from my raw data to the results. For most functions in R, I type the function name in the console and get the full code used. There are a few important functions in DESeq2 that are lacking this information, and I'm not finding it from a simple Google search. Would greatly appreciate it if someone can point me in the right direction for this information.

As an example, trying to get the code written for estimateDispersions produces the following (note that running BiocGenerics::estimateDispersions gives the same results). 

> estimateDispersions
standardGeneric for "estimateDispersions" defined from package "BiocGenerics"

function (object, ...)
<bytecode: 0x1080286f8>
<environment: 0x108020a40>
Methods may be defined for arguments: object
Use  showMethods("estimateDispersions")  for currently available ones.

ADD COMMENTlink modified 5 weeks ago by Dario Strbenac1.4k • written 5 weeks ago by jarrod.ellingson0
gravatar for Dario Strbenac
5 weeks ago by
Dario Strbenac1.4k
Dario Strbenac1.4k wrote:

The getMethod function does that for S4 methods. The code of estimateDispersions is displayed in the R console by executing getMethod(estimateDispersions, "DESeqDataSet"). The options you could specify for the second argument of getMethod are listed in the output of showMethods("estimateDispersions"). For estimateDispersions, DESeqDataSet is the only viable option.

ADD COMMENTlink written 5 weeks ago by Dario Strbenac1.4k

Ah--thank you. I had used showMethods("estimateDispersions"), based on what I was seeing in the console, but was unaware of getMethod. This gives me what I need for this issue, and I'll keep it in mind for future use.

ADD REPLYlink written 5 weeks ago by jarrod.ellingson0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 302 users visited in the last hour