SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 error : the model matrix is not full rank NicoBxl Bioinformatics 9 03-14-2017 06:10 AM
DESeq2 Simon Anders Bioinformatics 123 07-06-2015 01:45 AM
wired assembly result--abysmally short contigs chjp0632 Bioinformatics 6 09-09-2013 06:55 AM
Error Message in nbinomLRT in DESeq2 ToddB Bioinformatics 13 09-05-2013 06:22 AM
wired error profile mingkunli Bioinformatics 0 04-19-2013 01:46 AM

Reply
 
Thread Tools
Old 10-25-2013, 03:44 AM   #1
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default Wired DESeq2 error

Probably me who has done something wrong, but I dont see what.

First run was fine, but when I ran it again some days later I get this error.

> sampleFiles <- list.files(path="/Volumes/timemachine/HTseq_DEseq2",pattern="*.txt")
> sampleFiles <- sort(sampleFiles,decreasing=TRUE)
> status <- factor(c(rep("Healthy",26), rep("Diabetic",22)), levels=c("Healthy", "Diabetic"))
> timepoints = as.factor(c(2,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1))
> patients <- factor(c(88,87,86,83,82,81,75,72,71,70,13,08,01,88,87,86,83,82,81,75,72,71,70,13,08,01,79,77,76,73,67,62,61,55,53,21,04,79,77,76,73,67,62,61,55,53,21,04))
> sampleTable <- data.frame(sampleName = sampleFiles, fileName = sampleFiles, status=status, timepoints=timepoints, patients=patients)
> directory <- c("/Volumes/timemachine/HTseq_DEseq2/")
> design <- formula(~patients + timepoints*status)
> ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable, directory = directory, design = design)

Error in validObject(.Object) :
invalid class "DESeqDataSet" object: the model matrix is not full rank, i.e. one or more variables in the design formula are linear combinations of the others
sindrle is offline   Reply With Quote
Old 10-25-2013, 04:20 AM   #2
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

The error appears when:

timepoints = as.factor(c(1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2));

is changed to:

timepoints = as.factor(c(2,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1))
sindrle is offline   Reply With Quote
Old 10-25-2013, 04:36 AM   #3
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

I changed to:

timepoints = factor(c("After","After","After","After","After","After","After","After","After","After","After","After","After","Before","Before","Before","Before","Before","Before","Before","Before","Before","Before","Before","Before","Before","After","After","After","After","After","After","After","After","After","After","After","Before","Before","Before","Before","Before","Before","Before","Before","Before","Before","Before"), levels=c("Before","After"));

But when I change:

patients <-
as.factor(c(04,21,53,55,61,62,67,73,76,77,79,04,21,53,55,61,62,67,73,76,77,79, 01,08,13,70,71,72,75,81,82,83,86,87,88,01,08,13,70,71,72,75,81,82,83,86,87,88));

To:

patients <-
as.factor(c(88,87,86,83,82,81,75,72,71,70,13,08,01,88,87,86,83,82,81,75,72,71,70,13,08,01,79,77,76,73,67,62,61,55,53,21,04,79,77,76,73,67,62,61,55,53,21,04));

The error occurs again...

I really cant grasp this..
sindrle is offline   Reply With Quote
Old 10-25-2013, 04:40 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

This ends up being the same issue that you posted on the bioconductor email list, but there with edgeR (the solution is the same).
dpryan is offline   Reply With Quote
Old 10-25-2013, 05:20 AM   #5
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

Yes, Im sorry for that..

But I rally can't understand why:

timepoints = as.factor(c(1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2));

Changed to:

timepoints = as.factor(c(1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2));

Suddenly gives the error...
sindrle is offline   Reply With Quote
Old 10-25-2013, 05:29 AM   #6
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

Sorry for spamming my own thread.. But now it worked again!

Its indeed wired.
sindrle is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:07 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO