This function returns a "maximal" set of disjoint pairs of linked markers,
given a genetic marker map and a subset of the markers included in the map.
The pairs are identified in a greedy manner, successively choosing the
closest markers on each chromosome.