Question

How to reassign DESeqDataSet design to include newly added column to columnData

0

Entering edit mode

rnavabjj • 0

@user-25125

Last seen 3.0 years ago

United States

I would like to test for group-specific effects in my dataset. This contains a group, time and individual (ID) factors as is in the DESEQ2 Vignette under "Group-specific condition effects, individuals nested within groups". After adding a 'id.n' column to my dds, the vignette states:

"Now we can reassign our DESeqDataSet a design of ~ grp + grp:ind.n + grp:cnd, before we call DESeq. This new design will result in the following model matrix:"

It provides no information on reassigning the formula. I tried the following but the formula is linked to ibjects and this replaces it with language.

dds@design[[2]] <- 'group + group:ind.n+group:time'

Thank you for you help.

DESeq2 • 3.0k views

ADD COMMENT • link updated 3.3 years ago by Michael Love 43k • written 3.9 years ago by rnavabjj • 0

score 2 · Accepted Answer · 2021-03-26

2

Entering edit mode

Michael Love 43k

@mikelove

Last seen 1 day ago

United States

In Bioconductor, often there are dedicated "setter" and "getter" functions:

design(dds) <- ...

Another pointer, there is a way to see what functions are available for a specific object _and_ defined in a specific package, might be useful for you:

> showMethods(classes="DESeqDataSet", where=getNamespace("DESeq2"))

...

Function: design<- (package BiocGenerics)
object="DESeqDataSet", value="formula"
object="DESeqDataSet", value="matrix"

...

I was working at one point on trying to simplify retrieving this type of information, but haven't worked on it in a while:

https://rpubs.com/mikelove/pkgmethods

ADD COMMENT • link 3.9 years ago Michael Love 43k

0

Entering edit mode

Your timely response just saved my dissertation, Dr. Love. Thank you very much.

ADD REPLY • link 3.9 years ago rnavabjj • 0

0

Entering edit mode

Hi Dr. Love, I'm revisiting this data and am getting the "full model matrix is less than full rank" error when I try to call DEseq2.

My ColData:

> dds.primary@colData
DataFrame with 28 rows and 4 columns
                         time     group          ID     id.n
                     <factor>  <factor> <character> <factor>
trained.11post.txt       post   trained          11       11
trained.11pre.txt        pre    trained          11       11
trained.16post.txt       post   trained          16       16
trained.16pre.txt        pre    trained          16       16
trained.17post.txt       post   trained          17       17
...                       ...       ...         ...      ...
untrained.19pre.txt      pre  untrained          19       19
untrained.5.post.txt     post untrained           5       5 
untrained.5.pre.txt      pre  untrained           5       5 
untrained.6.post.txt     post untrained           6       6 
untrained.6.pre.txt      pre  untrained           6       6

My design:

design(dds.primary)<- ~ group + group:id.n + group+time

ADD REPLY • link 3.3 years ago rnavabjj • 0

0

Entering edit mode

Sorry, I don't have time to help with statistical consultation on the support site. I have to save my time for software related questions. I'd recommend to find a local statistician or someone familiar with linear models in R.

ADD REPLY • link 3.3 years ago Michael Love 43k