nestsplit

0th

Percentile

Nested Split

Applies two splitting operations to a point pattern, producing a list of lists of patterns.

Keywords
manip, spatial
Usage
nestsplit(X, …)
Arguments
X

Point pattern to be split. Object of class "ppp".

Data determining the splitting factors or splitting regions. See Details.

Details

This function splits the point pattern X into several sub-patterns using split.ppp, then splits each of the sub-patterns into sub-sub-patterns using split.ppp again. The result is a hyperframe containing the sub-sub-patterns and two factors indicating the grouping.

The arguments determine the two splitting factors or splitting regions. Each argument may be:

  • a factor (of length equal to the number of points in X)

  • the name of a column of marks of X (provided this column contains factor values)

  • a tessellation (class "tess")

  • a pixel image (class "im") with factor values

  • a window (class "owin")

  • identified by name (in the form name=value) as one of the formal arguments of quadrats or tess

The arguments will be processed to yield a list of two splitting factors/tessellations. The splits will be applied to X consecutively to produce the sub-sub-patterns.

Value

A hyperframe with three columns. The first column contains the sub-sub-patterns. The second and third columns are factors which identify the grouping according to the two splitting factors.

See Also

split.ppp, quantess

Aliases
  • nestsplit
Examples
# NOT RUN {
  # factor and tessellation
  Nft <- nestsplit(amacrine, marks(amacrine), quadrats(amacrine, 3, 1))
  Ntf <- nestsplit(amacrine, quadrats(amacrine, 3, 1), marks(amacrine))
  Ntf

  # two factors
  big <- with(marks(betacells), area > 300)
  Nff <- nestsplit(betacells, "type", factor(big))

  # two tessellations
  Tx <- quantess(redwood, "x", 4)
  Td <- dirichlet(runifpoint(5, Window(redwood)))
  Ntt <- nestsplit(redwood, Td, Tx)
  Ntt2 <- nestsplit(redwood, Td, ny=3)
# }
Documentation reproduced from package spatstat, version 1.56-1, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.