spatstat (version 1.46-1)

residualspaper: 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 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) Figure 4 (a,b,c): \adrian

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

## Not run: 
#   data(residualspaper)
#   
#   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")
# ## End(Not run)