copper {spatstat} | R Documentation |
These data come from an intensive geological survey of a 70 x 158 km region in central Queensland, Australia. They consist of 67 points representing copper ore deposits, and 146 line segments representing geological `lineaments'. Lineaments are linear features, visible on a satellite image, that are believed to consist largely of geological faults (Berman, 1986, p. 55). It would be of great interest to predict the occurrence of copper deposits from the lineament pattern, since the latter can easily be observed on satellite images.
These data were introduced and analysed by Berman (1986). They have also been studied by Berman and Turner (1992), Baddeley and Turner (2000) and Foxall and Baddeley (2002).
Many analyses have been performed on the southern half of the data only. This subset is also provided.
data(copper)
copper
is a list with the following entries:
"ppp"
representing the full point pattern of copper deposits.
See ppp.object
for details of the format.
"owin"
.
"ppp"
.
"im"
.
Dr J. Huntington. Coordinates kindly provided by Dr. Mark Berman, CSIRO, Sydney, Australia.
Baddeley, A. and Turner, R. (2000) Practical maximum pseudolikelihood for spatial point patterns. Australian and New Zealand Journal of Statistics 42, 283–322.
Berman, M. (1986). Testing for spatial association between a point process and another stochastic process. Applied Statistics 35, 54–62.
Berman, M. and Turner, T.R. (1992) Approximating point process likelihoods with GLIM. Applied Statistics 41, 31–38.
Foxall, R. and Baddeley, A. (2002) Nonparametric measures of association between a spatial point process and a random set, with geological applications. Applied Statistics 51, 165–182.
data(copper) # Plot full dataset plot(copper$points) cl <- copper$lines segments(cl[,1], cl[,2], cl[,3], cl[,4]) # Plot southern half of data plot(copper$SouthPoints) cl <- copper$SouthLines segments(cl[,1], cl[,2], cl[,3], cl[,4]) ## Not run: Z <- copper$SouthDistance() plot(Z) X <- copper$SouthPoints ppm(X, ~D, covariates=list(D=Z)) ## End(Not run)