gammalist: A list of objects produced by either gammaKpar or gammaCKpar.
nobs.a: number of observations in dataset 1
nobs.b: number of observations in dataset 2
em: parameters obtained from the Expectation-Maximization algorithm under the MAR assumption. These estimates are produced by emlinkMARmov
thresh: is the interval of posterior zeta values for the agreements that we want to examine closer. Ranges between 0 and 1. Can be a vector of length 1 (from specified value to 1) or 2 (from first specified value to second specified value).
n.cores: Number of cores to parallelize over. Default is NULL.
Returns
matchesLink returns an nmatches X 2 matrix with the indices of the matches rows in dataset A and dataset B.