deseq2 and dexseq for multiple tissue from wildtype only
1
0
Entering edit mode
badribio • 0
@badribio-12724
Last seen 5.4 years ago

I have to analyze following experimental design 

tissue wild type

sequence library

tissue1

paired-end

tissue1

paired-end

tissue1

paired-end

tissue2

paired-end

tissue2

paired-end

tissue2

paired-end

tissue3

paired-end

tissue3

paired-end

tissue3

paired-end

tissue4

paired-end

tissue4

paired-end

tissue4

paired-end

tissue5

paired-end

tissue5

paired-end

tissue5

paired-end

tissue6

paired-end

tissue6

paired-end

tissue6

paired-end

tissue7

paired-end

tissue7

paired-end

tissue7

paired-end

tissue8

paired-end

tissue8

paired-end

tissue8

paired-end

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. 

deseq2 dexseq rnaseq • 1.2k views
ADD COMMENT
2
Entering edit mode
@mikelove
Last seen 1 day ago
United States
Take a look at the software vignette, in DESeq2 you can use the 'contrast' argument to compare any two levels of a factor for example.
ADD COMMENT
0
Entering edit mode

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? 

ADD REPLY
1
Entering edit mode
Just use ~tissue. You can't add redundant information to a linear model
ADD REPLY
0
Entering edit mode

Thanks, is there a sample code that you would have come across in the past? It will help me understand better appreciate it!

ADD REPLY

Login before adding your answer.

Traffic: 734 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6