pcf_icpp: PCF for iterated cluster point process

Description Usage Arguments Details Value

Description

Closure function (ie. returns a function) that: Calculats the pcf for the n'th generation iterated Neyman-Scott type process with or without noise. It is asusmed that nothing except for the intensity depends on the generation number. The offspring density is assumed to be Gaussian. The retention probability in 0, meaning that points do not survive to the next generation.

Usage

1
2
3
4
5
6
pcf_icpp(n, init_int, clsiz, c = 1, sd, initial = "Pois",
  kernel_sd_det_init = NULL, kernel_sd_per_init = NULL,
  weight_det_init = 1, weight_per_init = 1/2, noise = "non",
  noise_int = 0, kernel_sd_det_noise = NULL,
  kernel_sd_per_noise = NULL, weight_det_noise = 1,
  weight_per_noise = 1/2)

Arguments

n

Positive integer. The number of iterations.

init_int

Non-negative numeric. Intensity of the initial generation.

clsiz

Non-negative numeric. The expected cluster size.

c

Numeric. Constant related to the over-/underdispersion of number of point in a cluster. Default is 1, corresponding to the number of points in a cluster follow a poisson distriution. See article on iterated cluster point processes.

sd

Non-negative number. Offspring density standard deviation.

initial

Vector of strings or a single string. Decribing the type of the initial process. Can be either "Pois" (default), "det", or "per".

kernel_sd_det_init

Non-negative numeric. Kernel standard deviation. (interaction 'radius' for the initial process). Ignored if initial is "Pois".

kernel_sd_per_init

Non-negative numeric. Kernel standard deviation. (interaction 'radius' for the initial process). Ignored if initial is "Pois".

weight_det_init

A positive integer. The alpha weight for an alpha weighted DPP. Default is 1 corresponding to the DPP (most repulsive gaussian DPP). Ignored if initial is "Pois".

weight_per_init

A positive half integer. The alpha weight for an alpha weighted PPP. Default is 1/2 corresponding to the most clustered gaussian PPP which is a Cox process.. Ignored if initial is "Pois".

noise

Vector of strings or a single string. Decribing the type of the noise process. Can be either "non" (default) "det", "pois", or "per".

noise_int

Non-negative numeric or the string "sr". If "sr" same reproduction system is assumed and noise_int = init_int * (1 - clsiz). Ignored if noise is "non".

kernel_sd_det_noise

Ignored if noise is "non" or Pois".

kernel_sd_per_noise

Ignored if noise is "non" or Pois".

weight_det_noise

A positive integer. The alpha weight for an alpha weighted DPP. Default is 1 corresponding to the DPP (most repulsive gaussian DPP). Ignored if noise is "non" or Pois".

weight_per_noise

A positive half integer. The alpha weight for an alpha weighted PPP. Default is 1/2 corresponding to the most clustered gaussian PPP which is a Cox process. Ignored if noise is "non" or Pois".

Details

The function is a closure function. The benefit of making a closure function here is that I can run the output function multiple times without running initial checks and preliminary calculations. We want this here because we have to compile the code for maybe 1000 values of r.

The initial processes and noise processes may be either Poisson, weighted determinantal or weighted permanetal point process.

Value

A function that takes inter-point distance (commonly referred to as r) as the only argument.


adchSTATS/icpp documentation built on May 28, 2019, 3:54 p.m.