factorial2x2 (version 0.1.0)

fac2x2design: Power for the 2/3-1/3, 1/3-1/3-1/3, 1/2-1/2 procedures

Description

Power for the 2/3-1/3, 1/3-1/3-1/3, 1/2-1/2 procedures

Usage

fac2x2design(n, rateC, hrA, hrB, hrAB, mincens, maxcens, dig = 2,
  alpha = 0.05, niter = 5, abseps = 0.001, corAa = 1/sqrt(2),
  corAab = 1/sqrt(2), coraab = 1/2)

Arguments

n

total sample size

rateC

group C one year event rate

hrA

group A to group C hazard ratio

hrB

group B to group C hazard ratio

hrAB

group AB to group C hazard ratio

mincens

minimum censoring time

maxcens

maximum censoring time

dig

number of decimal places to roundDown the critical value to

alpha

two-sided significance level

niter

number of times we call pmvnorm to average out its randomness

abseps

abseps setting in the pmvnorm call

corAa

correlation between the overall A and simple A log hazard ratio estimates

corAab

correlation between the overall A and simple AB log hazard ratio estimates

coraab

correalation between the simple A and simple AB log hazard ratio estimates

Value

powerA

power to detect the overall A effect at the two-sided alpha level

powerB

power to detect the overall B effect at the two-sided alpha level

power23.13

power to detect the overall A or simple AB effects using the 2/3-1/3 procedure

power13.13.13

power to detect the overall A, simple A, or simple AB effects using the 1/3-1/3-1/3 procedure

power12.12

power to detect the simple A or simple AB effects using the 1/2-1/2 procedure

events

expected number of events

evtprob

event probabilities for the C, A, B, and AB groups, respectively

References

Leifer, E.S., Troendle, J.F., Kolecki, A., Follmann, D. Joint testing of overall and simple effect for the two-by-two factorial design. (2019). Submitted.

Slud, E.V. Analysis of factorial survival experiments. Biometrics. 1994; 50: 25-38.

See Also

eventProb, crit2x2, strLgrkPower power23_13, power13_13_13, power12_12

Examples

# NOT RUN {
# Corresponds to scenario 5 in Table 2 from Leifer, Troendle, et al. (2019).
n <- 4600
rateC <- 0.0445
hrA <- 0.80
hrB <- 0.80
hrAB <- 0.72
mincens <- 4.0
maxcens <- 8.4

fac2x2design(n, rateC, hrA, hrB, hrAB, mincens, maxcens, dig = 2, alpha = 0.05, niter = 1)
# $powerA
# [1] 0.7182932

# $powerB
# [1] 0.7182932

# $power23.13
# [1] 0.9290062

# $power13.13.13
# [1] 0.9302078

# $power12.12
# [1] 0.9411688

# $events
# [1] 954.8738

# $evtprob
# probC     probA     probB    probAB
# 0.2446365 0.2012540 0.2012540 0.1831806

# }