Nested Split

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

manip, spatial
nestsplit(X, ...)
Point pattern to be split. Object of class "ppp".
Data determining the splitting factors or splitting regions. See 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 inX)
  • the name of a column of marks ofX(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 formname=value) as one of the formal arguments ofquadratsortess
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.


  • 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

  • nestsplit
# factor and tessellation
  Nft <- nestsplit(amacrine, marks(amacrine), quadrats(amacrine, 3, 1))
  Ntf <- nestsplit(amacrine, quadrats(amacrine, 3, 1), marks(amacrine))

  # 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.41-1, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.