lambda1 <- 0.001982
tl <- 30
tu <- 30*11
N <- 200
HR <- 0.55
tao <- 365*3
A <- 365
t.fixed <- (tl+tu)/2
shape1 <- 5
shape2 <- 5
pow.SEPPLE.random.DE(lambda1, tl, tu, N, HR, tao, A, t.fixed, dist="beta",
shape1=shape1, shape2=shape2, nsim=1000)
Run the code above in your browser using DataLab