| by.ppp {spatstat} | R Documentation |
Splits a point pattern into sub-patterns, and applies the function to each sub-pattern.
## S3 method for class 'ppp': by(data, INDICES=marks(data), FUN, ...)
data |
Point pattern (object of class "ppp"). |
INDICES |
Grouping variable. Either a factor, a pixel image with factor values, or a tessellation. |
FUN |
Function to be applied to subsets of data. |
... |
Additional arguments to FUN. |
This is a method for the generic function by
for point patterns (class "ppp").
The point pattern data is first divided into subsets
according to INDICES. Then the function FUN
is applied to each subset. The results of each computation are
returned in a list.
The argument INDICES may be
data.
The levels of INDICES
determine the destination of each point in data.
The ith point of data will be placed in the sub-pattern
split.ppp(data)$l where l = f[i].
"im") with factor values.
The pixel value of INDICES
at each point of data will be used as the classifying variable.
"tess").
Each point of data will be classified according to
the tile of the tessellation into which it falls.
If INDICES is missing, then data must be a multitype point pattern
(a marked point pattern whose marks vector is a factor).
Then the effect is that the points of each type
are separated into different point patterns.
A list (also of class "listof") containing the results returned
from FUN for each of the subpatterns.
Adrian Baddeley adrian@maths.uwa.edu.au http://www.maths.uwa.edu.au/~adrian/ and Rolf Turner r.turner@auckland.ac.nz
ppp,
split.ppp,
cut.ppp,
tess,
im.
# multitype point pattern, broken down by type
data(amacrine)
by(amacrine, FUN=density)
by(amacrine, FUN=function(x) { min(nndist(x)) } )
# passing arguments to FUN
by(amacrine, FUN=clarkevans, correction=c("Donnelly","cdf"))
# point pattern broken down by tessellation
data(swedishpines)
tes <- tess(xgrid=seq(0,96,length=5),ygrid=seq(0,100,length=5))
B <- by(swedishpines, tes, clarkevans, correction="Donnelly")
unlist(lapply(B, as.numeric))