Induced smoothing for rank regression with censored survival times

Adaptions of weighted rank regression to the accelerated failure time model for censored survival data have been successful in yielding asymptotically normal estimates and flexible weighting schemes to increase statistical efficiencies. However, for only one simple weighting scheme, Gehan or Wilcoxon weights, are estimating equations guaranteed to be monotone in parameter components, and even in this case are step functions, requiring the equivalent of linear programming for computation. The lack of smoothness makes standard error or covariance matrix estimation even more difficult. An induced smoothing technique overcame these difficulties in various problems involving monotone but pure jump estimating equations, including conventional rank regression. The present paper applies induced smoothing to the Gehan-Wilcoxon weighted rank regression for the accelerated failure time model, for the more difficult case of survival time data subject to censoring, where the inapplicability of permutation arguments necessitates a new method of estimating null variance of estimating functions. Smooth monotone parameter estimation and rapid, reliable standard error or covariance matrix estimation is obtained.

A parametric model for the interval censored survival times of Acacia mangium plantation in a spacing trial

Survival times for the Acacia mangium plantation in the Segaliud Lokan Project, Sabah, East Malaysia were analysed based on 20 permanent sample plots (PSPs) established in 1988 as a spacing experiment. The PSPs were established following a complete randomized block design with five levels of spacing randomly assigned to units within four blocks at different sites. The survival times of trees in years are of interest. Since the inventories were only conducted annually, the actual survival time for each tree was not observed. Hence, the data set comprises censored survival times. Initial analysis of the survival of the Acacia mangium plantation suggested there is block by spacing interaction; a Weibull model gives a reasonable fit to the replicate survival times within each PSP; but a standard Weibull regression model is inappropriate because the shape parameter differs between PSPs. In this paper we investigate the form of the non-constant Weibull shape parameter. Parsimonious models for the Weibull survival times have been derived using maximum likelihood methods. The factor selection for the parameters is based on a backward elimination procedure. The models are compared using likelihood ratio statistics. The results suggest that both Weibull parameters depend on spacing and block.

SRL 172 (killed Mycobacterium vaccae) in addition to standard chemotherapy improves quality of life without affecting survival, in patients with advanced non-small-cell lung cancer : phase III results

Background: This open-label, randomised phase III study was designed to further investigate the clinical activity and safety of SRL172 (killed Mycobacterium vaccae suspension) with chemotherapy in the treatment of non-small-cell lung cancer (NSCLC). Patients and methods: Patients were randomised to receive platinum-based chemotherapy, consisting of up to six cycles of MVP (mitomycin, vinblastine and cisplatin or carboplatin) with (210 patients) or without (209 patients) monthly SRL172. Results: There was no statistical difference between the two groups in overall survival (primary efficacy end point) over the course of the study (median overall survival of 223 days versus 225 days; P = 0.65). However, a higher proportion of patients were alive at the end of the 15-week treatment phase in the chemotherapy plus SRL172 group (90%), than in the chemotherapy alone group (83%) (P = 0.061). At the end of the treatment phase, the response rate was 37% in the combined group and 33% in the chemotherapy alone group. Patients in the chemotherapy alone group had greater deterioration in their Global Health Status score (-14.3) than patients in the chemotherapy plus SRL172 group (-6.6) (P = 0.02). Conclusion: In this non-placebo controlled trial, SRL172 when added to standard cancer chemotherapy significantly improved patient quality of life without affecting overall survival times. © 2004 European Society for Medical Oncology.

A boat hitchhiker's guide to survival: Cabomba caroliniana desiccation resistance and survival ability

Cabomba caroliniana is a submersed macrophyte that has become a serious invader. Cabomba predominantly spreads by stem fragments, in particular through unintentional transport on boat trailers ('hitch hiking'). Desiccation resistance affects the potential dispersal radius. Therefore, knowledge of maximum survival times allows predicting future dispersal. Experiments were conducted to assess desiccation resistance and survival ability of cabomba fragments under various environmental scenarios. Cabomba fragments were highly tolerant of desiccation. However, even relatively low wind speeds resulted in rapid mass loss, indicating a low survival rate of fragments exposed to air currents, such as fragments transported on a boat trailer. The experiments indicated that cabomba could survive at least 3 h of overland transport if exposed to wind. However, even small clumps of cabomba could potentially survive up to 42 h. Thus, targeting the transport of clumps of macrophytes should receive high priority in management. The high resilience of cabomba to desiccation demonstrates the risk of continuing spread. Because of the high probability of fragment viability on arrival, preventing fragment uptake on boat trailers is paramount to reduce the risk of further spread. These findings will assist improving models that predict the spread of aquatic invasive macrophytes.

A boat hitchhiker’s guide to survival: Cabomba caroliniana desiccation resistance and survival ability

Cabomba caroliniana is a submersed macrophyte that has become a serious invader. Cabomba predominantly spreads by stem fragments, in particular through unintentional transport on boat trailers (‘hitch hiking’). Desiccation resistance affects the potential dispersal radius. Therefore, knowledge of maximum survival times allows predicting future dispersal. Experiments were conducted to assess desiccation resistance and survival ability of cabomba fragments under various environmental scenarios. Cabomba fragments were highly tolerant of desiccation. However, even relatively low wind speeds resulted in rapid mass loss, indicating a low survival rate of fragments exposed to air currents, such as fragments transported on a boat trailer. The experiments indicated that cabomba could survive at least 3 h of overland transport if exposed to wind. However, even small clumps of cabomba could potentially survive up to 42 h. Thus, targeting the transport of clumps of macrophytes should receive high priority in management. The high resilience of cabomba to desiccation demonstrates the risk of continuing spread. Because of the high probability of fragment viability on arrival, preventing fragment uptake on boat trailers is paramount to reduce the risk of further spread. These findings will assist improving models that predict the spread of aquatic invasive macrophytes.

Tumour markers as prognostic factors of survival of gastric cancer patients

The incidence of gastric cancer in the last decades has declined rapidly in the industrialised countries. Worldwide, however, gastric cancer is still the second most common cause of cancer death. Although surgery is currently the most effective treatment, the rapid progress in adjuvant chemotherapy and radiation therapy requires a re-evaluation of prognosis assessment. The TNM staging system of the UICC is ubiquitously used; it groups patients by decreasing survival times from stage I to stage IV based on the spread of disease, i.e. depth of tumour penetration (T), extent of spread to lymph nodes (N), and the presence or absence of distant (M) metastases. This is by far the most consistent prognostic classification system today. However, even within the stage groups there are patients that follow a varying course of disease. Our knowledge of the molecular differences between tumours of the same stage and morphology has been accumulating over the years and methods for a more accurate assessment of the phenotype of neoplasias are of value when evaluating the prognosis of individual patients with gastric cancer. In this study, the immunohistochemical expression of tumour markers involved in different phases in tumourigenesis was examined. The aim was to find new markers which could provide prognostic information in addition to what is provided by the TNM variables. A total of 337 specimens from the primary tumour of patients who underwent surgery for gastric cancer were collected and the immunohistochemical expression of seven different biomarkers was analysed. DNA ploidy and S-phase fraction (SPF) was assessed by flow cytometry. Finally, all biomarkers and clinicopathological prognostic factors were combined and evaluated by a multivariate Cox regression model to elucidate which specific factors provide independent prognostic information. By univariate survival analysis the following variables were significant prognostic factors: epithelial and stromal syndecan-1 expression, stromal tenascin-C expression, expression of tumour-associated trypsin inhibitor (TATI) in cancer cells, nuclear p53 expression, nuclear p21 expression, DNA ploidy, and SPF. By multivariate survival analysis adjusted for all available clinicopathological and biomolecular variables, p53 expression, p21 expression, and DNA ploidy emerged as independent prognostic biomarkers, together with penetration depth of the tumour, presence of nodal metastases, surgical cure of the cancer, and age of the patient at the time of diagnosis.

A Ten-microRNA Expression Signature Predicts Survival in Glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor with very poor patient median survival. To identify a microRNA (miRNA) expression signature that can predict GBM patient survival, we analyzed the miRNA expression data of GBM patients (n = 222) derived from The Cancer Genome Atlas (TCGA) dataset. We divided the patients randomly into training and testing sets with equal number in each group. We identified 10 significant miRNAs using Cox regression analysis on the training set and formulated a risk score based on the expression signature of these miRNAs that segregated the patients into high and low risk groups with significantly different survival times (hazard ratio HR] = 2.4; 95% CI = 1.4-3.8; p < 0.0001). Of these 10 miRNAs, 7 were found to be risky miRNAs and 3 were found to be protective. This signature was independently validated in the testing set (HR = 1.7; 95% CI = 1.1-2.8; p = 0.002). GBM patients with high risk scores had overall poor survival compared to the patients with low risk scores. Overall survival among the entire patient set was 35.0% at 2 years, 21.5% at 3 years, 18.5% at 4 years and 11.8% at 5 years in the low risk group, versus 11.0%, 5.5%, 0.0 and 0.0% respectively in the high risk group (HR = 2.0; 95% CI = 1.4-2.8; p < 0.0001). Cox multivariate analysis with patient age as a covariate on the entire patient set identified risk score based on the 10 miRNA expression signature to be an independent predictor of patient survival (HR = 1.120; 95% CI = 1.04-1.20; p = 0.003). Thus we have identified a miRNA expression signature that can predict GBM patient survival. These findings may have implications in the understanding of gliomagenesis, development of targeted therapy and selection of high risk cancer patients for adjuvant therapy.

Simulating Coxian Phase-type Distributions for Patient Survival

Coxian phase-type distributions are a special type of Markov model that can be used to represent survival times in terms of phases through which an individual may progress until they eventually leave the system completely. Previous research has considered the Coxian phase-type distribution to be ideal in representing patient survival in hospital. However, problems exist in fitting the distributions. This paper investigates the problems that arise with the fitting process by simulating various Coxian phase-type models for the representation of patient survival and examining the estimated parameter values and eigenvalues obtained. The results indicate that numerical methods previously used for fitting the model parameters do not always converge. An alternative technique is therefore considered. All methods are influenced by the choice of initial parameter values. The investigation uses a data set of 1439 elderly patients and models their survival time, the length of time they spend in a UK hospital.

Intrinsic subtypes of gastric cancer, based on gene expression pattern, predict survival and respond differently to chemotherapy

BACKGROUND & AIMS:
Gastric cancer (GC) is a heterogeneous disease comprising multiple subtypes that have distinct biological properties and effects in patients. We sought to identify new, intrinsic subtypes of GC by gene expression analysis of a large panel of GC cell lines. We tested if these subtypes might be associated with differences in patient survival times and responses to various standard-of-care cytotoxic drugs.
METHODS:
We analyzed gene expression profiles for 37 GC cell lines to identify intrinsic GC subtypes. These subtypes were validated in primary tumors from 521 patients in 4 independent cohorts, where the subtypes were determined by either expression profiling or subtype-specific immunohistochemical markers (LGALS4, CDH17). In vitro sensitivity to 3 chemotherapy drugs (5-fluorouracil, cisplatin, oxaliplatin) was also assessed.
RESULTS:
Unsupervised cell line analysis identified 2 major intrinsic genomic subtypes (G-INT and G-DIF) that had distinct patterns of gene expression. The intrinsic subtypes, but not subtypes based on Lauren's histopathologic classification, were prognostic of survival, based on univariate and multivariate analysis in multiple patient cohorts. The G-INT cell lines were significantly more sensitive to 5-fluorouracil and oxaliplatin, but more resistant to cisplatin, than the G-DIF cell lines. In patients, intrinsic subtypes were associated with survival time following adjuvant, 5-fluorouracil-based therapy.
CONCLUSIONS:
Intrinsic subtypes of GC, based on distinct patterns of expression, are associated with patient survival and response to chemotherapy. Classification of GC based on intrinsic subtypes might be used to determine prognosis and customize therapy.

Low-Dose Aspirin Use After Diagnosis of Colorectal Cancer Does Not Increase Survival: A Case–Control Analysis of a Population-Based Cohort

BACKGROUND & AIMS: Individuals who began taking low-dose aspirin before they were diagnosed with colorectal cancer were reported to have longer survival times than patients who did not take this drug. We investigated survival times of patients who begin taking low-dose aspirin after a diagnosis of colorectal cancer in a large population-based cohort study.

METHODS: We performed a nested case-control analysis using a cohort of 4794 patients diagnosed with colorectal cancer from 1998 through 2007, identified from the UK Clinical Practice Research Datalink and confirmed by cancer registries. There were 1559 colorectal cancer-specific deaths, recorded by the Office of National Statistics; these were each matched with up to 5 risk-set controls. Conditional logistic regression was used to calculate odds ratios (OR) and 95% confidence intervals (CI), based on practitioner-recorded aspirin usage.

RESULTS: Overall, low-dose aspirin use after a diagnosis of colorectal cancer was not associated with colorectal cancer-specific mortality (adjusted OR = 1.06; 95% CI: 0.92-1.24) or all-cause mortality (adjusted OR = 1.06; 95% CI: 0.94-1.19). A dose-response association was not apparent; for example, low-dose aspirin use for more than 1 year after diagnosis was not associated with colorectal cancer-specific mortality (adjusted OR = 0.98; 95% CI: 0.82-1.19). There was also no association between low-dose aspirin usage and colon cancer-specific mortality (adjusted OR = 1.02; 95% CI: 0.83-1.25) or rectal cancer-specific mortality (adjusted OR = 1.10; 95% CI: 0.88-1.38).

CONCLUSIONS: In a large population-based cohort, low-dose aspirin usage after diagnosis of colorectal cancer did not increase survival time.

Parametric accelerated failure time models with random effects and an application to kidney transplant survival

Accelerated failure time models with a shared random component are described, and are used to evaluate the effect of explanatory factors and different transplant centres on survival times following kidney transplantation. Different combinations of the distribution of the random effects and baseline hazard function are considered and the fit of such models to the transplant data is critically assessed. A mixture model that combines short- and long-term components of a hazard function is then developed, which provides a more flexible model for the hazard function. The model can incorporate different explanatory variables and random effects in each component. The model is straightforward to fit using standard statistical software, and is shown to be a good fit to the transplant data. Copyright (C) 2004 John Wiley Sons, Ltd.

A Bayesian Analysis in the Presence of Covariates for Multivariate Survival Data: An example of Application

In this paper, we introduce a Bayesian analysis for survival multivariate data in the presence of a covariate vector and censored observations. Different ""frailties"" or latent variables are considered to capture the correlation among the survival times for the same individual. We assume Weibull or generalized Gamma distributions considering right censored lifetime data. We develop the Bayesian analysis using Markov Chain Monte Carlo (MCMC) methods.

A bivariate regression model for matched paired survival data: local influence and residual analysis

The use of bivariate distributions plays a fundamental role in survival and reliability studies. In this paper, we consider a location scale model for bivariate survival times based on the proposal of a copula to model the dependence of bivariate survival data. For the proposed model, we consider inferential procedures based on maximum likelihood. Gains in efficiency from bivariate models are also examined in the censored data setting. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and compared to the performance of the bivariate regression model for matched paired survival data. Sensitivity analysis methods such as local and total influence are presented and derived under three perturbation schemes. The martingale marginal and the deviance marginal residual measures are used to check the adequacy of the model. Furthermore, we propose a new measure which we call modified deviance component residual. The methodology in the paper is illustrated on a lifetime data set for kidney patients.