Question: Dynamic allocations possible when building DESeq object?
0
gravatar for bas_work
16 days ago by
bas_work0
bas_work0 wrote:

Hi,

As we are trying to build a pipeline in KNIME, I was wondering if it is at all possible to dynamically allocate a string to the design when building a DESeq. To clarify:

In dds <- DESeqDataSetFromMatrix(countData = counts, colData = sample, design= ~ condition)

we would like to allow other users to change the design parameter to for example "genotype", assuming that there is a column "genotype" in the sample metadata. There are ways to in KNIME to dynamically change input, but we do not get it to work in this specific case. Is it at all possible? Many thanks in advance for your input.

Kind regards, Bas

deseq2 • 72 views
ADD COMMENTlink modified 16 days ago by Michael Love26k • written 16 days ago by bas_work0
Answer: Dynamic allocations possible when building DESeq object?
0
gravatar for Michael Love
16 days ago by
Michael Love26k
United States
Michael Love26k wrote:

Are you going to call DESeqDataSetFromMatrix() within another function?

I'd recommend using a design matrix in the design slot in that case and not a formula.

If you create a formula variable inside of a function, it's difficult (impossible?) to prevent the function from bundling up all the objects into the formula (and so bloating the dataset). I've looked into this but I haven't found a solution, aside from pushing the formula creation to the user (so outside of the function and in the global environment). This is what happens in DESeq2, and so there is no dataset bloat.

In order to allow other developers to programmatically run DESeqDataSetFromMatrix() inside of a function, I modified DESeq2 a while ago so that design can be a matrix as created by model.matrix().

ADD COMMENTlink written 16 days ago by Michael Love26k

No, we are not trying to call DESeqDataSetFromMatrix() within another function and we have come to the conclusion that it is not really possible. Anyway, thanks for your thoughtful response.

Kind regards, Bas

ADD REPLYlink written 11 days ago by bas_work0
Please log in to add an answer.

Help
Access

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