Influence analysis for skew‐normal semiparametric joint models of multivariate longitudinal and multivariate survival data

The normality assumption of measurement error is a widely used distribution in joint models of longitudinal and survival data, but it may lead to unreasonable or even misleading results when longitudinal data reveal skewness feature. This paper proposes a new joint model for multivariate longitudinal and multivariate survival data by incorporating a nonparametric function into the trajectory function and hazard function and assuming that measurement errors in longitudinal measurement models follow a skew-normal distribution. A Monte Carlo Expectation-Maximization (EM) algorithm together with the penalized-splines technique and the Metropolis-Hastings algorithm within the Gibbs sampler is developed to estimate parameters and nonparametric functions in the considered joint models. Case deletion diagnostic measures are proposed to identify the potential influential observations, and an extended local influence method is presented to assess local influence of minor perturbations. Simulation studies and a real example from a clinical trial are presented to illustrate the proposed methodologies. Copyright © 2017 John Wiley & Sons, Ltd.

[1]  R. Dennis Cook,et al.  Detection of Influential Observation in Linear Regression , 2000, Technometrics.

[2]  Hongtu Zhu,et al.  Diagnostic measures for empirical likelihood of general estimating equations , 2008 .

[3]  E Lesaffre,et al.  Local influence in linear mixed models. , 1998, Biometrics.

[4]  Geoffrey J. McLachlan,et al.  On mixtures of skew normal and skew $$t$$-distributions , 2012, Adv. Data Anal. Classif..

[5]  Joseph G Ibrahim,et al.  Joint Models for Multivariate Longitudinal and Multivariate Survival Data , 2006, Biometrics.

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

[7]  G. C. Wei,et al.  A Monte Carlo Implementation of the EM Algorithm and the Poor Man's Data Augmentation Algorithms , 1990 .

[8]  Bo-Cheng Wei,et al.  Case-deletion measures for models with incomplete data , 2001 .

[9]  Joseph G. Ibrahim,et al.  BAYESIAN METHODS FOR JOINT MODELING OF LONGITUDINAL AND SURVIVAL DATA WITH APPLICATIONS TO CANCER VACCINE TRIALS , 2004 .

[10]  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.

[11]  Sik-Yum Lee,et al.  Local influence analysis of nonlinear structural equation models , 2004 .

[12]  Yangxin Huang,et al.  Bayesian Semiparametric Nonlinear Mixed‐Effects Joint Models for Data with Skewness, Missing Responses, and Measurement Errors in Covariates , 2012, Biometrics.

[13]  Wing K. Fung,et al.  Influence diagnostics and outlier tests for semiparametric mixed models , 2002 .

[14]  Yat Sun Poon,et al.  Conformal normal curvature and assessment of local influence , 1999 .

[15]  Paul H. C. Eilers,et al.  Flexible smoothing with B-splines and penalties , 1996 .

[16]  S. Lang,et al.  Bayesian P-Splines , 2004 .

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

[18]  S. Weisberg,et al.  Residuals and Influence in Regression , 1982 .

[19]  J. Ibrahim,et al.  A Flexible B‐Spline Model for Multiple Longitudinal Biomarkers and Survival , 2005, Biometrics.

[20]  M. Genton,et al.  Multivariate extended skew-t distributions and related families , 2010 .

[21]  R. Cook.  Assessment of Local Influence , 1986 .

[22]  Sik-Yum Lee,et al.  Local influence for incomplete data models , 2001 .

[23]  Hongtu Zhu,et al.  Bayesian Influence Measures for Joint Models for Longitudinal and Survival Data , 2012, Biometrics.

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