mscorev: Computes the M-scores used by senmw.

In sensitivitymw: Sensitivity Analysis for Observational Studies Using Weighted M-Statistics

Description

Of limited interest to most users, mscorev computes the M-scores used by senmw for sensitivity analysis in observational studies using Huber-Maritz M-tests. As discussed in Rosenbaum (2014, section 4.1), the ranges of these M-scores within matched sets are ranked.

Usage

mscorev(ymat, inner = 0, trim = 2.5, qu = 0.5, TonT = FALSE)

Arguments

ymat	ymat is a matrix as described in the documentation for senmw or senmv.
inner	inner is the parameter described in the documentation for senmw or senmv. See Rosenbaum (2013).
trim	trim is the parameter described in the documentation for senmw or senmv.
qu	qu is the lambda parameter described in the documentation for senmw or senmv.
TonT	The TonT parameter is not used in the sensitivitymw package. It is included so that the same function mscorev works for both the sensitivitymw and senstivitymv packages. TonT is relevant only when the number of controls varies from one set to another. See the help file for senmv in the sensitivitymv package for discussion of the TonT parameter.

Value

Generally, a matrix with the same dimensions as ymat containing the M-scores.

Note

The example uses the mercury data set. In particular, the first two rows of the mercury data were discussed in the final paragraph of Rosenbaum (2014, Section 2).

Author(s)

Paul R. Rosenbaum

References

Rosenbaum, P. R. (2007) Sensitivity analysis for m-estimates, tests and confidence intervals in matched observational studies. Biometrics, 2007, 63, 456-464. <doi:10.1111/j.1541-0420.2006.00717.x>

Rosenbaum, P. R. (2013) Impact of multiple matched controls on design sensitivity in observational studies. Biometrics, 2013, 69, 118-127. <doi:10.1111/j.1541-0420.2012.01821.x>

Rosenbaum, P. R. (2014) Weighted M-statistics with superior design sensitivity in matched observational studies with multiple controls. Journal of the American Statistical Association, 2014. <doi:10.1080/01621459.2013.879261>

Examples

data(mercury)
mercury[1:2,]
mscorev(mercury,trim=3)[1:2,]

sensitivitymw documentation built on Jan. 4, 2022, 1:09 a.m.