testtwice-package {testtwice}

Testing One Hypothesis Twice in Observational Studies

Description

Tests one hypothesis with several test statistics, correcting for multiple testing. The central function in the package is testtwice(). In a sensitivity analysis, the method has the largest design sensitivity of its component tests. The package implements the method and examples in Rosenbaum, P. R. (2012) <doi:10.1093/biomet/ass032> Testing one hypothesis twice in observational studies. Biometrika, 99(4), 763-774.

Details

The DESCRIPTION file:

Package:	testtwice
Type:	Package
Title:	Testing One Hypothesis Twice in Observational Studies
Version:	1.0.3
Author:	Paul R. Rosenbaum
Maintainer:	Paul R. Rosenbaum <[e-mail]>
Description:	Tests one hypothesis with several test statistics, correcting for multiple testing. The central function in the package is testtwice(). In a sensitivity analysis, the method has the largest design sensitivity of its component tests. The package implements the method and examples in Rosenbaum, P. R. (2012) <doi:10.1093/biomet/ass032> Testing one hypothesis twice in observational studies. Biometrika, 99(4), 763-774.
License:	GPL-2
Encoding:	UTF-8
LazyData:	true
Imports:	stats, mvtnorm

Index of help topics:

bmhranks                One-step and Two-step Signed Ranks.
multrnk                 A family of Signed Ranks for Matched Pair
                        Differences.
smokerlead              Pair Differences in Blood Lead Levels in 679
                        Matched Pairs of One Smoker and One Control.
testtwice               Computes the P-value and Sensitivity Bound for
                        Testing Twice.
testtwice-package       Testing One Hypothesis Twice in Observational
                        Studies
tt                      Computes the P-value and Sensitivity Bound for
                        Testing Twice From a User Supplied Matrix.

Author(s)

Paul R. Rosenbaum

Maintainer: Paul R. Rosenbaum <[e-mail]>

References

Rosenbaum, P. R. (2012) <doi:10.1093/biomet/ass032> Testing one hypothesis twice in observational studies. Biometrika, 99(4), 763-774.

Examples

data(smokerlead)
attach(smokerlead)

# The following example reproduces parts of the first column of Table 3 in Rosenbaum (2012).
w<-rank(abs(lead))
wd<-rank(abs(lead))*dose
wl<-rank(abs(llead))
H<-cbind(w,wd,wl)
tt(lead,H,gamma=2.8)
tt(lead,H,gamma=3.2)

rm(w,wd,wl,H)
detach(smokerlead)