Description Usage Arguments Details Value Author(s) References Examples

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

Description

1 2 3 4 |

`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 |

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.

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

Satrajit Roychoudhury, [email protected]

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

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.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.