knitr::opts_chunk$set(echo = TRUE, cache = TRUE, cache.rebuild = TRUE, autodep = TRUE)
set.seed(1) library(StartNetwork) library(parallel) library(ggplot2) n = 5 replicates = 100 datapoints = 10
mech_net_triangles <- purrr::partial(mech_net_triangles_n, n = !!n) true_value <- 0.5 theta_p <- rep(seq(0.4, 0.6, by = 0.02), datapoints) theta_s <- log(true_value/(1 - true_value)) cl <- parallel::makeCluster(parallel::detectCores()) g <- parallel::parLapply(cl, theta_p, StartNetwork::KL_ss, theta_s = theta_s, replicates = replicates, sorted = TRUE, mech_net = mech_net_triangles, lstat = function(x){length(igraph::triangles(x))/3}, mirror = TRUE, type = c("Bianconi", "Liebenau")) parallel::stopCluster(cl)
g_tidy <- StartNetwork::tidy_g(g) ggplot(g_tidy) + aes(x = parameter, y = value, group = interaction(parameter, type), col = type) + geom_boxplot(outlier.shape = NULL, position = "identity") + facet_wrap( ~ key, scales = "free_y") + ggplot2::geom_vline(mapping = aes(xintercept = true_value))
mech_net_triangles <- purrr::partial(mech_net_triangles_n, n = !!n) true_value <- 0.4 theta_p <- rep(seq(0.3, 0.5, by = 0.02), datapoints) theta_s <- log(true_value/(1 - true_value)) cl <- parallel::makeCluster(parallel::detectCores()) g <- parallel::parLapply(cl, theta_p, StartNetwork::KL_ss, theta_s = theta_s, replicates = replicates, sorted = TRUE, mech_net = mech_net_triangles, lstat = function(x){length(igraph::triangles(x))/3}, mirror = TRUE, type = c("Bianconi", "Liebenau")) parallel::stopCluster(cl)
g_tidy <- StartNetwork::tidy_g(g) ggplot(g_tidy) + aes(x = parameter, y = value, group = interaction(parameter, type), col = type) + geom_boxplot(outlier.shape = NULL, position = "identity") + facet_wrap( ~ key, scales = "free_y") + ggplot2::geom_vline(mapping = aes(xintercept = true_value))
mech_net_triangles <- purrr::partial(mech_net_triangles_n, n = !!n) true_value <- 0.3 theta_p <- rep(seq(0.04, 0.4, by = 0.02), datapoints) theta_s <- log(true_value/(1 - true_value)) cl <- parallel::makeCluster(parallel::detectCores()) g <- parallel::parLapply(cl, theta_p, StartNetwork::KL_ss, theta_s = theta_s, replicates = replicates, sorted = TRUE, mech_net = mech_net_triangles, lstat = function(x){length(igraph::triangles(x))/3}, mirror = TRUE, type = c("Bianconi", "Liebenau")) parallel::stopCluster(cl)
g_tidy <- StartNetwork::tidy_g(g) ggplot(g_tidy) + aes(x = parameter, y = value, group = interaction(parameter, type), col = type) + geom_boxplot(outlier.shape = NULL, position = "identity") + facet_wrap( ~ key, scales = "free_y") + ggplot2::geom_vline(mapping = aes(xintercept = true_value))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.