Bivariate frailty model for the analysis of multivariate survival time

Because of limitations of the univariate frailty model in analysis of multivariate survival data, a bivariate frailty model is introduced for the analysis of bivariate survival data. This provides tremendous flexibility especially in allowing negative associations between subjects within the same cluster. The approach involves incorporating into the model two possibly correlated frailties for each cluster. The bivariate lognormal distribution is used as the frailty distribution. The model is then generalized to multivariate survival data with two distinguished groups and also to alternating process data. A modified EM algorithm is developed with no requirement of specification of the baseline hazards. The estimators are generalized maximum likelihood estimators with subject-specific interpretation. The model is applied to a mental health study on evaluation of health policy effects for inpatient psychiatric care.     

Bayesian analysis of multivariate survival data using Monte Carlo methods

This paper deals with the analysis of multivariate survival data from a Bayesian perspective using Markov-chain Monte Carlo methods. The Metropolis along with the Gibbs algorithm is used to calculate some of the marginal posterior distributions. A multivariate survival model is proposed, since survival times within the same group are correlated as a consequence of a frailty random block effect. The conditional proportional-hazards model of Clayton and Cuzick is used with a martingale structured prior process (Arjas and Gasbarra) for the discretized baseline hazard. Besides the calculation of the marginal posterior distributions of the parameters of interest, this paper presents some Bayesian EDA diagnostic techniques to detect model adequacy. The methodology is exemplified with kidney infection data where the times to infections within the same patients are expected to be correlated.     

A Bayesian nonparametric estimator of a multivariate survival function

Using reinforced processes related to beta-Stacy process and generalized Pólya urn scheme jointly with a structure assumption about dependence, a Bayesian nonparametric prior and a predictive estimator for a multivariate survival function are provided. This estimator can be computed through an easy implementation of a Gibbs sampler algorithm. Moreover consistency of the estimator is studied.     

Urn Sampling and the Proportional Hazard Model

This paper investigates the urn sampling analogue for the score statistic relating survival to covariates assuming a proportional hazard model. The exact permutation distribution can be calculated as well as the exact low order moments for arbitrary censoring patterns. The asymptotic distribution of the score statistic is an easy consequence. The method is naturally extended to deal with the multivariate case, time varying covariates and interval censoring. Finally the relationship between the censoring process, the survival times and covariates are studied considering different reference sets for the distribution of the score statistic. Some assumptions about the censoring process are investigated and as a consequence the effect of censoring is clarified.     

The multivariate Markov and multiple Chebyshev inequalities

Abstract

A sharp probability inequality named the multivariate Markov inequality is derived for the intersection of the survival functions for non-negative random variables as an extension of the Markov inequality for a single variable. The corresponding result in Chebyshev’s inequality is also obtained as a special case of the multivariate Markov inequality, which is called the multiple Chebyshev inequality to distinguish from the multivariate Chebyshev inequality for a quadratic form of standardized uncorrelated variables. Further, the results are extended to the inequalities for the union of the survival functions and those with lower bounds.     

A Composite Likelihood Approach to Multivariate Survival Data

This paper is about the statistical analysis of multivariate survival data. We discuss the additive and multiplicative frailty models which have been the most popular models for multivariate survival data. As an alternative to the additive and multiplicative frailty models, we propose basing inference on a composite likelihood function that only requires modelling of the marginal distribution of pairs of failure times. Each marginal distribution of a pair of failure times is here assumed to follow a shared frailty model. The method is illustrated with a real-life example.     

A Bayesian positive dependence of survival times based on the multivariate arrangement increasing property

We introduce a new notion of positive dependence of survival times of system components using the multivariate arrangement increasing property. Following the spirit of Barlow and Mendel (J. Amer. Statist. Assoc. 87, 1116–1122), who introduced a new univariate aging notion relative to exchangeable populations of components, we characterize a multivariate positive dependence with respect to exchangeable multicomponent systems. Closure properties of such a class of distributions under some reliability operations are discussed. For an infinite population of systems our definition of multivariate positive dependence can be considered in the frequentist’s paradigm as multivariate totally positive of order 2 with an independence condition. de Finetti(-type) representations for a particular class of survival functions are also given.     

Study of the bivariate survival data using frailty models based on Lévy processes

Frailty models allow us to take into account the non-observable inhomogeneity of individual hazard functions. Although models with time-independent frailty have been intensively studied over the last decades and a wide range of applications in survival analysis have been found, the studies based on the models with time-dependent frailty are relatively rare. In this paper, we formulate and prove two propositions related to the identifiability of the bivariate survival models with frailty given by a nonnegative bivariate Lévy process. We discuss parametric and semiparametric procedures for estimating unknown parameters and baseline hazard functions. Numerical experiments with simulated and real data illustrate these procedures. The statements of the propositions can be easily extended to the multivariate case.

     

Bivariate modelling of longitudinal measurements of two human immunodeficiency type 1 disease progression markers in the presence of informative drop-outs

Summary.  The main statistical problem in many epidemiological studies which involve repeated measurements of surrogate markers is the frequent occurrence of missing data. Standard likelihood-based approaches like the linear random-effects model fail to give unbiased estimates when data are non-ignorably missing. In human immunodeficiency virus (HIV) type 1 infection, two markers which have been widely used to track progression of the disease are CD4 cell counts and HIV–ribonucleic acid (RNA) viral load levels. Repeated measurements of these markers tend to be informatively censored, which is a special case of non-ignorable missingness. In such cases, we need to apply methods that jointly model the observed data and the missingness process. Despite their high correlation, longitudinal data of these markers have been analysed independently by using mainly random-effects models. Touloumi and co-workers have proposed a model termed the joint multivariate random-effects model which combines a linear random-effects model for the underlying pattern of the marker with a log-normal survival model for the drop-out process. We extend the joint multivariate random-effects model to model simultaneously the CD4 cell and viral load data while adjusting for informative drop-outs due to disease progression or death. Estimates of all the model's parameters are obtained by using the restricted iterative generalized least squares method or a modified version of it using the EM algorithm as a nested algorithm in the case of censored survival data taking also into account non-linearity in the HIV–RNA trend. The method proposed is evaluated and compared with simpler approaches in a simulation study. Finally the method is applied to a subset of the data from the 'Concerted action on seroconversion to AIDS and death in Europe' study.     

SPC for short-run multivariate autocorrelated processes

This paper discusses the development of a multivariate control charting technique for short-run autocorrelated data manufacturing environment. The proposed approach is a combination of the multivariate residual charts for autocorrelated data and the multivariate transformation technique for i.i.d. process observations of short lengths. The proposed approach consists in fitting adequate multivariate time-series model of various process outputs and computes the residuals, transforming them into standard normal N(0, 1) data and then using standardized data as inputs to plot conventional univariate i.i.d. control charts. The objective for applying multivariate finite horizon techniques for autocorrelated processes is to allow continuous process monitoring, since all process outputs are controlled trough the use of a single control chart with constant control limits. Throughout simulated examples, it is shown that the proposed short-run process monitoring technique provides approximately similar shifts detection properties as VAR residual charts.     

Control chart for monitoring multivariate COM-Poisson attributes

Statistical process control of multi-attribute count data has received much attention with modern data-acquisition equipment and online computers. The multivariate Poisson distribution is often used to monitor multivariate attributes count data. However, little work has been done so far on under- or over-dispersed multivariate count data, which is common in many industrial processes, with positive or negative correlation. In this study, a Shewhart-type multivariate control chart is constructed to monitor such kind of data, namely the multivariate COM-Poisson (MCP) chart, based on the MCP distribution. The performance of the MCP chart is evaluated by the average run length in simulation. The proposed chart generalizes some existing multivariate attribute charts as its special cases. A real-life bivariate process and a simulated trivariate Poisson process are used to illustrate the application of the MCP chart.     

A Multivariate Model for Repeated Failure Time Measurements

A parametric multivariate failure time distribution is derived from a frailty-type model with a particular frailty distribution. It covers as special cases certain distributions which have been used for multivariate survival data in recent years. Some properties of the distribution are derived: its marginal and conditional distributions lie within the parametric family, and association between the component variates can be positive or, to a limited extent, negative. The simple closed form of the survivor function is useful for right-censored data, as occur commonly in survival analysis, and for calculating uniform residuals. Also featured is the distribution of ratios of paired failure times. The model is applied to data from the literature     

Dependence Estimation for Marginal Models of Multivariate Survival Data

We have previously(Segal and Neuhaus, 1993) devised methods for obtaining marginal regression coefficients and associated variance estimates for multivariate survival data, using a synthesis of the Poisson regression formulation for univariate censored survival analysis and generalized estimating equations (GEE's). The method is parametric in that a baseline survival distribution is specified. Analogous semiparametric models, with unspecified baseline survival, have also been developed (Wei, Lin and Weissfeld, 1989; Lin, 1994).Common to both these approaches is the provision of robust variances for the regression parameters. However, none of this work has addressed the more difficult area of dependence estimation. While GEE approaches ostensibly provide such estimates, we show that there are problems adopting these with multivariate survival data. Further, we demonstrate that these problems can affect estimation of the regression coefficients themselves. An alternate, ad hoc approach to dependence estimation, based on design effects, is proposed and evaluated via simulation and illustrative examples. This revised version was published online in July 2006 with corrections to the Cover Date.     

Multivariate process capability using principal component analysis in the presence of measurement errors

Various different definitions of multivariate process capability indices have been proposed in the literature. Most of the research works related to multivariate process capability indices assume no gauge measurement errors. However, in industrial applications, despite the use of highly advanced measuring instruments, account needs to be taken of gauge imprecision. In this paper we are going to examine the effects of measurement errors on multivariate process capability indices computed using the principal components analysis. We show that measurement errors alter the results of a multivariate process capability analysis, resulting in either a decrease or an increase in the capability of the process. In order to achieve accurate process capability assessments, we propose a method useful for overcoming the effects of gauge measurement errors.     

On improvement of prediction in arma processes

Necessary and sufficient conditions are derived in the paper that enable to decide whether an additional multivariate process will improve the prediction in a given multivariate discrete stationary process. The both processes are assumed to form together a process ARMAm n Further it was investigated wnen one can asser t that the both processes are uncorrelated provided the additional process did not improve the prediction in the original process, Some hints for the actual construction of predictors in a multivariate ARMA. (m n) process can be found in the paper.     

Continuous Time Wishart Process for Stochastic Risk

Risks are usually represented and measured by volatility-covolatility matrices. Wishart processes are models for a dynamic analysis of multivariate risk and describe the evolution of stochastic volatility-covolatility matrices, constrained to be symmetric positive definite. The autoregressive Wishart process (WAR) is the multivariate extension of the Cox, Ingersoll, Ross (CIR) process introduced for scalar stochastic volatility. As a CIR process it allows for closed-form solutions for a number of financial problems, such as term structure of T-bonds and corporate bonds, derivative pricing in a multivariate stochastic volatility model, and the structural model for credit risk. Moreover, the Wishart dynamics are very flexible and are serious competitors for less structural multivariate ARCH models.     

Diagnostic Plots for Assessing the Frailty Distribution in Multivariate Survival Data

In biomedical studies, frailty models arecommonly used in analyzing multivariate survival data, wherethe objective of the study is to estimate both the covariateeffect and the dependence between the multivariate survival times.However, inference based on these models are dependent on thedistributional assumption of frailty. We propose a diagnosticplot for assessing the frailty assumption. The proposed methodis based on the cross-ratio function and the diagnostic plotsuggested by Oakes (1989). We use kernel regression smoothingwith bandwidth choice by cross-validation, to obtain the proposedplot. The resulting plot is capable of differentiating betweenthe gamma and positive stable frailty models when strong associationis present. We illustrate the feasibility of our method usingsimulation studies under known frailty distributions. The approachis applied to data on blindness for each eye of diabetic patientswith adult onset diabetes and a reasonable fit to the gamma frailtymodel is found.     

Continuous Time Wishart Process for Stochastic Risk

Risks are usually represented and measured by volatility–covolatility matrices. Wishart processes are models for a dynamic analysis of multivariate risk and describe the evolution of stochastic volatility–covolatility matrices, constrained to be symmetric positive definite. The autoregressive Wishart process (WAR) is the multivariate extension of the Cox, Ingersoll, Ross (CIR) process introduced for scalar stochastic volatility. As a CIR process it allows for closed-form solutions for a number of financial problems, such as term structure of T-bonds and corporate bonds, derivative pricing in a multivariate stochastic volatility model, and the structural model for credit risk. Moreover, the Wishart dynamics are very flexible and are serious competitors for less structural multivariate ARCH models.     

Stochastic comparisons of multivariate frailty models

A multivariate frailty model in which survival function depends on baseline distributions of components and the frailty random variable is considered. Since misspecification in choice of frailty distribution and/or baseline distribution may affect the distribution of multivariate frailty model, using theory of stochastic orders, we compare multivariate frailty models arising from different choices of frailty distribution.     

Multivariate Synthetic |S| Control Chart with Variable Sampling Interval

A variable sampling interval (VSI) feature is introduced to the multivariate synthetic generalized sample variance |S| control chart. This multivariate synthetic control chart is a combination of the |S| sub-chart and the conforming run length sub-chart. The VSI feature enhances the performance of the multivariate synthetic control chart. The comparative results show that the VSI multivariate synthetic control chart performs better than other types of multivariate control charts for detecting shifts in the covariance matrix of a multivariate normally distributed process. An example is given to illustrate the operation of the VSI multivariate synthetic chart.