EdgeR: is more than one control per batch an issue?
1
0
Entering edit mode
Ekarl2 ▴ 80
@ekarl2-7202
Last seen 6.4 years ago
Sweden

The current dataset I am looking at has seven environments (control and three different concentrations of exposure A separately and exposure B separately) and each exposure/concentration combination has three replicates (one in each batch).

My goal is to compare each of the six experimental conditions with the control. Here my plan is to use one treatment factor with seven levels and one batch factor with three levels.

The complication, however, is that there are more than one control per batch. Batch 1 has one control, whereas batch 2 and 3 have two controls each. In those batches that has more than one control, these controls are _not_ technical replicates.

Will this be an issue for an edgeR-based differential expression analysis? It appears that "control" is unevenly represented in the different batches. Or can the kind of analysis I described above handle this? Are there any other ways can this be handled?

edgeR batch effect • 859 views
ADD COMMENT
3
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 6 hours ago
The city by the bay

There is no issue, the GLM machinery in edgeR will handle such designs without any problems.

ADD COMMENT
0
Entering edit mode

A quick follow-up. The reason for the fact that batch 1 only has one control is that one of them failed. To compensate, they made another control in batch 4, but that is the only sample in that batch. My instinct is this sample should be excluded because there is only one sample for batch 4, so there is no way to know what the effects of being in batch 4 is.

What do you think about this reasoning?

ADD REPLY
1
Entering edit mode

You are correct. The only way to include such information would be to use duplicateCorrelation with voom, but that has its own issues. To me, it doesn't seem worth the trouble just to include one extra sample.

ADD REPLY

Login before adding your answer.

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