rbounds (version 2.1)

data.prep: Format Match() object for sensitivity test.

Description

This function reshapes the output from Match() to create the necessary objects for mcontrol().

Usage

# Default Method 
data.prep(obj, Y = NULL, group.size = 3)

Arguments

obj

An object from the Match() function.

Y

Internal argument to the Match object.

group.size

The size of the matched groups. Three for one treated unit and two control units.

Value

Y

The matched outcomes

id

A vector which identifies the matched groups: 1, 1, 1 for matched group one; 2, 2, 2 for match group 2, etc.

treat

A vector with 1's for treated units and 0's for control units

Details

This functions takes a Match() object and formats it for use with the mcontrol() function. The output is a list with the three objects needed for the arguments of the mcontrol() function.

References

Rosenbaum, Paul R. (2002) Observational Studies. Springer-Verlag.

See Also

See also binarysens, psens, hlsens, Match, mcontrol

Examples

# Load Matching Software and Data
library(Matching)
data(lalonde)

# Estimate Propensity Score
DWglm <- glm(treat~age + I(age^2) + educ + I(educ^2) + black + hisp +
             married + nodegr + re74 + I(re74^2) + re75 + I(re75^2) +
             u74 + u75, family=binomial, data=lalonde)

# Save data objects
Y  <- lalonde$re78    # the outcome of interest
Tr <- lalonde$treat   # the treatment of interest

# Match      
mDW  <- Match(Y = Y, Tr = Tr, X = DWglm$fitted, M = 2)

# One should check balance, but let's skip that step for now.

# Create Data Object
tmp <- data.prep(mDW, group.size=3)

# Sensitivity Test
mcontrol(tmp$Y, tmp$id, tmp$treat, group.size = 3)