alpha <- 0.025
alpha0 <- 0.01
iuse0 <- 4
iuse1 <- 4
phi0 <- -4
phi1 <- -4
tc0 <- c(3,6,9,12,18)
tc1 <- c(6,12,18,36)
t0 <- (1:5)/5
t1 <- (1:4)/4
rho <- 0.5
delta0 <- 1
delta1 <- 3
sPwRnaiveBonf(alpha=alpha,
alpha0=alpha0,
t1=t1,
delta1=delta1,
iuse1=iuse1,
phi1=phi1)
Run the code above in your browser using DataLab