A BAYESIAN APPROACH FOR JOINT MODELING OF SKEW-NORMAL LONGITUDINAL MEASUREMENTS AND TIME TO EVENT DATA

• Joint modeling of longitudinal measurements and survival time has an important role in analyzing medical data sets. For example, in HIV data sets, a biological marker such as CD4 count measurements is considered as a predictor of survival. Usually, longitudinal responses of these studies are severely skew. An ordinary method for reducing the skewness is the use of square root or logarithm transformations of responses. In most of the HIV data sets, because of high rate of missingness, skewness is remained even after using the transformations. Therefore, a general form of distributions for considering skewness in the model should be used. In this paper, we have used multivariate skew-normal distribution to allow a flexible model for considering non-symmetrically of the responses. We have used a skew-normal mixed effect model for longitudinal measurements and a Cox proportional hazard model for time to event variable. These two models share some random effects. A Bayesian approach using Markov chain Monte Carlo is adopted for parameter estimation. Some simulation studies are performed to investigate the performance of the proposed method. Also, the method is illustrated using a real HIV data set. In these data, longitudinal outcomes are skew and death is considered as the event of interest. Different model structures are developed for analyzing this data set, where model selection is performed using some Bayesian criteria.

[1]  Dipankar Bandyopadhyay,et al.  Linear and Nonlinear Mixed‐Effects Models for Censored HIV Viral Loads Using Normal/Independent Distributions , 2011, Biometrics.

[2]  Yangxin Huang,et al.  Skew‐normal Bayesian nonlinear mixed‐effects models with application to AIDS studies , 2010, Statistics in medicine.

[3]  Hulin Wu,et al.  Joint inference for nonlinear mixed-effects models and time to event at the presence of missing data. , 2007, Biostatistics.

[4]  Yangxin Huang,et al.  Simultaneous Bayesian Inference for Linear, Nonlinear and Semiparametric Mixed-Effects Models with Skew-Normality and Measurement Errors in Covariates , 2011, The international journal of biostatistics.

[5]  N. Henze.  A Probabilistic Representation of the 'Skew-normal' Distribution , 1986 .

[6]  Likelihood-Based Inference for Multivariate Skew-Normal Regression Models , 2007 .

[7]  Bradley P Carlin,et al.  Separate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer Packages , 2004 .

[8]  R. Arellano-Valle,et al.  Bayesian Inference for Skew-normal Linear Mixed Models , 2007 .

[9]  R Henderson,et al.  Joint modelling of longitudinal measurements and event time data. , 2000, Biostatistics.

[10]  Amanda G Chetwynd,et al.  Joint modelling of repeated measurements and time‐to‐event outcomes: The fourth Armitage lecture , 2008, Statistics in medicine.

[11]  T. Baghfalaki,et al.  An EM Estimation Approach for Analyzing Bivariate Skew Normal Data with Non monotone Missing Values , 2011 .

[12]  A. Azzalini,et al.  The multivariate skew-normal distribution , 1996 .

[13]  F. Hsieh,et al.  Joint modelling of accelerated failure time and longitudinal data , 2005 .

[14]  A. Gelman.  Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) , 2004 .

[15]  A. Azzalini.  A class of distributions which includes the normal ones , 1985 .

[16]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[17]  R. Lempicki,et al.  General Immunologic Evaluation of Patients with Human Immunodeficiency Virus Infection , 2006 .

[18]  M. D. Branco,et al.  Bivariate random effect model using skew‐normal distribution with application to HIV‐RNA , 2007, Statistics in medicine.

[19]  M. Wulfsohn,et al.  Modeling the Relationship of Survival to Longitudinal Data Measured with Error. Applications to Survival and CD4 Counts in Patients with AIDS , 1995 .

[20]  V. H. Lachos,et al.  A nonlinear regression model with skew-normal errors , 2010 .

[21]  Hsiu J. Ho,et al.  Analysis of multivariate skew normal models with incomplete data , 2009, J. Multivar. Anal..

[22]  G. Casella,et al.  The Effect of Improper Priors on Gibbs Sampling in Hierarchical Linear Mixed Models , 1996 .

[23]  S. Sahu,et al.  A new class of multivariate skew distributions with applications to bayesian regression models , 2003 .

[24]  R. Hashemi,et al.  A Latent Process Model for Joint Modeling of Events and Marker , 2003, Lifetime data analysis.

[25]  Bradley P. Carlin,et al.  BAYES AND EMPIRICAL BAYES METHODS FOR DATA ANALYSIS , 1996, Stat. Comput..

[26]  Heleno Bolfarine,et al.  Skew normal measurement error models , 2005 .

[27]  Yangxin Huang,et al.  A bayesian approach to joint mixed-effects models with a skew-normal distribution and measurement errors in covariates. , 2011, Biometrics.

[28]  M. Wulfsohn,et al.  A joint model for survival and longitudinal data measured with error. , 1997, Biometrics.

[29]  M. Genton,et al.  On fundamental skew distributions , 2005 .

[30]  T. Baghfalaki,et al.  An ECM Estimation Approach for Analyzing Multivariate Skew-Normal Data with Dropout , 2012, Commun. Stat. Simul. Comput..

[31]  R. Arellano-Valle,et al.  LIKELIHOOD BASED INFERENCE FOR SKEW-NORMAL INDEPENDENT LINEAR MIXED MODELS , 2010 .

[32]  B P Carlin,et al.  Response of CD4 lymphocytes and clinical consequences of treatment using ddI or ddC in patients with advanced HIV infection. , 1996, Journal of acquired immune deficiency syndromes and human retrovirology : official publication of the International Retrovirology Association.

[33]  Yangxin Huang,et al.  Bayesian inference on joint models of HIV dynamics for time‐to‐event and longitudinal data with skewness and covariate measurement errors , 2011, Statistics in medicine.

[34]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[35]  D. Dey,et al.  A General Class of Multivariate Skew-Elliptical Distributions , 2001 .

[36]  A. Azzalini,et al.  Statistical applications of the multivariate skew normal distribution , 2009, 0911.2093.

[37]  Rolando De la Cruz.  Bayesian ideas and data analysis: An introduction for scientists and statisticians , 2012 .

[38]  Tsung-I Lin,et al.  Estimation and prediction in linear mixed models with skew‐normal random effects for longitudinal data , 2008, Statistics in medicine.

[39]  Inês Sousa,et al.  A review on joint modelling of longitudinal measurements and time-to-event , 2011 .

[40]  Dimitris Rizopoulos,et al.  JM: An R package for the joint modelling of longitudinal and time-to-event data , 2010 .

[41]  Anastasios A. Tsiatis,et al.  A semiparametric estimator for the proportional hazards model with longitudinal covariates measured with error , 2001 .