fetch_assay_data | R Documentation |
The (fetch|with)_assay_data
functions are some of the main workhose
functions of the facile ecosystem. These calls enable you to retrieve
raw and noramlized assay data from a FacileData container.
fetch_assay_data(
x,
features,
samples = NULL,
assay_name = ndefault_assay(x),
normalized = FALSE,
batch = NULL,
main = NULL,
as.matrix = FALSE,
drop_samples = TRUE,
...,
subset.threshold = 700,
aggregate = FALSE,
aggregate.by = "ewm",
verbose = FALSE
)
## S3 method for class 'facile_frame'
with_assay_data(
x,
features,
assay_name = NULL,
normalized = TRUE,
aggregate = FALSE,
aggregate.by = "ewm",
spread = TRUE,
with_assay_name = FALSE,
...,
verbose = FALSE,
.fds = fds(x)
)
x |
A |
features |
a feature descriptor (data.frame with assay and feature_id columms) |
samples |
a samples descriptor |
assay_name |
the name of the assay to fetch data from. Defaults to the
value of |
normalized |
return normalize or raw data values, defaults to |
batch |
The column names in |
main |
The name of a covaraite in |
as.matrix |
by default, the data is returned in a long-form tbl-like
result. If set to |
... |
parameters to pass to normalization methods |
subset.threshold |
sometimes fetching all the genes is faster than trying to subset. We have to figure out why that is, but I've previously tested random features of different lengths, and around 700 features was the elbow. |
aggregate.by |
do you want individual level results or geneset scores? Use 'ewm' for eigenWeightedMean, and that's all. |
.fds |
A |
feature_ids |
character vector of feature_ids |
with_symbols |
Do you want gene symbols returned, too? |
fetch_assay_data(x, ...)
will return the data in long form.
with_assay_data(x, ...)
is most typically used when you already have
a dataset x
(a facile_frame
) that you want to decorate with more assay
data. The assay data asked for will be appended on to x
in wide format.
Because fetch
is (most often) used at a lower level of granularity,
normalize
is by default set to FALSE
, while it is set to TRUE
in
with_assay_data
.
A tibble
(lazy or not) with assay data.
a tbl-like result
When normalized data is returned, we assume these data are log-like, and you
have the option to regress out batch effects using our
remove_batch_effect()
wrapper to limma::removeBatchEffect()
.
samples <- exampleFacileDataSet() |>
filter_samples(indication == "BLCA", sample_type == "tumor")
features <- c(PRF1='5551', GZMA='3001', CD274='29126')
dat <- with_assay_data(samples, features, normalized = TRUE, batch = "sex")
dat <- with_assay_data(samples, features, normalized = TRUE,
batch = c("sex", "stage"))
dat <- with_assay_data(samples, features, normealized = TRUE,
batch = c("sex", "stage"), main = "sample_type")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.