Skip to contents

This function is used to calculate the driver activities of one single group from the sparse eset obj and networks generated by SJARACNe. To calculate driver activities of multiple groups from a scMINER directory, please use getActivity_inBatch().

Usage

getActivity_individual(
  input_eset,
  network_file.tf = NULL,
  network_file.sig = NULL,
  driver_type = "TF_SIG",
  activity_method = "mean",
  do.z_normalization = TRUE
)

Arguments

input_eset

The group-specific expression set object which has been filtered, normalized and log-transformed

network_file.tf

The path to the TF network file generated by SJARACNe

network_file.sig

The path to the SIG network file generated by SJARACNe

driver_type

Character, type of the pre-defined driver list to use: "TF" for transcriptional factors only, "SIG" for signaling genes only, or "TF_SIG" for both. Default: "TF_SIG".

activity_method

Character, method used to calculate the activity: "mean" (the default), "weightedmean", "absmean" or "maxmean".

do.z_normalization

Logical, whether to do the z-normalization on the gene expression values in each sample. Set if to FALSE only when the expression values has been scaled in each cell. Default: TRUE.

Value

A expression set object of the group-of-interest. The assayData is the activity matrix of all cells of all groups, drivers by cells. The phenoData and featureData are exactly save with the input eset.

Examples

data(pbmc14k_expression.eset)
activity_group.eset <- getActivity_individual(input_eset = pbmc14k_expression.eset[, pbmc14k_expression.eset$trueLabel == "B"],
                                              network_file.tf = system.file("extdata/demo_pbmc14k/SJARACNe/B/TF/bt100_pc001/consensus_network_ncol_.txt", package = "scMINER"),
                                              network_file.sig = system.file("extdata/demo_pbmc14k/SJARACNe/B/SIG/bt100_pc001/consensus_network_ncol_.txt", package = "scMINER"),
                                              driver_type = "TF_SIG")