Question: Error in DESeq "object" not found
0
gravatar for A
18 months ago by
A40
A40 wrote:

Hi all, 

 

Any help would be much appreciated. I have searched for this error and cannot find it anywhere but whats more weird is that I have never experienced this before. 

I have run the following code:

dds<-DESeqDataSetFromMatrix(countData = countdata, colData = sampledata, design = ~ Time_point + Organ)

dds1 <- DESeq(dds, reduced = Organ, test = "LRT")... I am now getting this error:

Error in DESeq(ddsAMIR, reduced = Organ, test = "LRT") : 
  object 'Organ' not found

This is really confusing as the first function (DESeqDataSetFromMatrix works absolutely fine and results check out. The DESeq function then throws up this error. 

A few notes: Sample data... Time_point and Organ are both factors. Countdata are integer and the column names and sample data are in precisely the correct order (again, this has worked before without any problem, so I am not sure what has provoked this issue)

Any help would be much appreciated! Many thanks!

 

 

deseq2 lrt reduced model • 261 views
ADD COMMENTlink modified 18 months ago by James W. MacDonald51k • written 18 months ago by A40
Answer: Error in DESeq "object" not found
0
gravatar for James W. MacDonald
18 months ago by
United States
James W. MacDonald51k wrote:

From ?DESeq:

reduced: for  test="LRT" , a reduced formula to compare against, i.e.,
          the full formula with the term(s) of interest removed.
          alternatively, it can be a model matrix constructed by the
          user

And from the examples:

 # an alternate analysis: likelihood ratio test
     ddsLRT <- DESeq(dds, test="LRT", reduced= ~ 1)
     resLRT <- results(ddsLRT)

Which seems pretty clear?

ADD COMMENTlink written 18 months ago by James W. MacDonald51k

I'd like the term of interest Organ removed though... So it should recognise that I am trying to do this.. Especially as it does in the full model design for the first function?

ADD REPLYlink written 18 months ago by A40

I agree that it would be sweet if R would be able to read your mind and know what you want it to do, but it doesn't work that way. Hence the help pages.

Do you know what is meant by the word 'formula' in the help page, particularly given the example that contains one? If not, you should look at ?formula, the help page for formulas.

ADD REPLYlink written 18 months ago by James W. MacDonald51k
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: 298 users visited in the last hour