powDT: To calculate the testing power for step-down or step-up...

In DunnettTests: Software implementation of step-down and step-up Dunnett test procedures

Description

For the problem of comparing means of k treatment groups to the mean of one control group. The implementation of the function needs the following three assumptions: 1. The k treatment groups have identical treatment effect size. 2. The samples to be assigned to each of the k treatment groups are expected to be equal at size n. 3. The alternative hypotheses are one-sided. With the violations of either of the assumptions, simulation-based power evaluation is suggested.

Usage

powDT(r,k,mu,mu0,n,n0,contrast,sigma=NA,df=Inf,alpha=0.05,mcs=1e+05,testcall)

Arguments

r	The least number of null hypotheses to be rejected, e.g.,when r=1, the disjunctive power is returned and when r=k, the conjunctive power is returned.
k	Number of hypotheses to be tested, k≥ 2 and k≤ 16.
mu	Assumed population mean in each of the k treatment groups.
mu0	Assumed population mean in the control group.
n	Sample size in each of the k treatment groups
n0	Sample size in the control group
contrast	If mu and mu0 are concerned of mean of a continous outcome, specify contrast="means"; if mu and mu0 are concerned of proportion of binary outcome, specify contrast="props".
sigma	The population error variance, which should be specified when contrast="means"; if contrast="props", set sigma=NA as default and it will be calculated based on mu and mu0 specified within the function.
df	Degree of freedom of the t-test statistics. When (approximately) normally distributed test statistics are applied, set df=Inf (default).
alpha	The pre-specified overall significance level, default=0.05.
mcs	The number of monte-carlo sample points to numerically approximate the power for a given sample size, refer to Equation (4.3) and Equation (4.5) in Dunnett and Tamhane (1992).
testcall	The applied Dunnett test procedure: "SD"=step-down Dunnett test; "SU"=step-up Dunnett test.

Value

Return the power of rejecting at least r out of the k false null hypotheses.

Author(s)

FAN XIA <[e-mail]>

References

Charles W. Dunnett and Ajit C. Tamhane. A step-up multiple test procedure. Journal of the American Statistical Association, 87(417):162-170, 1992.

See Also

nvDT

Examples

#Compare group means of four treatment arms to a control arm (upper one-sided tests)
# Setting
k <- 4
mu <- 22 #assumed mean of each treatment arm
mu0 <- 20 #assumed mean of the control arm
n <- 100
n0 <- 80
sigma <- 5 #assumed population error variance
df <- n*k+n0-k-1 #consider the t distribution

# at one-sided significance level 0.05
# compare the testing powers of SD and SU Dunnett for rejecting at least one nulls 
powSD <- powDT(r=1,k,mu,mu0,n,n0,"means",sigma=sigma,df=df,testcall="SD")
powSU <- powDT(r=1,k,mu,mu0,n,n0,"means",sigma=sigma,df=df,testcall="SU")

DunnettTests documentation built on May 2, 2019, 9:13 a.m.