Data and Code From JRSS Discussion Paper on Residuals
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.
data(residualspaper)
residualspaper is a list with the following components:
The locations of Japanese pine seedlings and saplings
from Figure 1 of the paper.
A point pattern (object of class "ppp").
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.
The synthetic dataset in Figure 4 (a) of the paper.
The synthetic dataset in Figure 4 (b) of the paper.
The synthetic dataset in Figure 4 (c) of the paper.
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.
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.
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 Adrian.Baddeley@curtin.edu.au)
Figure 4 (a,b,c): Adrian Baddeley Adrian.Baddeley@curtin.edu.au
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.
if(FALSE) {
data(residualspaper)
if(require(spatstat.core)) {
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")
}
}Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.