Conducts optimal subset matching as described in the reference.
elastic(mdist, n = 0, val = 0)pairmatchelastic(mdist, n = 0, val = 0)
elastic returns an augmented version of the input matrix
mdist. pairmatchelastic returns a matrix of 1 column whose
values are the column numbers of matched controls and whose rownames are
the row numbers of matched treated units.
distance matrix with rows corresponding to treated units and columns corresponding to controls.
maximum number of treated units that can be excluded.
cost of excluding a treated unit (i.e. we prefer to exclude a
treated unit if it increases the total matched distance by more than
val).
Paul R. Rosenbaum (original forms), modifications by Luke Keele and Sam Pimentel
pairmatchelastic is the main function, which conducts an entire match.
elastic is a helper function which augments the original distance
matrix as described in the reference.
The original versions of these functions were written by Paul Rosenbaum and distributed in the supplemental material to the paper: "Optimal Matching of an Optimally Chosen Subset in Observational Studies," Paul R. Rosenbaum, Journal of Computational and Graphical Statistics, Vol. 21, Iss. 1, 2012.
Rosenbaum, Paul R. (2012) "Optimal Matching of an Optimally Chosen Subset in Observational Studies." Journal of Computational and Graphical Statistics, 21.1, 57-71.