Question: appropriate design formula
0
gravatar for mr6518849
15 days ago by
mr65188490
mr65188490 wrote:

Hello, I am not sure I am building the proper design formula for the question I want to answer

I have the following samples with three factors; clone, the structure and the condition.

clone structure diabetic
1 07        2D      Dia
2 21        2D     Ctrl
3 23        2D      Dia
4 32        2D     Ctrl
5 34        2D      Dia
6 43        2D     Ctrl
7 07        3D      Dia
8 21        3D     Ctrl
9 23        3D      Dia
10 32        3D     Ctrl
11 34        3D      Dia
12 43        3D     Ctrl

I want to pull the differentially expressed genes for structure (2D vs 3D) blocking by clone, and the one differentially expressed by condition (regardless of the structure variable).

I am not sure how to build the design in one call, therefore I have produced the first object to pull 2D vs 3D blocking by clone using the following.

des_structure <- DESeqDataSetFromHTSeqCount(sampleTable = samples,
                                       directory = "../data/htseq_geneCounts/",
                                       design= ~ clone+structure)

for the condition controlled by structure I have used

des_condidtion <- ddsHTSeq_diabetic <- DESeqDataSetFromHTSeqCount(sampleTable = samples,
                                                 directory = "../data/htseq_geneCounts/",
                                                 design= ~ structure+diabetic)

I am not sure whether this is the best way to build the designs.

deseq2 • 81 views
ADD COMMENTlink modified 14 days ago by Michael Love24k • written 15 days ago by mr65188490
Answer: appropriate design formula
1
gravatar for Michael Love
14 days ago by
Michael Love24k
United States
Michael Love24k wrote:

Yes those are the appropriate designs for DESeq2. For the second design, we ignore the fact that there are correlations between pairs of samples from the same clone, because DESeq2 cannot incorporate random effects. Another approach would be limma-voom which has a duplicateCorrelation function to inform about the sample pairings that are nested within the comparison of interest (diabetic status).

ADD COMMENTlink written 14 days ago by Michael Love24k

thank you for your feedback.

ADD REPLYlink written 14 days ago by mr65188490
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 341 users visited in the last hour