factorial2x2: Design and Analysis of 2x2 Factorial Trial

Description Usage Arguments Value References See Also Examples

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

1
2
3

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)

`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

`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

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.

eventProb, crit2x2, strLgrkPower power23_13, power13_13_13, power12_12

# 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