Difference of Disc Areas
Computes the area of that part of a disc that is not covered by other discs.
areaLoss(X, r, ngrid=spatstat.options("ngrid.disc"), subset=NULL)
- Locations of the centres of discs.
A point pattern (object of class
- Disc radius, or vector of disc radii.
- Dimensions of the
ngrid * ngridarray of test points that will be used to estimate the uncovered area.
- Optional. Index identifying a subset of the points of
Xfor which the area difference should be computed.
This function computes, for each point
and for each radius
the area of that part of the disc of radius
r centred at the
X[i] that is not covered by any of the
other discs of radius
r centred at the points
j not equal to
This area is important in some calculations related to
the area-interaction model
The result is a matrix, with one row for each point in
and one column for each entry of
The area is calculated by discrete approximation, using an
ngrid * ngrid grid of test points.
- A matrix with one row for each point in
X[subset]) and one column for each value in
data(cells) areaLoss(cells, 0.1, 32)