psens {rbounds}

Rosenbaum Bounds for Sign Rank

Description

Function to calculate Rosenbaum bounds for continuous or ordinal outcomes based on Wilcoxon sign rank test.

Usage

# Default Method 
psens(x, y, Gamma = 6, GammaInc = 1)

Arguments

x

Treatment group outcomes in same order as treatment group.

y

Control group outcomes in same order as treatment group.

Gamma

Upper-bound on gamma parameter.

GammaInc

To set user-specified increments for gamma parameter.

Author(s)

Luke Keele, University of Pennsylvania, [e-mail]

References

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

See Also

See also binarysens, hlsens, mcontrol

Examples

 
# Replication of Rosenbaum Sensitivity Tests From Chapter 4 of
# Observational Studies

# Data:  Matched Data of Lead Blood Levels in Children
trt <- c(38, 23, 41, 18, 37, 36, 23, 62, 31, 34, 24, 14, 21, 17, 16, 20,
         15, 10, 45, 39, 22, 35, 49, 48, 44, 35, 43, 39, 34, 13, 73, 25,
         27)
ctrl <- c(16, 18, 18, 24, 19, 11, 10, 15, 16, 18, 18, 13, 19, 10, 16,
          16, 24, 13, 9, 14, 21, 19, 7, 18, 19, 12, 11, 22, 25, 16, 13,
          11, 13)

psens(trt, ctrl)