single.network {inferCSN}

Construct network for single gene

Description

Construct network for single gene

Usage

single.network(
  matrix,
  regulators,
  target,
  cross_validation = FALSE,
  seed = 1,
  penalty = "L0",
  algorithm = "CD",
  regulators_num = (ncol(matrix) - 1),
  n_folds = 10,
  percent_samples = 1,
  r_threshold = 0,
  verbose = FALSE,
  ...
)

Arguments

matrix

An expression matrix, cells by genes.

regulators

A character vector with the regulators to consider for CSN inference.

target

Target gene.

cross_validation

Check whether cross validation is used.

seed

The seed used in randomly shuffling the data for cross-validation.

penalty

The type of regularization. This can take either one of the following choices: L0 and L0L2. For high-dimensional and sparse data, such as single-cell sequencing data, L0L2 is more effective.

algorithm

The type of algorithm used to minimize the objective function. Currently CD and CDPSI are supported. The CDPSI algorithm may yield better results, but it also increases running time.

regulators_num

The number of non-zore coefficients, this value will affect the final performance. The maximum support size at which to terminate the regularization path. Recommend setting this to a small fraction of min(n,p) (e.g. 0.05 * min(n,p)) as L0 regularization typically selects a small portion of non-zeros.

n_folds

The number of folds for cross-validation.

percent_samples

The percent of all samples used for sparse.regression. Default set to 1.

r_threshold

Threshold of R^2 or correlation coefficient.

verbose

Print detailed information.

...

Parameters for other methods.

Value

The weight data table of sub-network

Examples

data("example_matrix")
single_network <- single.network(
  example_matrix,
  regulators = colnames(example_matrix),
  target = "g1"
)
head(single_network)

single.network(
  example_matrix,
  regulators = "g1",
  target = "g2"
)