Handling biological replicates in DESeq2
2
1
Entering edit mode
sk • 0
@fa435831
Last seen 1 day ago
Canada

I have 4 biological replicates of mRNA samples. For my dds, I am comparing type but have also added in the replicate factor to account for the biological replicates.

Is this correct? I have a column specifying the replicate number (1, 2, 3, or 4) per sample as they do not have the same names.

dds <- DESeqDataSetFromMatrix(
  countData = counts_data,
  colData = metadata,
  design = ~ Replicate + Type  # accounting for biological replicates
)
DESeq2 • 346 views
ADD COMMENT
1
Entering edit mode
ATpoint ★ 4.5k
@atpoint-13662
Last seen 44 minutes ago
Germany

Generally, you don't have to include the replication information since the software will consider all elements of dds$Type that have the same factor level as replicates for this group. It can sometimes make sense to include the replicate information, for example if every replicate is a different batch, so essentially replicate would be the same as batch, and there is evidence, for example in PCA, that batch is driving unwanted variation. If not, don't include.

ADD COMMENT
0
Entering edit mode
swbarnes2 ★ 1.4k
@swbarnes2-14086
Last seen 9 hours ago
San Diego

Go read the vignette. They don't include replicate.

ADD COMMENT

Login before adding your answer.

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