SEQanswers Help with complicated DESeq2 design formula
 Register FAQ Members List Calendar Search Today's Posts Mark Forums Read

 Similar Threads Thread Thread Starter Forum Replies Last Post bryand RNA Sequencing 3 07-24-2015 02:22 AM pinuscatcus Bioinformatics 1 04-10-2015 05:53 AM dmosk Bioinformatics 35 01-05-2015 02:42 PM Wei-HD Bioinformatics 11 09-04-2014 11:13 AM KristenC RNA Sequencing 1 05-29-2014 11:05 AM

 12-01-2015, 09:35 AM #1 aroomacanvas Junior Member   Location: USA Join Date: Dec 2015 Posts: 1 Help with complicated DESeq2 design formula Hello. I'm familiar with DESeq2 but the most complex design equations I've used are having blocked pairs. I need to perform a more complicated design but I'm having trouble formulating the equation. I have a design table that I've read in that contains all the relevant information. Here are the specifications: condition = sample/control (what I'm testing over) week = a timepoint comparison of when the sample was generated grade = a comparison of what grade of disease the sample has (so 4 > 3 > 2 > 1 > 0 ie control) There are also 3 possible organ sites but some samples are a combination, so I've split it into 3 binary variables: is_SiteA, is_SiteB, and is_SiteC, all of which are 0 or 1 This is what the table looks like: Code: ```sample condition week grade is_SiteA is_SiteB is_SiteC SAMPA sample 14 3 0 1 0 SAMPB control 7 0 0 0 1 SAMPC control 14 0 1 0 1``` Right now my equation is this: ddsHTSeq <- DESeqDataSetFromMatrix(counts_table, sample_table, design = ~ week : grade : is_GI + is_Liver + is_Skin + condition) Is this correct at all? I really appreciate any help. FYI I do have more than 3 samples, before anyone asks! Last edited by aroomacanvas; 12-01-2015 at 11:03 AM.

 Tags deseq2