Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 - paired formula design with replicates bryand RNA Sequencing 3 07-24-2015 02:22 AM
DESeq2 help.. design pinuscatcus Bioinformatics 1 04-10-2015 05:53 AM
DESeq2 model matrix formula dmosk Bioinformatics 35 01-05-2015 02:42 PM
Deseq2 design Wei-HD Bioinformatics 11 09-04-2014 11:13 AM
Paired design versus unpaired design in DESeq2 KristenC RNA Sequencing 1 05-29-2014 11:05 AM

Thread Tools
Old 12-01-2015, 09:35 AM   #1
Junior Member
Location: USA

Join Date: Dec 2015
Posts: 1
Question 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:
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.
aroomacanvas is offline   Reply With Quote


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 04:23 PM.

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