Skip to contents

inferring cell-type specific gene regulatory network

Source: R/inferCSN.R
inferCSN.Rd

inferring cell-type specific gene regulatory network

Usage

inferCSN(
  object,
  penalty = c("L0", "L0L1", "L0L2"),
  cross_validation = FALSE,
  seed = 1,
  n_folds = 5,
  subsampling_method = c("sample", "meta_cells", "pseudobulk"),
  subsampling_ratio = 1,
  r_squared_threshold = 0,
  regulators = NULL,
  targets = NULL,
  cores = 1,
  verbose = TRUE,
  ...
)

# S4 method for class 'matrix'
inferCSN(
  object,
  penalty = c("L0", "L0L1", "L0L2"),
  cross_validation = FALSE,
  seed = 1,
  n_folds = 5,
  subsampling_method = c("sample", "meta_cells", "pseudobulk"),
  subsampling_ratio = 1,
  r_squared_threshold = 0,
  regulators = NULL,
  targets = NULL,
  cores = 1,
  verbose = TRUE,
  ...
)

# S4 method for class 'sparseMatrix'
inferCSN(
  object,
  penalty = c("L0", "L0L1", "L0L2"),
  cross_validation = FALSE,
  seed = 1,
  n_folds = 5,
  subsampling_method = c("sample", "meta_cells", "pseudobulk"),
  subsampling_ratio = 1,
  r_squared_threshold = 0,
  regulators = NULL,
  targets = NULL,
  cores = 1,
  verbose = TRUE,
  ...
)

Arguments

object

The input data for inferring network.

penalty

The type of regularization. This can take either one of the following choices: "L0", "L0L1", and "L0L2". For high-dimensional and sparse data, "L0L2" is more effective. Default is "L0".

cross_validation

Whether to use cross-validation. Default is FALSE.

seed

The random seed for cross-validation. Default is 1.

n_folds

The number of folds for cross-validation. Default is 5.

subsampling_method

The method to use for subsampling. Options are "sample", "pseudobulk" or "meta_cells".

subsampling_ratio

The percent of all samples used for fit_srm. Default is 1.

r_squared_threshold

Threshold of R² coefficient. Default is 0.

regulators

The regulator genes for which to infer the regulatory network.

targets

The target genes for which to infer the regulatory network.

cores

The number of cores to use for parallelization with foreach::foreach. Default is 1.

verbose

Whether to print progress messages. Default is TRUE.

...

Parameters for other methods.

Value

A data table of regulator-target regulatory relationships, which has three columns: regulator, target, and weight.

Examples

data(example_matrix)
network_table <- inferCSN(
  example_matrix
)
#> ℹ [2026-01-23 02:15:56] Inferring network for <matrix/array>...
#> ◌ [2026-01-23 02:15:56] Checking parameters...
#> ℹ [2026-01-23 02:15:56] Using L0 sparse regression model
#> ℹ [2026-01-23 02:15:56] Using 1 core
#> ⠙ [2026-01-23 02:15:56] Running for g1 [1/18] ■■■                              …
#> ✔ [2026-01-23 02:15:56] Completed 18 tasks in 174ms
#> 
#> ℹ [2026-01-23 02:15:56] Building results
#> ✔ [2026-01-23 02:15:56] Inferring network done
#> ℹ [2026-01-23 02:15:56] Network information:
#> ℹ                         Edges Regulators Targets
#> ℹ                       1   306         18      18

head(network_table)
#>   regulator target     weight
#> 1       g18     g1 -0.9223177
#> 2       g17    g18  0.8770468
#> 3        g4     g3  0.8103230
#> 4       g16    g15  0.7659245
#> 5       g17    g16  0.7558764
#> 6       g12    g11  0.7444053

inferCSN(
  example_matrix,
  regulators = c("g1", "g2"),
  targets = c("g3", "g4")
)
#> ℹ [2026-01-23 02:15:56] Inferring network for <matrix/array>...
#> ◌ [2026-01-23 02:15:56] Checking parameters...
#> ℹ [2026-01-23 02:15:56] Using L0 sparse regression model
#> ℹ [2026-01-23 02:15:56] Using 2 regulators
#> ℹ [2026-01-23 02:15:56] Using 2 targets
#> ℹ [2026-01-23 02:15:56] Using 1 core
#> ℹ [2026-01-23 02:15:56] Building results
#> ✔ [2026-01-23 02:15:56] Inferring network done
#> ℹ [2026-01-23 02:15:56] Network information:
#> ℹ                         Edges Regulators Targets
#> ℹ                       1     4          2       2
#>   regulator target     weight
#> 1        g2     g3  0.9848781
#> 2        g2     g4  0.9230387
#> 3        g1     g4 -0.3847071
#> 4        g1     g3 -0.1732490

inferCSN(
  example_matrix,
  regulators = c("g1", "g2"),
  targets = c("g3", "g0")
)
#> ℹ [2026-01-23 02:15:56] Inferring network for <matrix/array>...
#> ◌ [2026-01-23 02:15:56] Checking parameters...
#> ℹ [2026-01-23 02:15:56] Using L0 sparse regression model
#> ℹ [2026-01-23 02:15:56] Using 2 regulators
#> ! [2026-01-23 02:15:56] The following target is not in the matrix: "g0"
#> ! [2026-01-23 02:15:57] 1 out of 2 targets are in the input matrix
#> ℹ [2026-01-23 02:15:57] Using 1 core
#> ℹ [2026-01-23 02:15:57] Building results
#> ✔ [2026-01-23 02:15:57] Inferring network done
#> ℹ [2026-01-23 02:15:57] Network information:
#> ℹ                         Edges Regulators Targets
#> ℹ                       1     2          2       1
#>   regulator target     weight
#> 1        g2     g3  0.9848781
#> 2        g1     g3 -0.1732490