Title: | Two Stage Hazard Rate Comparison |
---|---|
Description: | Two-stage procedure compares hazard rate functions, which may or may not cross each other. |
Authors: | Jun Sheng, Peihua Qiu <[email protected]>, and Charles J. Geyer <[email protected]> |
Maintainer: | Charles J. Geyer <[email protected]> |
License: | MIT + file LICENSE |
Version: | 0.1-6 |
Built: | 2025-03-10 02:16:20 UTC |
Source: | https://github.com/cran/TSHRC |
Rats were taken from 50 distinct litters and one rat from the litter was randomly selected and given the drug. For each litter two rats were selected as controls and were given a placebo. In the treatment group, there are 29 censored observations of the times to tumor.
data(rats)
data(rats)
A data frame with records for 150 rats.
Time to tumor or censoring.
Zero or one, zero indicates censoring.
Zero or one, one indicates treatment group.
Qiu, P. and Sheng, J. (2008). A two-stage procedure for comparing hazard rate functions. Journal of the Royal Statistical Society: Series B, 70:191-208.
Mantel, N., Bohidar, N. R. and Ciminera, J. L. (1977). Mantel-Haenszel analysis of litter-matched time-to-response data, with modifications for recovery of interlitter information. Cancer Research, 37:3863-3868.
library(TSHRC) data(rats)
library(TSHRC) data(rats)
Two-stage procedure for comparing hazard rate functions, especially suited for situation where hazard rate functions cross.
twostage(time, delta, group, nboot, alpha = 0.05, eps = 0.1)
twostage(time, delta, group, nboot, alpha = 0.05, eps = 0.1)
time |
a vector of type |
delta |
a vector of type |
group |
a vector of type |
nboot |
number of bootstrap samples. Must be positive integer. |
alpha |
significance level, must be between zero and one. |
eps |
small number used in defining the test statistic. Must be between zero and one. See cited paper for discussion. |
This procedure compares two hazard rate functions in two stages. In the first stage, the conventional log-rank test is performed, which is powerful only when the two hazard rate functions do not cross each other. When the two hazard rate functions cross each other, positive differences and negative differences between the two functions would be canceled out in the log-rank test statistic, resulting in ineffective comparison. If the log-rank test gives a significant result, then the entire two-stage procedure stops and we conclude that the two hazard rate functions are significantly different. Otherwise, the stage-II test is performed, which is designed specifically for detecting crossing difference between the two hazard rate functions and has the property that its test statistic is independent of the log-rank test statistic. The independence property is used in properly defining the p-value of the two-stage procedure, based on the p-values of the two individual tests.
vector of P-values. "LRPV"
: p-value of the log-rank test,
"MTPV"
: p-value of the suggested stage-II test,
"TSPV"
: p-value of the two-stage test.
Qiu, P. and Sheng, J. (2008). A two-stage procedure for comparing hazard rate functions. Journal of the Royal Statistical Society: Series B, 70:191-208.
library(TSHRC) data(rats) attach(rats) twostage(time, delta, group, nboot = 100)
library(TSHRC) data(rats) attach(rats) twostage(time, delta, group, nboot = 100)