Create a Point Matching
Creates an object of class
a matching of two planar point patterns (objects of class
pppmatching(X, Y, am, type = NULL, cutoff = NULL, q = NULL, mdist = NULL)
- Two point patterns (objects of class
Y$nmatrix with entries $\geq 0$ that specifies which points are matched and with what weight; alternatively, an object that can be coerced to this form by
- A character string giving the type of the matching.
NULLfor a generic or unknown matching.
- cutoff, q
- Numerical values specifying the cutoff value $> 0$ for interpoint distances and
the order $q \in [1,\infty]$ of the average that is applied to them.
NULLif not applicable or unknown.
- Numerical value for the distance to be associated with the matching.
am is interpreted as a "generalized adjacency matrix":
[i,j]-th entry is positive, then the
X and the
j-th point of
Y are matched and the
value of the entry gives the corresponding weight of the match. For
an unweighted matching all the weights should be set to $1$.
# a random unweighted complete matching X <- runifpoint(10) Y <- runifpoint(10) am <- r2dtable(1, rep(1,10), rep(1,10))[] # generates a random permutation matrix m <- pppmatching(X, Y, am) summary(m) m$matrix plot(m) # a random weighted complete matching X <- runifpoint(7) Y <- runifpoint(7) am <- r2dtable(1, rep(10,7), rep(10,7))[]/10 # generates a random doubly stochastic matrix m2 <- pppmatching(X, Y, am) summary(m2) m2$matrix # Note: plotting does currently not distinguish # between different weights plot(m2)