SDALGCPCreatePoint: SDALGCPCreatePoint function

In olatunjijohnson/SDALGCP: Spatially Discrete Approximation to Log-Gaussian Cox Processes for Aggregated Disease Count Data

Description

This wrapper function generate a random point pattern using Simple Sequential Inhibition (SSI) process, uniform sampling and regular grid point.

Usage

SDALGCPCreatePoint(
  my_shp,
  delta,
  weighted = FALSE,
  lambdamax = NULL,
  pop = NULL,
  pop_shp = NULL,
  n = NULL,
  method = 1,
  plot = FALSE,
  rho = NULL,
  giveup = NULL
)

Arguments

my_shp	A SpatialPolygons or SpatialPolygonsDataFrame object containing the polygons (i.e each regions).
delta	distance between points
weighted	To specify if you want to use the population density, default to FALSE, i.e population density is not used.
lambdamax	the maximum value of the population density in the polygon.
pop	the population density.
pop_shp	Optional, The raster of population density map for population weighted approach.
n	optional; the number of points to create in the polygon, if not supplied, it is computed as n = rho*|A|*4/(π*delta^2)
method	To specify which method to use to sample the points, the options are 1 for Simple Sequential Inhibition (SSI) process, 2 for Uniform sampling and 3 for regular grid. 1 is the default
plot	To display the plot of the points inside the polygon, default to TRUE.
rho	Optional, The packing density, default set to 0.55
giveup	Number of rejected proposals after which the algorithm should terminate.

Details

This algorithm generates points inside the polygon using three algorithms specified in the method.

Value

It returns a list of the coordinates of the points created in each polygon and it has an associated attribute weighted which is either TRUE or FALSE to indicate if the population density is used or not.

Author(s)

Olatunji O. Johnson [e-mail]

Emanuele Giorgi [e-mail]

Peter J. Diggle [e-mail]

See Also

SDALGCPSSIPoint, SDALGCPUniformPoint, SDALGCPRegularPoint

olatunjijohnson/SDALGCP documentation built on March 20, 2021, 4:24 a.m.