The bayesian bias correction method of the first-order approximation of nonlinear mixed-effects models for population...

Population pharmacokinetic analysis usually employs nonlinear mixed-effects models. To estimate the parameters, Beal and Sheiner (1982) proposed the first-order method that employs a first-order Taylor series expansion around the means of random individual parameters. Because of the small computational burden and the high convergence proportion of maximization of the log likelihood function, this method is often used in practice. However, it is known that the estimates are biased. This paper proposes a simple procedure to reduce the bias. The proposed method maximizes the nonapproximated log likelihood functions of each individual given estimates of the population parameters derived from the first-order method, and the derived Bayes estimates of the random individual parameters are utilized to improve the estimates of the population mean parameters. We confirmed that the proposed method reduced the bias using simulated data and actual erythropoietin concentration data.