Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

rdpp

Simulation of a Determinantal Point Process


Description

Generates simulated realisations from a determinantal point process.

Usage

rdpp(eig, index, basis = "fourierbasis",
       window = boxx(rep(list(0:1), ncol(index))),
       reject_max = 10000, progress = 0, debug = FALSE, ...)

Arguments

eig

vector of values between 0 and 1 specifying the non-zero eigenvalues for the process.

index

data.frame or matrix (or something acceptable to as.matrix) specifying indices of the basis functions.

basis

character string giving the name of the basis.

window

window (of class "owin", "box3" or "boxx") giving the domain of the point process.

reject_max

integer giving the maximal number of trials for rejection sampling.

progress

integer giving the interval for making a progress report. The value zero turns reporting off.

debug

logical value indicating whether debug informationb should be outputted.

...

Ignored.

Value

A point pattern (object of class "ppp").

Author(s)

Adrian Baddeley Adrian.Baddeley@curtin.edu.au, Rolf Turner r.turner@auckland.ac.nz and Ege Rubak rubak@math.aau.dk.

Examples

index <- expand.grid(-2:2,-2:2)
eig <- exp(-rowSums(index^2))
X <- rdpp(eig, index)
X
## To simulate a det. projection p. p. with the given indices set eig=1:
XX <- rdpp(1, index)
XX

spatstat.core

Core Functionality of the 'spatstat' Family

v2.1-2
GPL (>= 2)
Authors
Adrian Baddeley [aut, cre], Rolf Turner [aut], Ege Rubak [aut], Kasper Klitgaard Berthelsen [ctb], Achmad Choiruddin [ctb], Jean-Francois Coeurjolly [ctb], Ottmar Cronie [ctb], Tilman Davies [ctb], Julian Gilbey [ctb], Yongtao Guan [ctb], Ute Hahn [ctb], Kassel Hingee [ctb], Abdollah Jalilian [ctb], Marie-Colette van Lieshout [ctb], Greg McSwiggan [ctb], Tuomas Rajala [ctb], Suman Rakshit [ctb], Dominic Schuhmacher [ctb], Rasmus Plenge Waagepetersen [ctb], Hangsheng Wang [ctb]
Initial release
2021-04-17

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.