hlsens | R Documentation |
Function to calculate Rosenbaum bounds for continuous or ordinal outcomes based on Hodges-Lehmann point estimate.
# Default Method hlsens(x, y, pr = 0.1, Gamma = 6, GammaInc = 1)
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. |
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.
Luke Keele, University of Pennsylvania, luke.keele@gmail.com
Rosenbaum, Paul R. (2002) Observational Studies. Springer-Verlag.
See also binarysens
,
psens
,
mcontrol
# 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.