residualspaper {spatstat.data}

Data and Code From JRSS Discussion Paper on Residuals

Description

This dataset contains the point patterns used as examples in the paper of Baddeley et al (2005). [Figure 2 is already available in spatstat.data as the copper dataset.]

R code is also provided to reproduce all the Figures displayed in Baddeley et al (2005). The component plotfig is a function, which can be called with a numeric or character argument specifying the Figure or Figures that should be plotted. See the Examples.

Usage

data(residualspaper)

Format

residualspaper is a list with the following components:

Fig1

The locations of Japanese pine seedlings and saplings from Figure 1 of the paper. A point pattern (object of class "ppp").

Fig3

The Chorley-Ribble data from Figure 3 of the paper. A list with three components, lung, larynx and incin. Each is a matrix with 2 columns giving the coordinates of the lung cancer cases, larynx cancer cases, and the incinerator, respectively. Coordinates are Eastings and Northings in km.

Fig4a

The synthetic dataset in Figure 4 (a) of the paper.

Fig4b

The synthetic dataset in Figure 4 (b) of the paper.

Fig4c

The synthetic dataset in Figure 4 (c) of the paper.

Fig11

The covariate displayed in Figure 11. A pixel image (object of class "im") whose pixel values are distances to the nearest line segment in the copper data.

plotfig

A function which will compute and plot any of the Figures from the paper. The argument of plotfig is either a numeric vector or a character vector, specifying the Figure or Figures to be plotted. See the Examples.

Source

Figure 1: Prof M. Numata. Data kindly supplied by Professor Y. Ogata with kind permission of Prof M. Tanemura.

Figure 3: Professor P.J. Diggle (rescaled by Adrian Baddeley [e-mail])

Figure 4 (a,b,c): Adrian Baddeley [e-mail]

References

Baddeley, A., Turner, R., Moller, J. and Hazelton, M. (2005) Residual analysis for spatial point processes. Journal of the Royal Statistical Society, Series B 67, 617–666.

Examples

if(FALSE) {
  data(residualspaper)
  
  if(require(spatstat.model)) {

  X <- residualspaper$Fig4a
  summary(X)
  plot(X)

  # reproduce all Figures
  residualspaper$plotfig()

  # reproduce Figures 1 to 10
  residualspaper$plotfig(1:10)

  # reproduce Figure 7 (a)
  residualspaper$plotfig("7a")
  }
}