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

 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 Forum Rules

All times are GMT -8. The time now is 11:33 PM.