# amplify: Amplification of sensitivity analysis in observational... In sensitivitymv: Sensitivity Analysis in Observational Studies

## Description

Use of this function is illustrated in Rosenbaum (2015). Uses the method in Rosenbaum and Silber (2009) to interpret a value of the sensitivity parameter gamma, for instance the parameter in the senmv function. Each value of gamma amplifies to a curve (lambda,delta) in a two-dimensional sensitivity analysis, the inference being the same for all points on the curve. That is, a one-dimensional sensitivity analysis in terms of gamma has a two-dimensional interpretation.

## Usage

 `1` ```amplify(gamma, lambda) ```

## Arguments

 `gamma` gamma > 1 is the value of the sensitivity parameter, for instance the parameter in senmv. length(gamma)>1 will generate an error. `lambda` lambda is a vector of values > gamma. An error will result unless lambda[i] > gamma > 1 for every i.

## Details

A single value of gamma, say gamma = 3.5 in the example, corresponds to a curve of values of (lambda, delta), including (4, 26), (6,8), (8,6), and (11,5) in the example. An unobserved covariate that is associated with a lambda = 6 fold increase in the odds of treatment and a delta = 8 fold increase in the odds of a positive pair difference is equivalent to gamma = 3.5.

The curve is gamma = (lambda*delta + 1)/(lambda+delta). Amplify is given one gamma and a vector of lambdas and solves for the vector of deltas. The calculation is elementary.

This interpretation of gamma is developed in detail in Rosenbaum and Silber (2009), and it makes use of Wolfe's (1974) family of semiparametric deformations of an arbitrary symmetric distribuiton. Strictly speaking, the amplification describes matched pairs, not matched sets. The senmv function views a k-to-1 matched set with k controls matched to one treated individual as a collection of k correlated treated-minus-control matched pair differences; see Rosenbaum (2007). For matched sets, it is natural to think of the amplification as describing any one of the k matched pair differences in a k-to-1 matched set. The curve has asymptotes that the function amplify does not compute: gamma corresponds with (lambda,delta) = (gamma, Inf) and (Inf, gamma). A related though distict idea is developed in Gastwirth et al (1998). The two approaches agree when the outcome is binary, that is, for McNemar's test.

## Value

Returns a vector of values of delta of length(lambda) with names lambda.

## Note

The example expands the discussion of Table 1 in Rosenbaum (2007). The study is insensitive to a bias of gamma = 3.5. An unobserved covariate associated with a lambda = 6 fold increase in the odds of treatment and a delta= 8 fold increase in the odds of positive pair difference is equivalent to gamma = 3.5. Also, gamma = 3.5 is equivalent to (lambda,delta) = (4,26), (8,6) and (11,5).

## Author(s)

Paul R. Rosenbaum

## References

Gastwirth, J. L., Krieger, A. M., Rosenbaum, P. R. (1998) Dual and simultaneous sensitivity analysis for matched pairs. Biometrika, 85, 907-920.

Rosenbaum, P. R. (2007) Sensitivity analysis for m-estimates, tests and confidence intervals in matched observational studies. Biometrics, 2007, 63, 456-464.

Rosenbaum, P. R. (2015). Two R packages for sensitivity analysis in observational studies. Observational Studies, 1(1), 1-17. Free on-line at obsstudies.org.

Rosenbaum, P. R. and Silber, J. H. (2009) Amplification of sensitivity analysis in observational studies. Journal of the American Statistical Association, 104, 1398-1405. <doi:10.1198/jasa.2009.tm08470>

Rosenbaum, P. R. (2017) Observation and Experiment: An Introduction to Causal Inference. Cambridge, MA: Harvard University Press. Amplification is discussed in Chapter 9.

Wolfe, D. A. (1974) A charaterization of population weighted symmetry and related results. Journal of the American Statistical Association, 69, 819-822.

## Examples

 ```1 2 3 4``` ``` data(erpcp) senmv(erpcp,gamma=3.5,trim=1) amplify(3.5,6) amplify(3.5,c(4,6,8,11)) ```

sensitivitymv documentation built on May 2, 2019, 2:06 a.m.