View source: R/statistic-udt.R
run_udt | R Documentation |
Calculates regulatory activities by using UDT.
run_udt(
mat,
network,
.source = source,
.target = target,
.mor = mor,
.likelihood = likelihood,
sparse = FALSE,
center = FALSE,
na.rm = FALSE,
min_n = 20,
seed = 42,
minsize = 5
)
mat |
Matrix to evaluate (e.g. expression matrix).
Target nodes in rows and conditions in columns.
|
network |
Tibble or dataframe with edges and it's associated metadata. |
.source |
Column with source nodes. |
.target |
Column with target nodes. |
.mor |
Column with edge mode of regulation (i.e. mor). |
.likelihood |
Deprecated argument. Now it will always be set to 1. |
sparse |
Deprecated parameter. |
center |
Logical value indicating if |
na.rm |
Should missing values (including NaN) be omitted from the
calculations of |
min_n |
An integer for the minimum number of data points in a node that are required for the node to be split further. |
seed |
A single value, interpreted as an integer, or NULL for random number generation. |
minsize |
Integer indicating the minimum number of targets per source. |
UDT fits a single regression decision tree for each sample and regulator,
where the observed molecular readouts in mat are the response variable and
the regulator weights in net are the explanatory one. Target features with
no associated weight are set to zero. The obtained feature importance from
the fitted model is the activity udt
of a given regulator.
A long format tibble of the enrichment scores for each source across the samples. Resulting tibble contains the following columns:
statistic
: Indicates which method is associated with which score.
source
: Source nodes of network
.
condition
: Condition representing each column of mat
.
score
: Regulatory activity (enrichment score).
Other decoupleR statistics:
decouple()
,
run_aucell()
,
run_fgsea()
,
run_gsva()
,
run_mdt()
,
run_mlm()
,
run_ora()
,
run_ulm()
,
run_viper()
,
run_wmean()
,
run_wsum()
inputs_dir <- system.file("testdata", "inputs", package = "decoupleR")
mat <- readRDS(file.path(inputs_dir, "mat.rds"))
net <- readRDS(file.path(inputs_dir, "net.rds"))
run_udt(mat, net, minsize=0)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.