acceptreject {bSims} R Documentation

## Spatial point process simulator

### Description

Spatial point process simulator based on accept/reject algorithm.

### Usage

```acceptreject(n, f = NULL, x0 = 0, x1 = 1, y0 = 0, y1 = 1,
m = 0, maxit = 100, fail = FALSE)
```

### Arguments

 `n` number of points to generate. `f` a function returning probability (value between 0 and 1) given distance as the first and only argument. The function generates spatially uniform Poisson point process (complete spatial randomness) when `NULL`. `x0, x1, y0, y1` x and y ranges (bounding box). `m` margin width for avoiding edge effects. `maxit` maximum number of iterations per point to try if no acceptance happens. `fail` logical, what to do when there is a problem. `TRUE` gives error, the default `FALSE` gives only a warning.

### Value

A matrix with `n` rows and 2 columns for x and y coordinates.

Peter Solymos

### Examples

```## complete spatial randomness
plot(acceptreject(100), asp=1)

## more systematic
distance <- seq(0,1,0.01)
f <- function(d)
(1-exp(-d^2/0.1^2) + dlnorm(d, 0.2)/dlnorm(exp(0.2-1),0.2)) / 2
op <- par(mfrow = c(1, 2))
plot(distance, f(distance), type="l")
plot(acceptreject(100, f, m=1), asp=1)
par(op)

## more clustered
f <- function(d)
exp(-d^2/0.1^2) + 0.5*(1-exp(-d^2/0.4^2))
op <- par(mfrow = c(1, 2))
plot(distance, f(distance), type="l")
plot(acceptreject(100, f, m=1), asp=1)
par(op)
```

[Package bSims version 0.2-1 Index]