The joint modeling of mean and variance for longitudinal data is an active research area. This type of model has the advantage of accounting for heteroscedasticity commonly observed in between and within subject variations. Most of researches focus on improving the estimating efficiency but ignore many data features frequently encountered in practice. In this article, we develop a mixed-effects location scale joint model that concurrently accounts for longitudinal data with multiple features. Specifically, our joint model handles heterogeneity, skewness, limit of detection, measurement errors in covariates which are typically observed in the collection of longitudinal data from many studies. We employ a Bayesian approach for making inference on the joint model. The proposed model and method are applied to an AIDS study. Simulation studies are performed to assess the performance of the proposed method. Alternative models under different conditions are compared.
PourahmadiM. Joint mean-covariance models with applications to longitudinal data: unconstrained parameterisation. Biometrika1999; 86: 677–690.
2.
ElliottM. Identifying latent clusters of variability in longitudinal data. Biostatistics2007; 8: 756–771.
3.
HedekerDMermelsteinRJDemirtasH. An application of a mixed-effects location scale model for analysis of ecological momentary assessment (EMA) data. Biometrics2008; 64: 627–634.
4.
LiXHedekerD. A three-level mixed-effects location scale model with an application to ecological momentary assessment data. Stat Med2012; 31: 3192–3210.
5.
YaoWLiR. New local estimation procedure for a non-parametric regression function for longitudinal data. J R Stat Soc Ser B Stat Methodol2013; 75: 123–138.
6.
ZhangWLengCTangCY. A joint modelling approach for longitudinal studies. J R Stat Soc Ser B Stat Methodol2015; 77: 219–238.
7.
KapurKLiXBloodEAet al.Bayesian mixed-effects location and scale models for multivariate longitudinal outcomes: an application to ecological momentary assessment data. Stat Med2015; 34: 630–651.
8.
SahuSDeyDBrancoM. A new class of multivariate skew distributions with applications to Bayesian regression models. Can J Stat2003; 31: 129–150.
9.
LachosVBandyopadhyayDDeyD. Linear and nonlinear mixed-effects models for censored HIV viral loads using normal/independent distributions. Biometrics2011; 67: 1594–1604.
10.
HoHLinT. Robust linear mixed models using the skew-t distribution with application to schizophrenia data. Biom J2010; 52: 449–469.
11.
WuL. A joint model for nonlinear mixed-effects models with censoring and covariates measured with error. J Am Stat Assoc2002; 97: 955–964.
12.
HughesJ. Mixed effects models with censored data with application to HIV RNA levels. Biometrics1999; 55: 625–629.
13.
LiuWWuL. Simultaneous inference for semiparametric nonlinear mixed-effects models with covariate measurement errors and missing responses. Biometrics2007; 63: 342–350.
14.
HuangYDagneG. A Bayesian approach to joint mixed-effects models with a skew-normal distribution and measurement errors in covariates. Biometrics2011; 67: 260–269.
15.
HuangYDagneGWuL. Mixed-effects joint models with skew-normal distribution for HIV dynamic response with missing and mismeasured time-varying covariate. Int J Biostat2012; 8: 34–34.
16.
LuTHuangY. Bayesian inference on mixed-effects varying-coefficient joint models with skew-t distribution for longitudinal data with multiple features. Stat Methods Med Res. 2015; 26: 1146–1164.
17.
LuXHuangYZhuY. Finite mixture of nonlinear mixed-effects joint models in the presence of missing and mismeasured covariate, with application to AIDS studies. Comput Stat Data Anal2016; 93: 119–130.
18.
AzzaliniACapitanioA. Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t-distribution. J R Stat Soc Ser B Stat Methodol2003; 65: 367–389.
19.
Arellano-ValleRGentonM. On fundamental skew distributions. J Multivariate Anal2005; 96: 93–116.
20.
AzzaliniAGentonM. Robust likelihood methods based on the skew-t and related distributions. Int Stat Rev2008; 76: 106–129.
21.
JaraAQuintanaFMartinE. Linear mixed models with skew-elliptical distributions: a Bayesian approach. Comput Stat Data Anal2008; 52: 5033–5045.
22.
LunnDThomasABestNet al.WinBUGS – a Bayesian modelling framework: concepts, structure, an extensibility. Stat Comput2000; 10: 325–337.
23.
FonsecaTFerreiraMMigonH. Objective Bayesian analysis for the Student-t regression model. Biometrika2008; 95: 325–333.
24.
VillaCWalkerS. Objective prior for the number of degrees of freedom of a t distribution. Bayesian Anal2014; 9: 197–220.
25.
GelmanA. Prior distributions for variance parameters in hierarchical models. Bayesian Anal2006; 3: 515–533.
26.
GelmanARubinD. Inference from iterative simulation using multiple sequences. Stat Sci1992; 7: 457–511.
27.
SpiegelhalterDBestNCarlinBet al.Bayesian measures of model complexity and fit. J R Stat Soc Ser B Stat Methodol2002; 64: 583–639.
28.
GelmanACarlinJSternHet al.Bayesian data analysis, London: Chapman and Hall, 2003.
