Bayesian quantile regression-based nonlinear mixed-effects joint models for time-to-event and longitudinal data with multiple features

Semiparametric multivariate joint model for skewed-longitudinal and survival data: A Bayesian approach

Joint models and statistical inference for longitudinal and survival data have been an active area of statistical research and have mostly coupled a longitudinal biomarker-based mixed-effects model with normal distribution and an event time-based survival model. In practice, however, the following issues may standout: (i) Normality of model error in longitudinal models is a routine assumption, but it may be unrealistically violating data features of subject variations. (ii) Data collected are often featured by the mixed types of multiple longitudinal outcomes which are significantly correlated, ignoring their correlation may lead to biased estimation. Additionally, a parametric model specification may be inflexible to capture the complicated patterns of longitudinal data. (iii) Missing observations in the longitudinal data are often encountered; the missing measures are likely to be informative (nonignorable) and ignoring this phenomenon may result in inaccurate inference. Multilevel item response theory (MLIRT) models have been increasingly used to analyze the multiple longitudinal data of mixed types (ie, continuous and categorical) in clinical studies. In this article, we develop an MLIRT-based semiparametric joint model with skew-t distribution that consists of an extended MLIRT model for the mixed types of multiple longitudinal data and a Cox proportional hazards model, linked through random-effects. A Bayesian approach is employed for joint modeling. Simulation studies are conducted to assess performance of the proposed models and method. A real example from primary biliary cirrhosis clinical study is analyzed to estimate parameters in the joint model and also evaluate sensitivity of parameter estimates for various plausible nonignorable missing data mechanisms.

Study on Bayesian Skew-Normal Linear Mixed Model and Its Application in Fire Insurance

Fire insurance is a crucial component of property insurance, and its rating depends on the forecast of insurance loss claim data. Fire insurance loss claim data have complicated characteristics such as skewness and heavy tail. The traditional linear mixed model is commonly difficult to accurately describe the distribution of loss. Therefore, it is crucial to establish a scientific and reasonable distribution model of fire insurance loss claim data. In this study, the random effects and random errors in the linear mixed model are firstly assumed to obey the skew-normal distribution. Then, a skew-normal linear mixed model is established using the Bayesian MCMC method based on a set of U.S. property insurance loss claims data. Comparative analysis is conducted with the linear mixed model of logarithmic transformation. Afterward, a Bayesian skew-normal linear mixed model for Chinese fire insurance loss claims data is designed. The posterior distribution of claim data parameters and related parameter estimation are employed with the R language JAGS package to obtain the predicted and simulated loss claim values. Finally, the optimization model in this study is used to determine the insurance rate. The results demonstrate that the model established by the Bayesian MCMC method can overcome data skewness, and the fitting and correlation with the sample data are better than the log-normal linear mixed model. Hence, it can be concluded that the distribution model proposed in this paper is reasonable for describing insurance claims. This study innovates a new approach for calculating the insurance premium rate and expands the application of the Bayesian method in the fire insurance field.

Bayesian Joint Modeling of Multivariate Longitudinal and Survival Data With an Application to Diabetes Study

Joint models of longitudinal and time-to-event data have received a lot of attention in epidemiological and clinical research under a linear mixed-effects model with the normal assumption for a single longitudinal outcome and Cox proportional hazards model. However, those model-based analyses may not provide robust inference when longitudinal measurements exhibit skewness and/or heavy tails. In addition, the data collected are often featured by multivariate longitudinal outcomes which are significantly correlated, and ignoring their correlation may lead to biased estimation. Under the umbrella of Bayesian inference, this article introduces multivariate joint (MVJ) models with a skewed distribution for multiple longitudinal exposures in an attempt to cope with correlated multiple longitudinal outcomes, adjust departures from normality, and tailor linkage in specifying a time-to-event process. We develop a Bayesian joint modeling approach to MVJ models that couples a multivariate linear mixed-effects (MLME) model with the skew-normal (SN) distribution and a Cox proportional hazards model. Our proposed models and method are evaluated by simulation studies and are applied to a real example from a diabetes study.

Modelling the association between biomarkers and clinical outcome: an introduction to nonlinear joint models

Nonlinear joint models are a powerful tool to precisely analyze the association between a nonlinear biomarker and a time-to-event process, such as death. Here, we review the main methodological techniques required to build these models and to make inferences and predictions. We describe the main clinical applications and discuss the future developments of such models.

Bayesian joint modeling for partially linear mixed-effects quantile regression of longitudinal and time-to-event data with limit of detection, covariate measurement errors and skewness

Joint modeling analysis of longitudinal and time-to-event data has been an active area of statistical methodological study and biomedical research, but the majority of them are based on mean-regression. Quantile regression (QR) can characterize the entire conditional distribution of the outcome variable, and may be more robust to outliers/heavy tails and misspecification of error distribution. Additionally, a parametric specification may be insufficient and inflexible to capture the complicated longitudinal pattern of biomarkers. Thus, this study proposes novel QR-based partially linear mixed-effects joint models with three components (QR-based longitudinal response, longitudinal covariate, and time-to-event processes), and applies to Multicenter AIDS Cohort Study (MACS). Many common data features, including left-censoring due to a limit of detection, covariate measurement error, and asymmetric distribution, are considered to obtain reliable parameter estimates. Many interesting findings are discovered by the complicated joint models under Bayesian inference framework. Simulation studies are also implemented to assess the performance of the proposed joint models under different scenarios.

Bayesian joint modelling of longitudinal and time to event data: a methodological review

BACKGROUND

In clinical research, there is an increasing interest in joint modelling of longitudinal and time-to-event data, since it reduces bias in parameter estimation and increases the efficiency of statistical inference. Inference and prediction from frequentist approaches of joint models have been extensively reviewed, and due to the recent popularity of data-driven Bayesian approaches, a review on current Bayesian estimation of joint model is useful to draw recommendations for future researches.

METHODS

We have undertaken a comprehensive review on Bayesian univariate and multivariate joint models. We focused on type of outcomes, model assumptions, association structure, estimation algorithm, dynamic prediction and software implementation.

RESULTS

A total of 89 articles have been identified, consisting of 75 methodological and 14 applied articles. The most common approach to model the longitudinal and time-to-event outcomes jointly included linear mixed effect models with proportional hazards. A random effect association structure was generally used for linking the two sub-models. Markov Chain Monte Carlo (MCMC) algorithms were commonly used (93% articles) to estimate the model parameters. Only six articles were primarily focused on dynamic predictions for longitudinal or event-time outcomes.

CONCLUSION

Methodologies for a wide variety of data types have been proposed; however the research is limited if the association between the two outcomes changes over time, and there is also lack of methods to determine the association structure in the absence of clinical background knowledge. Joint modelling has been proved to be beneficial in producing more accurate dynamic prediction; however, there is a lack of sufficient tools to validate the prediction.

Quantile regression-based Bayesian joint modeling analysis of longitudinal-survival data, with application to an AIDS cohort study

In longitudinal studies, it is of interest to investigate how repeatedly measured markers are associated with time to an event. Joint models have received increasing attention on analyzing such complex longitudinal-survival data with multiple data features, but most of them are mean regression-based models. This paper formulates a quantile regression (QR) based joint models in general forms that consider left-censoring due to the limit of detection, covariates with measurement errors and skewness. The joint models consist of three components: (i) QR-based nonlinear mixed-effects Tobit model using asymmetric Laplace distribution for response dynamic process; (ii) nonparametric linear mixed-effects model with skew-normal distribution for mismeasured covariate; and (iii) Cox proportional hazard model for event time. For the purpose of simultaneously estimating model parameters, we propose a Bayesian method to jointly model the three components which are linked through the random effects. We apply the proposed modeling procedure to analyze the Multicenter AIDS Cohort Study data, and assess the performance of the proposed models and method through simulation studies. The findings suggest that our QR-based joint models may provide comprehensive understanding of heterogeneous outcome trajectories at different quantiles, and more reliable and robust results if the data exhibits these features.

Modelling and estimation of nonlinear quantile regression with clustered data

In regression applications, the presence of nonlinearity and correlation among observations offer computational challenges not only in traditional settings such as least squares regression, but also (and especially) when the objective function is nonsmooth as in the case of quantile regression. Methods are developed for the modelling and estimation of nonlinear conditional quantile functions when data are clustered within two-level nested designs. The proposed estimation algorithm is a blend of a smoothing algorithm for quantile regression and a second order Laplacian approximation for nonlinear mixed models. This optimization approach has the appealing advantage of reducing the original nonsmooth problem to an approximated L 2 problem. While the estimation algorithm is iterative, the objective function to be optimized has a simple analytic form. The proposed methods are assessed through a simulation study and two applications, one in pharmacokinetics and one related to growth curve modelling in agriculture.

Bayesian quantile regression-based partially linear mixed-effects joint models for longitudinal data with multiple features

In longitudinal AIDS studies, it is of interest to investigate the relationship between HIV viral load and CD4 cell counts, as well as the complicated time effect. Most of common models to analyze such complex longitudinal data are based on mean-regression, which fails to provide efficient estimates due to outliers and/or heavy tails. Quantile regression-based partially linear mixed-effects models, a special case of semiparametric models enjoying benefits of both parametric and nonparametric models, have the flexibility to monitor the viral dynamics nonparametrically and detect the varying CD4 effects parametrically at different quantiles of viral load. Meanwhile, it is critical to consider various data features of repeated measurements, including left-censoring due to a limit of detection, covariate measurement error, and asymmetric distribution. In this research, we first establish a Bayesian joint models that accounts for all these data features simultaneously in the framework of quantile regression-based partially linear mixed-effects models. The proposed models are applied to analyze the Multicenter AIDS Cohort Study (MACS) data. Simulation studies are also conducted to assess the performance of the proposed methods under different scenarios.