The equations on page 162 should be corrected.

The analytical solution for multi-compartment models with a non-zero initial condition is complex because of the inter-compartmental transfer. An elegant solution and its implementation in the ‘wnl' R package can be useful in solving examples of textbooks and developing software of therapeutic drug monitoring, pharmacokinetic simulation, and parameter estimation. This solution uses Laplace transformation, convolution, matrix inversion, and the fact that the general solution of an inhomogeneous ordinary differential equation is the sum of a homogenous and a particular solution, together.
Drug Monitoring

There are some errors in the published article. The authors would like to make corrections in the original version of the article.

The first-order conditional estimation (FOCE) method is more complex than the first-order (FO) approximation method because it estimates the empirical Bayes estimate (EBE) for each iteration. By contrast, it is a further approximation of the Laplacian (LAPL) method, which uses second-order expansion terms. FOCE without INTERACTION can only be used for an additive error model, while FOCE with INTERACTION (FOCEI) can be used for any error model. The formula for FOCE without INTERACTION can be derived directly from the extension of the FO method, while the FOCE with INTERACTION method is a slight simplification of the LAPL method. Detailed formulas and R scripts are presented here for the reproduction of objective function values by NONMEM.
Bays ; Methods ; Reproduction*

SAS® is commonly used for bioequivalence (BE) data analysis. R is a free and open software for general purpose data analysis, and is less frequently used than SAS® for BE data analysis. This tutorial explains how R can be used for BE data analysis to generate comparable results with SAS® . The main SAS® procedures for BE data analysis are PROC GLM and PROC MIXED, and the corresponding R main packages are “sasLM” and “nlme” respectively. For fixed effects only or balanced data, the SAS® PROC GLM and R “sasLM” provide good estimates; however, for a mixed-effects model with unbalanced data, the SAS® PROC MIXED and R “nlme” are better for providing estimates without bias. The SAS® and R scripts are provided for convenience.

SAS® is commonly used for bioequivalence (BE) data analysis. R is a free and open software for general purpose data analysis, and is less frequently used than SAS® for BE data analysis. This tutorial explains how R can be used for BE data analysis to generate comparable results with SAS® . The main SAS® procedures for BE data analysis are PROC GLM and PROC MIXED, and the corresponding R main packages are “sasLM” and “nlme” respectively. For fixed effects only or balanced data, the SAS® PROC GLM and R “sasLM” provide good estimates; however, for a mixed-effects model with unbalanced data, the SAS® PROC MIXED and R “nlme” are better for providing estimates without bias. The SAS® and R scripts are provided for convenience.

In healthcare situations, time-to-event (TTE) data are common outcomes. A parametric approach is often employed to handle TTE data because it is possible to easily visualize different scenarios via simulation. Not all pharmacometricians are familiar with the use of non-linear mixed effects models (NONMEMs) to deal with TTE data. Therefore, this tutorial simply explains how to analyze TTE data using NONMEM. We show how to write the code and evaluate the model. We also provide an example of a hands-on model for training.

Analysis of a 2 × 2 table for clinical data involves computing the point estimate and confidence interval for risk difference, relative risk, or odds ratio. While point estimates of these comparative parameters are uniquely defined, several statistical methods have been proposed to estimate the confidence interval for each parameter. The Miettinen-Nurminen (MN) score method is expected to be used increasingly over traditional interval estimation methods. The MN score method has not been previously implemented in R software for data with stratification. There is a need for a comprehensive software implementation of the MN score method. This article describes the implementation of the MN score method in the sasLM R software package. To demonstrate the usage of the sasLM functions introduced, recently published clinical data are provided as examples.

Wald confidence interval has been used as the conventional method of interval estimation for the parameters in nonlinear models. Because Wald confidence interval is symmetric around the point estimate, it does not reflect the asymmetry of the likelihood profile in nonlinear regression. In contrast, a likelihood interval is estimated directly from the likelihood profile and does reflect the shape of the likelihood profile. However, the lack of software for the estimation of likelihood intervals and visualization of likelihood profiles posed an obstacle to the use of likelihood intervals in nonlinear models. There was a need for software implementation to tackle these tasks. Likelihood interval estimation and likelihood profile plotting for nonlinear models had not been previously implemented in R software. This article describes the implementation of likelihood interval estimation and likelihood profile plotting in the wnl R software package. To demonstrate the usage of implemented functions, an example of fitting a nonlinear pharmacokinetic model to concentration-time data is presented.

BACKGROUND: Bivalirudin is a direct thrombin inhibitor for patients with unstable angina undergoing percutaneous coronary intervention. METHODS: A sensitive liquid chromatography-tandem mass spectrometry (LC/MS/MS) method was developed and validated for the determination of bivalirudin, in human plasma using nafarelin as internal standard (IS). Chromatographic separation was performed using a Shiseido MG3 mm column (2.0 x 50 mm) with a gradient mobile phase consisting of water and acetonitrile containing 0.1 % formic acid at a flow rate of 0.4 mL/min, and total run time was within 5 min. Detection and quantification was performed by the mass spectrometer using a multiple reaction-monitoring mode at m/z 1091.0 --> 650.3 for bivalirudin, and m/z 662.1 --> 249.3 for IS. RESULTS: The assay was linear over a concentration range of 10 - 10000 ng/mL with a lower limit of quantification of 10 ng/mL in human plasma. CONCLUSION: This method was successfully applied for pharmacokinetics study after intravenous administration of bivalirudin to healthy Korean male volunteers.
Administration, Intravenous ; Angina, Unstable ; Chromatography, Liquid ; Humans* ; Male ; Mass Spectrometry ; Methods ; Nafarelin ; Percutaneous Coronary Intervention ; Pharmacokinetics ; Plasma* ; Thrombin ; Water