Optimal Subset Matching without Balance Constraints

Share:

Description

Conducts optimal subset matching as described in the reference.

Usage

1
2
pairmatchelastic(mdist, n = 0, val = 0)
elastic(mdist, n = 0, val = 0)

Arguments

mdist

distance matrix with rows corresponding to treated units and columns corresponding to controls.

n

maximum number of treated units that can be excluded.

val

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).

Details

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 reference.

Value

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.

Author(s)

Paul R. Rosenbaum (original form), Luke Keele, Penn State University, ljk20@psu.edu

Sam Pimentel, University of Pennsylvania, spi@wharton.upenn.edu

References

Rosenbaum, Paul R. (2012) "Optimal Matching of an Optimally Chosen Subset in Observational Studies." Journal of Computational and Graphical Statistics, 21.1, 57-71.