Collapses the columns in
object by summing within levels
of a grouping factor
groupby. The purpose of this function
is to sum up read counts from technical replicates to create an object
with a single column of read counts for each sample.
Note: by "technical replicates", we mean multiple sequencing runs of the same
library, in constrast to "biological replicates" in which multiple
libraries are prepared from separate biological units.
Optionally renames the columns of returned object with the levels of the
Note: this function is written very simply and
can be easily altered to produce other behavior by examining the source code.
a grouping factor, as long as the columns of object
optional, the names of each unique column in object. if provided,
a new column
whether to rename the columns of the returned object using the levels of the grouping factor
object with as many columns as levels in
This object has assay/count data which is summed from the various
columns which are grouped together, and the
colData is subset using
the first column for each group in
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
dds <- makeExampleDESeqDataSet(m=12) # make data with two technical replicates for three samples dds$sample <- factor(sample(paste0("sample",rep(1:9, c(2,1,1,2,1,1,2,1,1))))) dds$run <- paste0("run",1:12) ddsColl <- collapseReplicates(dds, dds$sample, dds$run) # examine the colData and column names of the collapsed data colData(ddsColl) colnames(ddsColl) # check that the sum of the counts for "sample1" is the same # as the counts in the "sample1" column in ddsColl matchFirstLevel <- dds$sample == levels(dds$sample) stopifnot(all(rowSums(counts(dds[,matchFirstLevel])) == counts(ddsColl[,1])))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.