Question: Missing design when using Scran to create a DESeqDataSet
gravatar for berry14
5 months ago by
berry140 wrote:

The Scran vignette indicates the SCESet object can be convert to an object for DESeq2 analysis.

Once it converted, it's already an DESeqDataSet, yet without a design.  it's not obvious to me how to add a design to this DESeqDataSet.

Any suggestion is welcome. Thanks.

ADD COMMENTlink modified 5 months ago by Aaron Lun18k • written 5 months ago by berry140
gravatar for Aaron Lun
5 months ago by
Aaron Lun18k
Cambridge, United Kingdom
Aaron Lun18k wrote:

It seems like ?design would be a good place to start.

ADD COMMENTlink written 5 months ago by Aaron Lun18k

In DESeq2, the design is added when setting up the object using various methods, such as

DESeqDataSetFromTximport, DESeqDataSetFromHTSeqCount, DESeqDataSetFromMatrix, or DESeqDataSet

When Scran convert SCESet to DESeqDataSet, the count matrix is converted directly to DESeqDataSet. I was not able to find an option to add design?

ADD REPLYlink written 5 months ago by berry140

You can just add it afterwards with design(dds) <- formula(...), assuming your factors are in colData. I don't want to add more arguments to convertTo than the bare minimum, otherwise the interface becomes unwieldy (a handful of package-specific arguments for each output class adds up to a lot of arguments in the final function). It is easier for everyone if package-specific operations are performed outside of convertTo on the returned object.

ADD REPLYlink modified 5 months ago • written 5 months ago by Aaron Lun18k

Yes, you are right. That works. Thanks!

And I agree the bare minimum is a cleaner solution.

ADD REPLYlink written 5 months ago by berry140
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: 121 users visited in the last hour