rm.combo.WLRmax: Weighted Logrank combination with asymptotic p-value...

Description Usage Arguments Details Value Author(s) References Examples

View source: R/rm.combo.WLRmax.R

Description

Description

Usage

1
2
3
4
rm.combo.WLRmax(time = NULL, status = NULL, arm = NULL, wt = NULL,
  adjust.methods = c("holm", "hochberg", "hommel", "bonferroni", "asymp")[5],
  ties.method = c("exact", "breslow", "efron")[2], one.sided = FALSE,
  HT.est = FALSE, max = TRUE, alpha = 0.025)

Arguments

time

time-to-event variable

status

event indicator: 0=censored, 1=event

arm

treatment group indicator

wt

a list with pairs of non-negative numbers setting rho and gamma for each Fleming-Harrington weighting sheme to be used

adjust.methods

methods used to adjust Type I error for testing multiple weighted logrank tests

ties.method

methods for dealing with tied event times

one.sided

logical indicating 1-sided (TRUE) vs 2-sided (FALSE) testing

HT.est

logical used when adjust.methods is "asymp", max is TRUE to create an estimate of HR that is averaged across the specified tests

max

indicator of extent of output (TRUE for maximum output)

alpha

level for testing

Details

This routine provides several options for correcting p-values for the maximum of multiple weighted logrank test statistics. The most favorable of these (i.e., minimally sufficient p-value correction) is provided by the "asymp" option which uses a published asymptotic correlation between different weighted logrank tests (Karrison, 2016). However, there are also more conservative (less powerful) options to correct using Holm, Hochberg, Hommel or Bonferroni corrections.

Note that the user needs to load the 'dplyr' and 'Matrix' packages since these have NOT been added to the imports list for the package. The correlation between different weighted logrank test statistics noted in the Karrison (2016) is the same as the variance of a weighted logrank statistic with the sum of rho values and sum of gamma values as the respective rho and gamma parameters.

Value

The function returns a list with the follow components

pval

One-sided p-Value corresponding to Zmax

rho, gamma

rho and gamma corresponding to Zmax

Zmax

Zmax

Author(s)

Satrajit Roychoudhury, [email protected]

References

TG Karrison (2016), Versatile tests for comparing survival curves based on weighted log-rank statistics, Stata Journal, 16:3 pp 678-690

Examples

1
2
3
4
5
with(Ex1delayedEffect,rm.combo.WLRmax(time=month,status=evntd,arm=trt,adjust.methods="asymp",
                                   HT.est=TRUE,wt=list(a1=c(0,1),a2=c(1,0),a3=c(1,1)))) 
with(Ex1delayedEffect,rm.combo.WLRmax(time=month,status=evntd,arm=trt,
                                   wt=list(a1=c(0,1),a2=c(1,0),a3=c(1,1)))) 
                                   

keaven/nphsim documentation built on July 2, 2018, 1:10 a.m.