Visualize the GRN with regulators as green square nodes, and targets as gray discs.
draw_network(nodes, edges)dataframe containing nodes, from DIANE::network_data()$nodes
dataframe containing edges, from DIANE::network_data()$edges
if (FALSE) {
data("abiotic_stresses")
data("regulators_per_organism")
genes <- get_locus(abiotic_stresses$heat_DEGs)
# mat was inferred using the function network_inference
mat <- abiotic_stresses$heat_DEGs_regulatory_links
network <- DIANE::network_thresholding(mat, n_edges = length(genes))
data <- network_data(network, regulators_per_organism[["Arabidopsis thaliana"]])
DIANE::draw_network(data$nodes, data$edges)}