redwood.
Strauss (1975) divided the sampling region into two subregions I and II
demarcated by a diagonal line across the region. The spatial pattern
appears to be slightly regular in region I and strongly clustered in
region II. The dataset redwoodfull contains the full point pattern
of 195 trees.
The auxiliary information about the subregions is contained in
redwoodfull.extra, which is a list with entries
diag The coordinates of the diagonal boundary
between regions I and II
regionI Region I as a window object
regionII Region II as a window object
regionR Ripley's subrectangle (approximate)
plot Function to plot the full data and auxiliary markings
}
Ripley (1977) extracted a subset of these data, containing 62 points,
lying within a square subregion which overlaps regions I and II.
He rescaled the data to the unit square.
This has been re-analysed many times, and is the dataset usually known as
``the redwood data'' in the spatial statistics literature.
The exact dataset used by Ripley is supplied in the redwood.
There are some minor inconsistencies with
redwood since it originates from a different digitisation.
The approximate position of the square chosen by Ripley
within the redwoodfull pattern
is indicated by the window redwoodfull.extra$regionR.
data(redwoodfull)redwoodfull is an object of class "ppp"
representing the point pattern of tree locations.
The window has been rescaled to the unit square.
See ppp.object for details of the format of a
point pattern object. The dataset redwoodfull.extra is a list with entries
diag coordinates of endpoints of a line,
in format list(x=numeric(2),y=numeric(2))
regionI a window object
regionII a window object
regionR a window object
plot Function with no arguments
}
Ripley, B.D. (1977) Modelling spatial patterns (with discussion). Journal of the Royal Statistical Society, Series B 39, 172--212.
Strauss, D.J. (1975) A model for clustering. Biometrika 63, 467--475.
redwooddata(redwoodfull)
plot(redwoodfull)
redwoodfull.extra$plot()
# extract the pattern in region II
redwoodII <- redwoodfull[, redwoodfull.extra$regionII]Run the code above in your browser using DataLab