Under-smoothed kernel confidence intervals for the hazard ratio based on censored data

In cancer clinical trials, it is often of interest in estimating the ratios of hazard rates at some specific time points during the study from two independent populations. In this paper, we consider nonparametric confidence interval procedures for the hazard ratio based on kernel estimates for the hazard rates with under-smoothing bandwidths. Two methods are used to derive the confidence intervals: one based on the asymptotic normality of the ratio of the kernel estimates for the hazard rates in two populations and another through Fieller's Theorem. The performances of the proposed confidence intervals are evaluated through Monte-Carlo simulations and applied to the analysis of data from a clinical trial on early breast cancer.     

Symmetric tests and confidence intervals for survival probabilities and quantiles of censored survival data

We describe existing tests and introduce two new tests concerning the value of a survival function. These tests may be used to construct a confidence interval for the survival probability at a given time or for a quantile of the survival distribution. Simulation studies show that error rates can differ substantially from their nominal values, particularly at survival probabilities close to zero or one. We recommend our new constrained bootstrap test for its good overall performance.     

Estimation of Bayesian confidence intervals for the product of conditional survival probabilities

A model is derived to estimate the survival probability of a time interval when censorings occur. The time interval is divided into partial intervals in order to obtain the conditional survival probabilities, each of which is a parameter of a Binomial distributed random variable. To allow for the dependence between the events in the different intervals these parameters are transformed. Corresponding a priori density functions are formulated regarding both the Bayesian uniform distribution and the special model. The a posteriori density function is derived for the product of the conditional survival probabilities, and formulae for the BAYE sian confidence interval and the expectation are given. Lower and upper bounds for the confidence interval and the expectation are derived. Some examples are given to compare the results with other methods.     

A software package for cDNA microarray data normalization and assessing confidence intervals

DNA microarray data are affected by variations from a number of sources. Before these data can be used to infer biological information, the extent of these variations must be assessed. Here we describe an open source software package, lcDNA, that provides tools for filtering, normalizing, and assessing the statistical significance of cDNA microarray data. The program employs a hierarchical Bayesian model and Markov Chain Monte Carlo simulation to estimate gene-specific confidence intervals for each gene in a cDNA microarray data set. This program is designed to perform these primary analytical operations on data from two-channel spotted, or in situ synthesized, DNA microarrays.     

A method for analyzing censored survival phenotype with gene expression data

Background  

Survival time is an important clinical trait for many disease studies. Previous works have shown certain relationship between patients' gene expression profiles and survival time. However, due to the censoring effects of survival time and the high dimensionality of gene expression data, effective and unbiased selection of a gene expression signature to predict survival probabilities requires further study.     

Bayesian inference on age-specific survival for censored and truncated data

1. Traditional estimation of age-specific survival and mortality rates in vertebrates is limited to individuals with known age. Although this subject has been studied extensively using effective capture-recapture and capture-recovery models, inference remains challenging because of large numbers of incomplete records (i.e. unknown age of many individuals) and because of the inadequate duration of the studies. 2. Here, we present a hierarchical model for capture-recapture/recovery (CRR) data sets with large proportions of unknown times of birth and death. The model uses a Bayesian framework to draw inference on population-level age-specific demographic rates using parametric survival functions and applies this information to reconstruct times of birth and death for individuals with unknown age. 3. We simulated a set of CRR data sets with varying study span and proportions of individuals with known age, and varying recapture and recovery probabilities. We used these data sets to compare our method to a traditional CRR model, which requires knowledge of individual ages. Subsequently, we applied our method to a subset of a long-term CRR data set on Soay sheep. 4. Our results show that this method performs better than the common CRR model when sample sizes are low. Still, our model is sensitive to the choice of priors with low recapture probability and short studies. In such cases, priors that overestimate survival perform better than those that underestimate it. Also, the model was able to estimate accurately ages at death for Soay sheep, with an average error of 0.94 years and to identify differences in mortality rate between sexes. 5. Although many of the problems in the estimation of age-specific survival can be reduced through more efficient sampling schemes, most ecological data sets are still sparse and with a large proportion of missing records. Thus, improved sampling needs still to be combined with statistical models capable of overcoming the unavoidable limitations of any fieldwork. We show that our approach provides reliable estimates of parameters and unknown times of birth and death even with the most incomplete data sets while being flexible enough to accommodate multiple recapture probabilities and covariates.     

Testing equality of survival functions based on both paired and unpaired censored data

We introduce two test procedures for comparing two survival distributions on the basis of randomly right-censored data consisting of both paired and unpaired observations. Our procedures are based on generalizations of a pooled rank test statistic previously proposed for uncensored data. One generalization adapts the Prentice-Wilcoxon score, while the other adapts the Akritas score. The use of these particular scoring systems in pooled rank tests with randomly right-censored paired data has been advocated by several researchers. Our test procedures utilize the permutation distributions of the test statistics based on a novel manner of permuting the scores. Permutation versions of tests for right-censored paired data and for two independent right-censored samples that use the proposed scoring systems are obtained as special cases of our test procedures. Simulation results show that our test procedures have high power for detecting scale and location shifts in exponential and log-logistic distributions for the survival times. We also demonstrate the advantages of our test procedures in terms of utilizing randomly occurring unpaired observations that are discarded in test procedures for paired data. The tests are applied to skin graft data previously reported elsewhere.     

A study on confidence intervals for incremental cost-effectiveness ratios

In health policy and economics studies, the incremental cost-effectiveness ratio (ICER) has long been used to compare the economic consequences relative to the health benefits of therapies. Due to the skewed distributions of the costs and ICERs, much research has been done on how to obtain confidence intervals of ICERs, using either parametric or nonparametric methods, with or without the presence of censoring. In this paper, we will examine and compare the finite sample performance of many approaches via simulation studies. For the special situation when the health effect of the treatment is not statistically significant, we will propose a new bootstrapping approach to improve upon the bootstrap percentile method that is currently available. The most efficient way of constructing confidence intervals will be identified and extended to the censored data case. Finally, a data example from a cardiovascular clinical trial is used to demonstrate the application of these methods.     

Nonparametric confidence intervals for the ratio of marginal hazard rates of paired survival times

Paired survival times with potential censoring are often observed from two treatment groups in clinical trials and other types of clinical studies. The ratio of marginal hazard rates may be used to quantify the treatment effect in these studies. In this paper, a recently proposed nonparametric kernel method is used to estimate the marginal hazard rate, and the method of variance estimates recovery (MOVER) is used for the construction of the confidence intervals of a time‐dependent hazard ratio based on the confidence limits of a single marginal hazard rate. Two methods are proposed: one uses the delta method and another adopts the transformation method to construct confidence limits for the marginal hazard rate. Simulations are performed to evaluate the performance of the proposed methods. Real data from two clinical trials are analyzed using the proposed methods.