vdraw_intensity_step_regular_R {nhppp}

Vectorized sampling from a non homogeneous Poisson Point Process (NHPPP) from an interval (thinning method) with piecewise constant_majorizers (R)

Description

Vectorized sampling from a non homogeneous Poisson Point Process (NHPPP) from an interval (thinning method) with piecewise constant_majorizers. The majorizers are step functions over equal-length time intevals.

Usage

vdraw_intensity_step_regular_R(
  lambda = NULL,
  lambda_args = NULL,
  Lambda_maj_matrix = NULL,
  lambda_maj_matrix = NULL,
  range_t = NULL,
  tol = 10^-6,
  atmost1 = FALSE,
  force_zt_majorizer = FALSE,
  ...
)

Arguments

lambda

(function) a vectorized intensity function, with one or two arguments. The first is time. The optional second is a named list with additional arguments.

lambda_args

(list) optional list of named arguments for lambda()

Lambda_maj_matrix

(matrix) for the majorizeintegrated intensity rates at the end of each interval

lambda_maj_matrix

(matrix) intensity rates, one per interval

range_t

(vector, or matrix) t_min and t_max, possibly vectorized

tol

(scalar, double) tolerance for the number of events

atmost1

boolean, draw at most 1 event time

force_zt_majorizer

boolean, force the use of the zero-truncated majorizer; - used for flexibility in calling from other functions. Keep the default FALSE unless you know what you are doing.

...

(any) other arguments (ignored – used for flexibility in calling from other functions)

Value

a matrix of event times (columns) per draw (rows) NAs are structural empty spots

Examples

Z <- vdraw_intensity_step_regular_R(
  lambda = function(x, lambda_args = NULL) 0.1 * x,
  range_t = c(1, 10),
  lambda_maj_matrix = matrix(rep(1, 5), nrow = 1)
)