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