Batch correction for a sets of genes
1
0
Entering edit mode
Fereshteh ▴ 30
@fereshteh-15803
Last seen 10 weeks ago
United Kingdom

Hi,

I have 2545 and 1402 genes in Oncology Biomarker Panel and Precision Immuno-oncology Panel respectively. I have 719 common genes between two panels. I will need to merge raw read counts from these panels for differential expression (we know experimental condition, chemistry for samples from the same patients makes this reasonable). However, I thought to do differential expression analysis for 719 common by taking panel effect as batch like below

    condition   batch
A1  treatment   1
A2  treatment   1
A3  treatment   1
A4  treatment   1
A5  treatment   1
A6  treatment   1
A7  control     1
A8  control     1
A9  control     1
A10 control     1
A11 control     1
A12 control     1
B1  treatment   2
B2  treatment   2
B3  treatment   2
B4  treatment   2
B5  treatment   2
B6  treatment   2
B7  control     2
B8  control     2
B9  control     2
B10 control     2
B11 control     2
B12 control     2 

dds <- DESeqDataSetFromMatrix(countData = countData,colData = colData,design = ~ batch + condition)

neu.dds.LRT <- DESeq(dds,betaPrior=FALSE, test="LRT",
                     full=~ batch + condition, reduced=~batch)
dpsc.res.LRT <- results(neu.dds.LRT , contrast=c("condition", "treatment", "control"))
)

So that I would have differentially expressed genes among 719 common genes for treatment Vs control considering panel effect. Then I should combine uncommon genes between panels and do differential expression analysis separately (no longer with batch correction). I would have two sets of results. My question is , is there any way to have a set of differential expression analysis with two panels raw read counts although they have 719 common enes? Or, is these any way to combine design after making dds to unified my results?

Thanks for any assistance

cancer deseq2 r rnaseq • 429 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 17 hours ago
United States

What technology are you using here?

I'm not sure at all about the idea of merging counts from different platforms, and whether taking the intersection and allowing for a shift is sufficient to deal with the systematic differences.

ADD COMMENT
0
Entering edit mode

Thank you, they are Oncology Biomarker Panel and Precision Immuno-oncology Panel in HTG assay

ADD REPLY
1
Entering edit mode

I don’t know anything about the technology so I’d recommend speaking with someone more familiar.

ADD REPLY

Login before adding your answer.

Traffic: 281 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