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])
Z <- copper$SouthDistance()
plot(Z)
X <- copper$SouthPoints
ppm(X, ~D, covariates=list(D=Z))
Run the code above in your browser using DataLab