SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 Anomaly/Issue? The_Corsair Bioinformatics 7 03-02-2015 07:07 PM
BluePippin issue matth431 General 0 11-25-2014 03:13 AM
maybe a Gviz issue??? emolinari Bioinformatics 2 06-24-2013 07:26 AM
Picard Issue ercfrtz Bioinformatics 1 06-05-2012 10:00 AM
CisGenome Issue Nicholas_ Bioinformatics 0 10-10-2011 02:24 AM

Reply
 
Thread Tools
Old 04-03-2017, 07:43 AM   #1
chammer
Junior Member
 
Location: Lausanne

Join Date: Feb 2014
Posts: 5
Question DESeq2 collinearity issue

Hi everybody,

I have been trying to use DESeq2 for analyzing RNA-seq data, and ran into a problem.

The conditions table looks like this:

Code:
sample	donor	virus	vpu
DonorA1_01	A1	none	mock
DonorA1_02	A1	CH293	wt
DonorA1_03	A1	CH293	stop
DonorA1_04	A1	CH293	R50K
DonorA1_05	A1	CH293	teth_count
DonorA1_06	A1	CH077	wt
DonorA1_07	A1	CH077	stop
DonorA1_08	A1	CH077	R50K
DonorA1_09	A1	CH077	teth_count
DonorA1_10	A1	STC01	wt
DonorA1_11	A1	STC01	stop
DonorA1_12	A1	STC01	R50K
DonorA1_13	A1	STC01	teth_count
DonorX_01	X	none	mock
DonorX_02	X	CH293	wt
DonorX_03	X	CH293	stop
DonorX_04	X	CH293	R50K
DonorX_05	X	CH293	teth_count
DonorX_06	X	CH077	wt
DonorX_07	X	CH077	stop
DonorX_08	X	CH077	R50K
DonorX_09	X	CH077	teth_count
DonorX_10	X	STC01	wt
DonorX_11	X	STC01	stop
DonorX_12	X	STC01	R50K
DonorX_13	X	STC01	teth_count
DonorY_01	Y	none	mock
DonorY_02	Y	CH293	wt
DonorY_03	Y	CH293	stop
DonorY_04	Y	CH293	R50K
DonorY_05	Y	CH293	teth_count
DonorY_06	Y	CH077	wt
DonorY_07	Y	CH077	stop
DonorY_08	Y	CH077	R50K
DonorY_09	Y	CH077	teth_count
DonorY_10	Y	STC01	wt
DonorY_11	Y	STC01	stop
DonorY_12	Y	STC01	R50K
DonorY_13	Y	STC01	teth_count
DonorZ_01	Z	none	mock
DonorZ_02	Z	CH293	wt
DonorZ_03	Z	CH293	stop
DonorZ_04	Z	CH293	R50K
DonorZ_05	Z	CH293	teth_count
DonorZ_06	Z	CH077	wt
DonorZ_07	Z	CH077	stop
DonorZ_08	Z	CH077	R50K
DonorZ_09	Z	CH077	teth_count
DonorZ_10	Z	STC01	wt
DonorZ_11	Z	STC01	stop
DonorZ_12	Z	STC01	R50K
DonorZ_13	Z	STC01	teth_count
When I specify the model for differential expression analysis as dds <- DESeqDataSetFromTximport(txi, samples, ~vpu+donor+virus), I get an error message:

Error in checkFullRank(modelMatrix) :
the model matrix is not full rank, so the model cannot be fit as specified.
One or more variables or interaction terms in the design formula are linear
combinations of the others and must be removed.

Same problem if the model is “vpu + virus” only.

I understand that this is because of collinearity among the variables but I am not sure how to resolve the issue.

Any help would be highly appreciated!

Thanks!
Chris
chammer is offline   Reply With Quote
Old 04-04-2017, 01:08 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

The problem is that "virus==none" is the same as "vpu==mock".
dpryan is offline   Reply With Quote
Old 04-04-2017, 02:28 AM   #3
chammer
Junior Member
 
Location: Lausanne

Join Date: Feb 2014
Posts: 5
Default

Thanks a lot!
With respect to the most important comparisons, which is VPU "wt" vs. "stop", leaving out the non-infected samples works.
But if I want to compare non-infected vs. infected samples, how could this be resolved? Using an additional dummy variable?
How would that look like?

Thanks!
Chris
chammer is offline   Reply With Quote
Old 04-04-2017, 02:34 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

My guess (since I haven't a clue about the background to your experiment) is that you want vpu="Wt" for the mock infected samples. The combination of vpu and virus would then be distinct.

*Edit*: Alternatively, set the virus="none" samples to what they're mock infected with (e.g., CH293), which I suspect will be clearer.
dpryan is offline   Reply With Quote
Old 04-04-2017, 11:30 AM   #5
chammer
Junior Member
 
Location: Lausanne

Join Date: Feb 2014
Posts: 5
Default

Perfect, thanks a lot for your help, Devon!

Best,
Chris
chammer is offline   Reply With Quote
Reply

Tags
collinearity, deseq2, rnaseq

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 07:54 PM.


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