Longitudinal data arise frequently in medical studies and it is a common practice to analyze such complex data with nonlinear mixed-effects (NLME) models, which enable us to account for between-subject and within-subject variations. To partially explain the variations, time-dependent covariates are usually introduced to these models. Some covariates, however, may be often measured with substantial errors and missing observations. It is often the case that model random error is assumed to be distributed normally, but the normality assumption may not always give robust and reliable results, particularly if the data exhibit skewness. In the literature, there has been considerable interest in accommodating either skewed response or covariate measured with error and missing data in such models, but there has been relatively little study concerning all these features simultaneously. This article is to address simultaneous impact of skewness in response and measurement error and missing data in covariate by jointly modeling the response and covariate processes under a framework of Bayesian semiparametric nonlinear mixed-effects models. In particular, we aim at exploring how mixed-effects joint models based on one-compartment model with one phase time-varying decay rate and two-compartment model with two phase time-varying decay rates contribute to modeling results and inference. The method is illustrated by an AIDS data example to compare potential models with different distributional specifications and various scenarios. The findings from this study suggest that the one-compartment model with a skewnormal distribution may provide more reasonable results if the data exhibit skewness in response and/or have measurement error and missing observations in covariates.