Description Usage Arguments Details Value Examples
adjust_abundance() takes as input a 'tbl' formatted as | <SAMPLE> | <TRANSCRIPT> | <COUNT> | <...> | and returns a 'tbl' with an edditional adjusted abundance column. This method uses scaled counts if present.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 | adjust_abundance(
.data,
.formula,
.sample = NULL,
.transcript = NULL,
.abundance = NULL,
log_transform = TRUE,
action = "add",
...
)
## S4 method for signature 'spec_tbl_df'
adjust_abundance(
.data,
.formula,
.sample = NULL,
.transcript = NULL,
.abundance = NULL,
log_transform = TRUE,
action = "add",
...
)
## S4 method for signature 'tbl_df'
adjust_abundance(
.data,
.formula,
.sample = NULL,
.transcript = NULL,
.abundance = NULL,
log_transform = TRUE,
action = "add",
...
)
## S4 method for signature 'tidybulk'
adjust_abundance(
.data,
.formula,
.sample = NULL,
.transcript = NULL,
.abundance = NULL,
log_transform = TRUE,
action = "add",
...
)
## S4 method for signature 'SummarizedExperiment'
adjust_abundance(
.data,
.formula,
.sample = NULL,
.transcript = NULL,
.abundance = NULL,
log_transform = TRUE,
action = "add",
...
)
## S4 method for signature 'RangedSummarizedExperiment'
adjust_abundance(
.data,
.formula,
.sample = NULL,
.transcript = NULL,
.abundance = NULL,
log_transform = TRUE,
action = "add",
...
)
|
.data |
A 'tbl' formatted as | <SAMPLE> | <TRANSCRIPT> | <COUNT> | <...> | |
.formula |
A formula with no response variable, representing the desired linear model where the first covariate is the factor of interest and the second covariate is the unwanted variation (of the kind ~ factor_of_intrest + batch) |
.sample |
The name of the sample column |
.transcript |
The name of the transcript/gene column |
.abundance |
The name of the transcript/gene abundance column |
log_transform |
A boolean, whether the value should be log-transformed (e.g., TRUE for RNA sequencing data) |
action |
A character string. Whether to join the new information to the input tbl (add), or just get the non-redundant tbl with the new information (get). |
... |
Further parameters passed to the function sva::ComBat |
maturing
This function adjusts the abundance for (known) unwanted variation. At the moment just an unwanted covariate is allowed at a time using Combat (DOI: 10.1093/bioinformatics/bts034)
Underlying method: sva::ComBat(data, batch = my_batch, mod = design, prior.plots = FALSE, ...)
A 'tbl' with additional columns for the adjusted counts as '<COUNT COLUMN>_adjusted'
A 'tbl' with additional columns for the adjusted counts as '<COUNT COLUMN>_adjusted'
A 'tbl' with additional columns for the adjusted counts as '<COUNT COLUMN>_adjusted'
A 'tbl' with additional columns for the adjusted counts as '<COUNT COLUMN>_adjusted'
A 'SummarizedExperiment' object
A 'SummarizedExperiment' object
1 2 3 4 5 6 7 8 9 | cm = tidybulk::counts_mini
cm$batch = 0
cm$batch[cm$sample %in% c("SRR1740035", "SRR1740043")] = 1
res =
cm %>%
tidybulk(sample, transcript, count) %>%
identify_abundant() %>%
adjust_abundance( ~ condition + batch )
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.