Jointly Modeling Event Time and Skewed-Longitudinal Data with Missing Response and Mismeasured Covariate for AIDS Studies

In longitudinal studies it is often of interest to investigate how a repeatedly measured marker in time is associated with a time to an event of interest. This type of research question has given rise to a rapidly developing field of biostatistics research that deals with the joint modeling of longitudinal and time-to-event data. Normality of model errors in longitudinal model is a routine assumption, but it may be unrealistically obscuring important features of subject variations. Covariates are usually introduced in the models to partially explain between- and within-subject variations, but some covariates such as CD4 cell count may be often measured with substantial errors. Moreover, the responses may encounter nonignorable missing. Statistical analysis may be complicated dramatically based on longitudinal-survival joint models where longitudinal data with skewness, missing values, and measurement errors are observed. In this article, we relax the distributional assumptions for the longitudinal models using skewed (parametric) distribution and unspecified (nonparametric) distribution placed by a Dirichlet process prior, and address the simultaneous influence of skewness, missingness, covariate measurement error, and time-to-event process by jointly modeling three components (response process with missing values, covariate process with measurement errors, and time-to-event process) linked through the random-effects that characterize the underlying individual-specific longitudinal processes in Bayesian analysis. The method is illustrated with an AIDS study by jointly modeling HIV/CD4 dynamics and time to viral rebound in comparison with potential models with various scenarios and different distributional specifications.

[1]  Lang Wu,et al.  Simultaneous inference and bias analysis for longitudinal data with covariate measurement error and missing responses. , 2011, Biometrics.

[2]  D. Rubin,et al.  Statistical Analysis with Missing Data. , 1989 .

[3]  Ernesto San Martín,et al.  Linear mixed models with skew-elliptical distributions: A Bayesian approach , 2008, Comput. Stat. Data Anal..

[4]  Hulin Wu,et al.  Identification of significant host factors for HIV dynamics modelled by non‐linear mixed‐effects models , 2002, Statistics in medicine.

[5]  T. Ferguson.  A Bayesian Analysis of Some Nonparametric Problems , 1973 .

[6]  Raymond J. Carroll,et al.  Measurement error in nonlinear models: a modern perspective , 2006 .

[7]  D. Dunson.  Bayesian dynamic modeling of latent trait distributions. , 2006, Biostatistics.

[8]  S. Lipsitz,et al.  Missing responses in generalised linear mixed models when the missing data mechanism is nonignorable , 2001 .

[9]  Roderick J. A. Little,et al.  Statistical Analysis with Missing Data: Little/Statistical Analysis with Missing Data , 2002 .

[10]  A A Ding,et al.  Relationships between antiviral treatment effects and biphasic viral decay rates in modeling HIV dynamics. , 1999, Mathematical biosciences.

[11]  W Liu,et al.  Joint Inference on HIV Viral Dynamics and Immune Suppression in Presence of Measurement Errors , 2010, Biometrics.

[12]  John K Kruschke,et al.  Bayesian data analysis. , 2010, Wiley interdisciplinary reviews. Cognitive science.

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

[14]  Andrew Thomas,et al.  WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility , 2000, Stat. Comput..

[15]  J. Ibrahim,et al.  A Bayesian semiparametric joint hierarchical model for longitudinal and survival data. , 2003, Biometrics.

[16]  A. Ding,et al.  Assessing antiviral potency of anti-HIV therapies in vivo by comparing viral decay rates in viral dynamic models. , 2001, Biostatistics.

[17]  Lang Wu,et al.  A Joint Model for Nonlinear Mixed-Effects Models With Censoring and Covariates Measured With Error, With Application to AIDS Studies , 2002 .

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

[19]  Hulin Wu,et al.  Hierarchical Bayesian Methods for Estimation of Parameters in a Longitudinal HIV Dynamic System , 2006, Biometrics.

[20]  M. Escobar.  Estimating Normal Means with a Dirichlet Process Prior , 1994 .

[21]  Tsung-I Lin,et al.  Robust linear mixed models using the skew t distribution with application to schizophrenia data , 2010, Biometrical journal. Biometrische Zeitschrift.

[22]  J G Ibrahim,et al.  A semiparametric Bayesian approach to the random effects model. , 1998, Biometrics.

[23]  Hulin Wu,et al.  The study of long‐term HIV dynamics using semi‐parametric non‐linear mixed‐effects models , 2002, Statistics in medicine.

[24]  M. Lederman,et al.  Immunologic responses associated with 12 weeks of combination antiretroviral therapy consisting of zidovudine, lamivudine, and ritonavir: results of AIDS Clinical Trials Group Protocol 315. , 1998, The Journal of infectious diseases.

[25]  D Commenges,et al.  Maximum Likelihood Estimation in Dynamical Models of HIV , 2007, Biometrics.

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

[27]  Wei Liu,et al.  Simultaneous Inference for Semiparametric Nonlinear Mixed‐Effects Models with Covariate Measurement Errors and Missing Responses , 2007, Biometrics.

[28]  A. Azzalini,et al.  Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t‐distribution , 2003, 0911.2342.

[29]  A. Perelson,et al.  HIV-1 Dynamics in Vivo: Virion Clearance Rate, Infected Cell Life-Span, and Viral Generation Time , 1996, Science.

[30]  Nian-Sheng Tang,et al.  A semiparametric Bayesian approach to generalized partial linear mixed models for longitudinal data , 2012, Comput. Stat. Data Anal..

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

[32]  V. De Gruttola,et al.  Modelling progression of CD4-lymphocyte count and its relationship to survival time. , 1994, Biometrics.

[33]  Lancelot F. James,et al.  Approximate Dirichlet Process Computing in Finite Normal Mixtures , 2002 .

[34]  H Wu,et al.  Population HIV‐1 Dynamics In Vivo: Applicable Models and Inferential Tools for Virological Data from AIDS Clinical Trials , 1999, Biometrics.

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

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