|Thread||Thread Starter||Forum||Replies||Last Post|
|DEXSeq with more than one condition||unicornich||Bioinformatics||1||03-21-2016 03:13 AM|
|DESeq2 and zero counts in one condition||gmatt||RNA Sequencing||1||01-09-2015 10:32 AM|
|Paired-test in DESeq or DESeq2||gcai||Bioinformatics||3||06-17-2014 01:57 AM|
|Paired design versus unpaired design in DESeq2||KristenC||RNA Sequencing||1||05-29-2014 11:05 AM|
|DESeq2: Difference between condition+type vs. 3 conditions||mevers||Bioinformatics||3||09-02-2013 01:17 AM|
|02-23-2015, 12:00 PM||#1|
Join Date: Mar 2011
DESeq2 paired plus condition
I am having some trouble getting my head around it today, its Monday.
Essentially, I have a design matrix like so:
So, according the the vignette, I need to do ~Patient + Condition. I really want condition differences, NOT patient differences.
I can view the results via results(slim_dds_1, contrast=c("Condition","3","2")). How do I know what was a patient difference? I suppose I could do the reduced formula and compare?
Last edited by bioBob; 02-23-2015 at 12:29 PM.
|02-24-2015, 07:09 AM||#3|
Join Date: Jul 2014
You can extract patient differences like so:
results(slim_dds_1, contrast=c("Patient", "58", "64")
|02-24-2015, 07:23 AM||#4|
Location: Freiburg, Germany
Join Date: Jul 2011
@bioBob: for future reference, we don't get easily notified when a post gets edited. For a speedier reply, post a comment after editing (or just post the edit as a new comment). Anyway, I suspect that fanli answered your question. If instead you really want to see what genes vary with patient in general (i.e., not between any individual pair of patients but generally across them), then you can use the LRT test with a reduced model of ~Condition.