combo.wlr: Weighted Logrank Combo Test and Estimation
In keaven/nphsim: Non proportional hazards sample size and simulation

Description Usage Arguments Details Value Examples

Description

1	combo.wlr(survival, cnsr, trt, stra = NULL, fparam)

`survival`	time-to-event variable
`cnsr`	censoring variable: 1=censoring, 0=event
`trt`	treatment varaible. Accepted values are either "experiment" or "control"
`stra`	stratification variable. Currently is set to `NULL`.
`fparam`	a list input for additional function call arguments `fparam$rgs` required: a list of rho gamma pair in the form of `list(c(rho1, gamma1), c(rho2, gamma2), etc...)` `fparam$draw` Number of resampling

Details are provided here

The function return a list with the follow components

pval: One-sided p-Value corresponding to Zmax
rho, gamma: rho and gamma corresponding to Zmax
Zmax: Zmax

medC = 6 
hr <- c(1, 0.6)
intervals <- 3 
gamma <- c(2.5, 5,  7.5,  10) ## a ramp-up enrollment
R     <- c(2  , 2,  2  ,  6 ) ## enrollment period: total of 12 months
eta <- -log(0.99) ## 1% monthly dropout rate
sim1 <- nphsim(nsim=2,lambdaC=log(2)/medC,lambdaE=log(2)/medC*hr, ssC=300,ssE=300,
               intervals=intervals,gamma=gamma, R=R,eta=eta)
               
rgs <- list(c(0, 0), c(0, 1), c(0, 2), c(0, 3), c(1, 0), c(2, 0), c(2, 2))    
draws <- 10        
test1<-simtest(x=sim1,anaD=c(250),method=combo.wlr,fparam=list(rgs=rgs,draws=draws))
test1$result[]

# direct function call (without cutoff)
combo.wlr(survival=sim1$simd$survival, cnsr=sim1$simd$cnsr, trt=sim1$simd$treatment,
          fparam=list(rgs=rgs,draws=draws))