Permutation and Bayesian tests for testing random effects in linear mixed-effects models

Smoothed simulated pseudo-maximum likelihood estimation for nonlinear mixed effects models with censored responses

Nonlinear mixed effects models have been widely applied to analyses of data that arise from biological, agricultural, and environmental sciences. Estimation of and inference on parameters in nonlinear mixed effects models are often based on the specification of a likelihood function. Maximizing this likelihood function can be complicated by the specification of the random effects distribution, especially in the presence of multiple random effects. The implementation of nonlinear mixed effects models can be further complicated by left-censored responses, representing measurements from bioassays where the exact quantification below a certain threshold is not possible. Motivated by the need to characterize the nonlinear human immunodeficiency virus RNA viral load trajectories after the interruption of antiretroviral therapy, we propose a smoothed simulated pseudo-maximum likelihood estimation approach to fit nonlinear mixed effects models in the presence of left-censored observations. We establish the consistency and asymptotic normality of the resulting estimators. We develop testing procedures for the correlation among random effects and for testing the distributional assumptions on random effects against a specific alternative. In contrast to the existing variants of expectation-maximization approaches, the proposed methods offer flexibility in the specification of the random effects distribution and convenience in making inference about higher-order correlation parameters. We evaluate the finite-sample performance of the proposed methods through extensive simulation studies and illustrate them on a combined dataset from six AIDS Clinical Trials Group treatment interruption studies.

MEGH: A parametric class of general hazard models for clustered survival data

In many applications of survival data analysis, the individuals are treated in different medical centres or belong to different clusters defined by geographical or administrative regions. The analysis of such data requires accounting for between-cluster variability. Ignoring such variability would impose unrealistic assumptions in the analysis and could affect the inference on the statistical models. We develop a novel parametric mixed-effects general hazard (MEGH) model that is particularly suitable for the analysis of clustered survival data. The proposed structure generalises the mixed-effects proportional hazards and mixed-effects accelerated failure time structures, among other structures, which are obtained as special cases of the MEGH structure. We develop a likelihood-based algorithm for parameter estimation in general subclasses of the MEGH model, which is implemented in our R package MEGH. We propose diagnostic tools for assessing the random effects and their distributional assumption in the proposed MEGH model. We investigate the performance of the MEGH model using theoretical and simulation studies, as well as a real data application on leukaemia.