hlsens: Rosenbaum Bounds for Hodges-Lehmann Point Estimate In rbounds: Perform Rosenbaum Bounds Sensitivity Tests for Matched and Unmatched Data

 hlsens R Documentation

Rosenbaum Bounds for Hodges-Lehmann Point Estimate

Description

Function to calculate Rosenbaum bounds for continuous or ordinal outcomes based on Hodges-Lehmann point estimate.

Usage

```# Default Method
hlsens(x, y, pr = 0.1, 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. `pr` Search precision parameter. `Gamma` Upper-bound on gamma parameter. `GammaInc` To set user specified increments for gamma parameter.

Details

For large data sets this function can be quite slow if pr is set to low. If the data set is larger, it is best to set pr to .5 before trying values such as .01. Generally, the results from the function are insensitive to the value for pr.

Author(s)

Luke Keele, University of Pennsylvania, luke.keele@gmail.com

References

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

See also `binarysens`, `psens`, `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)

hlsens(trt, ctrl)
```

rbounds documentation built on April 30, 2022, 1:07 a.m.