second_GBM_step: Re-iterate through the core GBM model building with optimal...

Description Usage Arguments Value Author(s) See Also

View source: R/proposed_steps.R


This function re-performs the core GBM model building (only one time) using the optimal set of transcription factors obtained from select_ideal_k followed by get_colids for individual target gene to return a regularized GRN.


second_GBM_step(E, K, df_colids, tfs, targets, Ntfs, Ntargets, lf,  M,  nu, s_f)



N-by-p expression matrix. Columns correspond to genes, rows correspond to experiments. E is expected to be already normalized using standard methods, for example RMA. Colnames of E is the set of all genes.


N-by-p initial perturbation matrix. It directly corresponds to E matrix, e.g. if K[i,j] is equal to 1, it means that gene j was knocked-out in experiment i. Single gene knock-out experiments are rows of K with only one value 1. Colnames of K is set to be the set of all genes. By default it's a matrix of zeros of the same size as E, e.g. unknown initial perturbation state of genes.


A matrix made up of column vectors where each column vector represents the optimal set of active Tfs which regulate each target gene and obtained from get_colids. Some column vectors are just made up of zeros indicating that corresponding target genes are isolated and not regulated by any Tf


List of names of transcription factors.


List of names of target genes.


Total number of transcription factors used in the experiment.


Total number of target genes used in the experiment


Loss Function: 1 -> Least Squares and 2 -> Least Absolute Deviation


Number of extensions in boosting model, e.g. number of iterations of the main loop of RGBM algorithm. By default it's 5000.


Shrinkage factor, learning rate, 0<nu<=1. Each extension to boosting model will be multiplied by the learning rate. By default it's 0.001.


Sampling rate of transcription factors, 0<s_f<=1. Fraction of transcription factors from E, as indicated by tfs vector, which will be sampled without replacement to calculate each extesion in boosting model. By default it's 0.3.


Returns a regularized GRN of the form Ntfs-by-Ntargets


Raghvendra Mall <>

See Also


RGBM documentation built on May 1, 2019, 7:59 p.m.

Related to second_GBM_step in RGBM...