Regression with bivariate grouped data

The problem of accounting for the grouping of continuous, bivariate data in regression analyses is considered. Reasons why grouping must be taken seriously are advanced, and a strategy for accounting for grouping is demonstrated. The specific model asserts that, in the absence of grouping, the data would be bivariate normal. This model is used to adjust estimates of parameters in a regression relating disease severity to a grouped exposure variable, using data on pneumoconiosis in English coal miners (Ashford, 1959, Biometrics 15, 573-581). The choice of computing methods is discussed and likelihood formulas are presented.     

Distribution-free regression analysis of grouped survival data

Methods based on regression models for logarithmic hazard functions, Cox models, are given for analysis of grouped and censored survival data. By making an approximation it is possible to obtain explicitly a maximum likelihood function involving only the regression parameters. This likelihood function is a convenient analog to Cox's partial likelihood for ungrouped data. The method is applied to data from a toxicological experiment.     

Regression with frailty in survival analysis

In studies of survival, the hazard function for each individual may depend on observed risk variables but usually not all such variables are known or measurable. This unknown factor of the hazard function is usually termed the individual heterogeneity or frailty. When survival is time to the occurrence of a particular type of event and more than one such time may be obtained for each individual, frailty is a common factor among such recurrence times. A model including frailty is fitted to such repeated measures of recurrence times.     

Analysis of doubly-censored survival data, with application to AIDS

This paper proposes nonparametric and weakly structured parametric methods for analyzing survival data in which both the time origin and the failure event can be right- or interval-censored. Such data arise in clinical investigations of the human immunodeficiency virus (HIV) when the infection and clinical status of patients are observed only at several time points. The proposed methods generalize the self-consistency algorithm proposed by Turnbull (1976, Journal of the Royal Statistical Society, Series B 38, 290-295) for singly-censored univariate data, and are illustrated with the results from a study of hemophiliacs who were infected with HIV by contaminated blood factor.     

Control of data quality for population-based cancer survival analysis

BackgroundPopulation-based cancer survival is an important measure of the overall effectiveness of cancer care in a population. Population-based cancer registries collect data that enable the estimation of cancer survival. To ensure accurate, consistent and comparable survival estimates, strict control of data quality is required before the survival analyses are carried out. In this paper, we present a basis for data quality control for cancer survival.MethodsWe propose three distinct phases for the quality control. Firstly, each individual variable within a given record is examined to identify departures from the study protocol; secondly, each record is checked and excluded if it is ineligible or logically incoherent for analysis; lastly, the distributions of key characteristics in the whole dataset are examined for their plausibility.ResultsData for patients diagnosed with bladder cancer in England between 1991 and 2010 are used as an example to aid the interpretation of the differences in data quality. The effect of different aspects of data quality on survival estimates is discussed.ConclusionsWe recommend that the results of data quality procedures should be reported together with the findings from survival analysis, to facilitate their interpretation.     

Dimension reduction methods for microarrays with application to censored survival data

MOTIVATION: Recent research has shown that gene expression profiles can potentially be used for predicting various clinical phenotypes, such as tumor class, drug response and survival time. While there has been extensive studies on tumor classification, there has been less emphasis on other phenotypic features, in particular, patient survival time or time to cancer recurrence, which are subject to right censoring. We consider in this paper an analysis of censored survival time based on microarray gene expression profiles. RESULTS: We propose a dimension reduction strategy, which combines principal components analysis and sliced inverse regression, to identify linear combinations of genes, that both account for the variability in the gene expression levels and preserve the phenotypic information. The extracted gene combinations are then employed as covariates in a predictive survival model formulation. We apply the proposed method to a large diffuse large-B-cell lymphoma dataset, which consists of 240 patients and 7399 genes, and build a Cox proportional hazards model based on the derived gene expression components. The proposed method is shown to provide a good predictive performance for patient survival, as demonstrated by both the significant survival difference between the predicted risk groups and the receiver operator characteristics analysis. AVAILABILITY: R programs are available upon request from the authors. SUPPLEMENTARY INFORMATION: http://dna.ucdavis.edu/~hli/bioinfo-surv-supp.pdf.     

Maximum likelihood estimation for a special type of grouped data with an application to a dose-response problem.

A short-cut method is given for calculating grouped maximum likelihood (ML) estimates when the data are relatively coarsely grouped in some directions, but more finely grouped in others. The algebraic details are then worked out for a dose-response problem that generates data of this kind. The situation envisaged is a variation on the usual quantal response problem in that dosage levels are taken to be random but grouped. Finally, the method is applied both to real and simulated response data conforming to this pattern and shown to work well in practice.     

Regression analysis of mean quality-adjusted lifetime with censored data

In clinical trials of chronic diseases such as acquired immunodeficiency syndrome, cancer, or cardiovascular diseases, the concept of quality-adjusted lifetime (QAL) has received more and more attention. In this paper, we consider the problem of how the covariates affect the mean QAL when the data are subject to right censoring. We allow a very general form for the mean model as a function of covariates. Using the idea of inverse probability weighting, we first construct a simple weighted estimating equation for the parameters in our mean model. We then find the form of the most efficient estimating equation, which yields the most efficient estimator for the regression parameters. Since the most efficient estimator depends on the distribution of the health history processes, and thus cannot be estimated nonparametrically, we consider different approaches for improving the efficiency of the simple weighted estimating equation using observed data. The applicability of these methods is demonstrated by both simulation experiments and a data example from a breast cancer clinical trial study.     

Neural network analysis of follow-up data in primary breast cancer

This paper reports on the performance of a recently developed neural network environment incorporating likelihood-based optimization and complexity reduction techniques in the analysis of breast cancer follow-up data with the goal of building up a clinical decision support system. The inputs to the neural network include classical factors such as grading, age, tumor size, estrogen and progesterone receptor measurements, as well as tumor biological markers such as PAI-1 and uPA. The network learns the structural relationship between these factors and the follow-up data. Examples of neural models for relapse-free survival are presented, which are based on data from 784 breast cancer patients who received their primary therapy at the Department of Obstetrics and Gynecology, Technische Universit?t München, Germany. The performance of the neural analysis as quantified by various indicators (likelihood, Kaplan-Meier curves, log-rank tests) was very high. For example, dividing the patients into two equally sized groups based on the neural score (i.e., cutoff = median score) leads to an estimated difference in relapse-free survival of 40% or better (80% vs. 40%) after 10 years in Kaplan-Meier analysis. Evidence for factor interactions as well as for time-varying impacts is presented. The neural network weights included in the models are significant at the 5% level. The use of neural network analysis and scoring in combination with strong tumor biological factors such as uPA and PAI-1 appears to result in a very effective risk group discrimination. Considerable additional comparison of data from different patient series will be required to establish the generalization capability more firmly. Nonetheless, the improvement of risk group discrimination represents an important step toward the use of neural networks for decision support in a clinical framework and in making the most of biological markers.     

Differential coexpression analysis using microarray data and its application to human cancer

MOTIVATION: Microarrays have been used to identify differential expression of individual genes or cluster genes that are coexpressed over various conditions. However, alteration in coexpression relationships has not been studied. Here we introduce a model for finding differential coexpression from microarrays and test its biological validity with respect to cancer. RESULTS: We collected 10 published gene expression datasets from cancers of 13 different tissues and constructed 2 distinct coexpression networks: a tumor network and normal network. Comparison of the two networks showed that cancer affected many coexpression relationships. Functional changes such as alteration in energy metabolism, promotion of cell growth and enhanced immune activity were accompanied with coexpression changes. Coregulation of collagen genes that may control invasion and metastatic spread of tumor cells was also found. Cluster analysis in the tumor network identified groups of highly interconnected genes related to ribosomal protein synthesis, the cell cycle and antigen presentation. Metallothionein expression was also found to be clustered, which may play a role in apoptosis control in tumor cells. Our results show that this model would serve as a novel method for analyzing microarrays beyond the specific implications for cancer.     

Proteomic analysis of human lysosomes: application to monocytic and breast cancer cells

To date, about fifty lysosomal hydrolases have been identified, and most of them are targeted towards the lysosomes through a specific mannose-6-phosphate (M-6-P) tag. As more lysosomal hydrolases were expected to be discovered, we performed a proteomic study of soluble lysosomal proteins. Human cells were induced to secrete M-6-P proteins which were affinity purified on immobilized M-6-P receptor. The purified proteins were resolved by two-dimensional electrophoresis and analyzed by mass spectrometry. Twenty-two proteins were identified, among which 16 were well-known lysosomal hydrolases. The remaining species distributed as follows: epididymis-specific alpha-mannosidase is a new mannosidase homolog, cystatin F and CREG (cellular repressor of E1A-stimulated genes) were previously identified as M-6-P proteins (Journet et al., Electrophoresis 2000, 21, 3411-3419), and the last three, which are not hydrolases, were up to now considered as nonlysosomal. This two-dimensional reference map of human U937 M-6-P proteins was afterwards used for comparison with M-6-P proteins purified either from U937 differentiated into macrophage-like cells, or from human breast cancer MCF7 cells. Phorbol ester induced differentiation of U937 cells led to limited proteolytic cleavage or maturation of a discrete number of hydrolases. Five additional lysosomal hydrolases were identified from MCF7 samples. These results prove the usefulness of such a procedure to analyze the lysosomal content of various cell lines, to discover new M-6-P proteins, as well as to point towards unknown biological processes.     

Regression analysis of multivariate panel count data

We consider panel count data which are frequently obtained in prospective studies involving recurrent events that are only detected and recorded at periodic assessment times. The data take the form of counts of the cumulative number of events detected at each inspection time, along with explanatory covariates. Examples arise in diverse areas such as epidemiological studies, medical follow-up studies, reliability studies, and tumorigenicity experiments. This article is concerned with regression analysis of multivariate panel count data which arise if more than one type of recurrent event is of interest and individuals are only observed intermittently. We present a class of marginal mean models which leave the dependence structures for related types of recurrent events completely unspecified. Estimating equations are developed for regression parameters, and the resulting estimates are shown to be consistent and asymptotically normal. Simulation studies show that the proposed estimation procedures work well for practical situations. The methodology is applied to a motivating study of patients with psoriatic arthritis in which the events of interest are the onset of joint damage according to 2 different criteria.