View source: R/grn_inference.R
grn_filter | R Documentation |
Filter a gene regulatory network based on optimal scale-free topology fit
grn_filter(edgelist, nsplit = 10, bp_param = BiocParallel::SerialParam())
edgelist |
A gene regulatory network represented as an edge list. |
nsplit |
Number of groups in which the edge list will be split. Default: 10. |
bp_param |
BiocParallel back-end to be used. Default: BiocParallel::SerialParam() |
The edge list will be split in n groups and the scale-free topology fit will be tested for each subset of the edge list. For instance, if an edge list of 10000 rows is used as input, the function will test SFT fit for the top 1000 edges, then top 2000 edges, and so on up to the whole edge list.
The edge list that best fits the scale-free topology.
data(filt.se)
tfs <- sample(rownames(filt.se), size=50, replace=FALSE)
grn_list <- grn_combined(filt.se, regulators=tfs, nTrees=2)
ranked_grn <- grn_average_rank(grn_list)
# split in only 2 groups for demonstration purposes
filtered_edges <- grn_filter(ranked_grn, nsplit=2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.