## -----------------------------------------------------------------------------
# Load libraries
library("netdist")
library("purrr")
## -----------------------------------------------------------------------------
# Maximum graphlet size to calculate counts and netdis statistic for.
max_graphlet_size <- 4
# Ego network neighbourhood size
neighbourhood_size <- 2
# Minimum size of ego networks to consider
min_ego_nodes <- 3
min_ego_edges <- 1
# Reference graph
ref_path <- system.file(file.path("extdata", "random", "ER_1250_10_1"),
package = "netdist")
ref_graph <- read_simple_graph(ref_path, format = "ncol")
## -----------------------------------------------------------------------------
source_dir <- system.file(file.path("extdata", "VRPINS"), package = "netdist")
graphs <- read_simple_graphs(source_dir, format = "ncol", pattern = "*")
## -----------------------------------------------------------------------------
# Calculate netdis statistics
results <- netdis_many_to_many(graphs,
ref_graph,
max_graphlet_size = max_graphlet_size,
neighbourhood_size = neighbourhood_size,
min_ego_nodes = min_ego_nodes,
min_ego_edges = min_ego_edges)
print(results$netdis)
print(results$comp_spec)
## -----------------------------------------------------------------------------
binning_fn <- purrr::partial(binned_densities_adaptive,
min_counts_per_interval = 10,
num_intervals = 50)
# Calculate netdis statistics
results <- netdis_many_to_many(graphs,
ref_graph,
max_graphlet_size = max_graphlet_size,
neighbourhood_size = neighbourhood_size,
min_ego_nodes = min_ego_nodes,
min_ego_edges = min_ego_edges,
binning_fn = binning_fn)
print(results$netdis)
print(results$comp_spec)
## -----------------------------------------------------------------------------
bin_counts_fn <- density_binned_counts_gp
exp_counts_fn <- purrr::partial(netdis_expected_counts,
scale_fn = NULL)
# Calculate netdis statistics
results <- netdis_many_to_many(graphs,
ref_graph = NULL,
max_graphlet_size = max_graphlet_size,
neighbourhood_size = neighbourhood_size,
min_ego_nodes = min_ego_nodes,
min_ego_edges = min_ego_edges,
bin_counts_fn = bin_counts_fn,
exp_counts_fn = exp_counts_fn)
print(results$netdis)
print(results$comp_spec)
## -----------------------------------------------------------------------------
binning_fn <- single_density_bin
bin_counts_fn <- density_binned_counts
exp_counts_fn <- netdis_expected_counts
# Calculate netdis statistics
results <- netdis_many_to_many(graphs,
ref_graph = NULL,
max_graphlet_size = max_graphlet_size,
neighbourhood_size = neighbourhood_size,
min_ego_nodes = min_ego_nodes,
min_ego_edges = min_ego_edges,
binning_fn = binning_fn,
bin_counts_fn = bin_counts_fn,
exp_counts_fn = exp_counts_fn)
print(results$netdis)
print(results$comp_spec)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.