## S3 method for class 'pp3':
closepairs(X, rmax, ordered=TRUE,
what=c("all", "indices"), ...)## S3 method for class 'pp3':
crosspairs(X, Y, rmax, what=c("all", "indices"), ...)
"pp3"
).ordered=TRUE
, each pair will appear twice
in the output, as (i,j)
and again as (j,i)
. If
ordered=FALSE
what="all"
(the default) then the
returned information includes the indices i,j
of each pair,
their x,y,z
coordinates, and the di
and j
,
and possibly other components as described under Details.rmax
. Floating-point numbers in a computer
are not mathematical Real Numbers: they are approximations using
finite-precision binary arithmetic.
The approximation is accurate to a tolerance of about
.Machine$double.eps
.
If the true interpoint distance $d$ and the threshold rmax
are equal, or if their difference is no more than .Machine$double.eps
,
the result may be incorrect.
closepairs(X,rmax)
identifies all pairs of neighbours
in the pattern X
and returns them. The result is
a list with the following components:
[object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object]
If what="indices"
then only the components i
and
j
are returned. This is slightly faster.
crosspairs(X,rmax)
identifies all pairs of neighbours
(X[i], Y[j])
between the patterns X
and Y
,
and returns them. The result is
a list with the same format as for closepairs
.
closepairs
X <- pp3(runif(10), runif(10), runif(10), box3(c(0,1)))
Y <- pp3(runif(10), runif(10), runif(10), box3(c(0,1)))
a <- closepairs(X, 0.1)
b <- crosspairs(X, Y, 0.1)
Run the code above in your browser using DataLab