pvalci: Common p-value-based procedures: Simultaneous confidence...

Description Usage Arguments Details Value Source References See Also Examples

View source: R/pvalci.R

Description

Computation of simultaneous confidence intervals for selected multiple testing procedures based on univariate p-values (Bonferroni, Holm and fixed-sequence procedures).

Usage

1
pvalci(rawp, est, stderror, weight, covprob, proc)

Arguments

rawp

Vector of raw p-values.

est

Vector of point estimates.

stderror

Vector of standard errors associated with the point estimates.

weight

Vector of hypothesis weights whose sum is equal to 1 (default is a vector of equal weights).

covprob

Simultaneous coverage probability (default is 0.975).

proc

Vector of character strings containing the procedure name. This vector should include any of the following: "Bonferroni", "Holm", "Fixed-sequence".

Details

This function computes one-sided simultaneous confidence limits for the Bonferroni, Holm (Holm, 1979) and fixed-sequence (Westfall and Krishen, 2001) procedures in in general one-sided hypothesis testing problems (equally or unequally weighted null hypotheses).

The simultaneous confidence intervals are computed using the methods developed in Hsu and Berger (1999), Strassburger and Bretz (2008) and Guilbaud (2008). For more information on the algorithms used in the function, see Dmitrienko et al. (2009, Section 2.6).

Value

A data frame result with columns for the raw p-values, point estimates, standard errors, weights, adjusted p-values, and simultaneous confidence limits for each of the procedures.

Source

http://multxpert.com/wiki/MultXpert_package

References

Dmitrienko, A., Bretz, F., Westfall, P.H., Troendle, J., Wiens, B.L., Tamhane, A.C., Hsu, J.C. (2009). Multiple testing methodology. Multiple Testing Problems in Pharmaceutical Statistics. Dmitrienko, A., Tamhane, A.C., Bretz, F. (editors). Chapman and Hall/CRC Press, New York.

Guilbaud, O. (2008). Simultaneous confidence regions corresponding to Holm's stepdown procedure and other closed-testing procedures. Biometrical Journal. 5, 678–692.

Holm, S. (1979). A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics. 6, 65–70.

Hsu, J.C., Berger, R.L. (1999). Stepwise confidence intervals without multiplicity adjustment for dose-response and toxicity studies. Journal of the American Statistical Association. 94, 468–482.

Strassburger, K., Bretz, F. (2008). Compatible simultaneous lower confidence bounds for the Holm procedure and other Bonferroni based closed tests. Statistics in Medicine. 27, 4914–4927.

Westfall, P. H., Krishen, A. (2001). Optimally weighted, fixed sequence, and gatekeeping multiple testing procedures. Journal of Statistical Planning and Inference. 99, 25–40.

See Also

parci

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
# Consider a clinical trial conducted to evaluate the effect of three
# doses of a treatment compared to a placebo with respect to a normally
# distributed endpoint 

# Three null hypotheses of no effect are tested in the trial:
# Null hypothesis H1: No difference between Dose 1 and Placebo
# Null hypothesis H2: No difference between Dose 2 and Placebo
# Null hypothesis H3: No difference between Dose 3 and Placebo

# Null hypotheses of no treatment effect are equally weighted
weight<-c(1/3,1/3,1/3)

# Treatment effect estimates (mean  dose-placebo differences)
est<-c(2.3,2.5,1.9)

# Pooled standard deviation
sd<-9.5

# Study design is balanced with 180 patients per treatment arm
n<-180

# Standard errors
stderror<-rep(sd*sqrt(2/n),3)

# T-statistics associated with the three dose-placebo tests
stat<-est/stderror

# Compute degrees of freedom
nu<-2*(n-1)

# Compute raw one-sided p-values
rawp<-1-pt(stat,nu)

# Compute lower one-sided simultaneous confidence limits 
# for the Bonferroni procedure 
pvalci(rawp,est,stderror,weight,covprob=0.975,proc="Bonferroni")

# Compute lower one-sided simultaneous confidence limits
# for the Holm and Fixed-sequence procedures
pvalci(rawp,est,stderror,weight,covprob=0.975,proc=c("Holm", "Fixed-sequence"))

multxpert documentation built on May 31, 2017, 4:09 a.m.

Related to pvalci in multxpert...