I have to analyze following experimental design
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I basically have various type of tissues from wildtype mouse ( 8 tissues numbered 1-8) I would like to check for gene expression and junction usage differences between the tissues using deseq2 and dexseq respectively, what should be the right approach for design model and what kind of functions should i be using to achieve this? (apologies for such an open ended question) historically have always performed wt vs mt , any help appreciated.
so if i understand it correct i can add group as a factor for the design something like this
tissue wild type
group
tissue1
group1
tissue1
group1
tissue1
group1
tissue2
group2
tissue2
group2
tissue2
group2
tissue3
group3
tissue3
group3
tissue3
group3
tissue4
group4
tissue4
group4
tissue4
group4
tissue5
group5
tissue5
group5
tissue5
group5
tissue6
group6
tissue6
group6
tissue6
group6
tissue7
group7
tissue7
group7
tissue7
group7
tissue8
group8
tissue8
group8
tissue8
group8
and `design = ~ tissue + group`
is this correct?
Thanks, is there a sample code that you would have come across in the past? It will help me understand better appreciate it!