Search
Question: Looking for code showing how estimateDispersions works.
0
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, ...)
standardGeneric("estimateDispersions")
<bytecode: 0x1080286f8>
<environment: 0x108020a40>
Methods may be defined for arguments: object
Use  showMethods("estimateDispersions")  for currently available ones.

modified 5 weeks ago by Dario Strbenac1.4k • written 5 weeks ago by jarrod.ellingson0
2
5 weeks ago by
Dario Strbenac1.4k
Australia
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.