Examples in other dimensions

This vignette contains a number of examples on how to use clugenr in 1D and 5D. Examples require the following setup code:

# Load the clugenr library
library(clugenr)

# Load functions for plotting examples
source("plot_examples_1d.R", local = knitr::knit_global())
source("plot_examples_nd.R", local = knitr::knit_global())

# Keep examples reproducible in newer R versions
RNGversion("3.6.0")

Basic 1D example with density plot

The 1D examples were plotted with the plot_examples_1d() function available here.

seed <- 222
# Custom proj_dist_fn: point projections placed using the Weibull distribution
proj_wbull <- function(len, n) rweibull(n, shape = 1.5, scale = len / 2) - len / 2
e082 <- clugen(1, 3, 2000, 1, 0, 10, 6, 1.5, 0, seed = seed)
e083 <- clugen(1, 3, 2000, 1, 0, 10, 6, 1.5, 0, seed = seed, proj_dist_fn = "unif")
e084 <- clugen(1, 3, 2000, 1, 0, 10, 6, 1.5, 0, seed = seed, proj_dist_fn = proj_wbull)
plot_examples_1d(list(e = e082, t = "e082: proj_dist_fn='norm'"),
                 list(e = e083, t = "e083: proj_dist_fn='unif'"),
                 list(e = e084, t = "e084: custom proj_dist_fn (Weibull)"),
                 pmargin = 0.07,
                 ymax = 0.575)

5D example with default optional arguments

The 5D examples were plotted with the plot_examples_nd() function available here.

nd <- 5
seed <- 321
e085 <- clugen(nd, 6, 1500, c(1, 1, 0.5, 0, 0), pi / 16, rep.int(30, nd), 30, 4, 3,
              seed = seed)
plot_examples_nd(e085, "e085: 5D with optional parameters set to defaults")

5D example with proj_dist_fn = "unif" and point_dist_fn = "n"

nd <- 5
seed <- 123
e086 <- clugen(nd, 6, 1500, c(0.1, 0.3, 0.5, 0.3, 0.1), pi / 12, rep.int(30, nd),
              35, 5, 3.5, seed = seed,
              proj_dist_fn = "unif", point_dist_fn = "n")
plot_examples_nd(e086, "e086: 5D with proj_dist_fn='unif' and point_dist_fn='n'")