Handling Data Below the Limit of Quantification in Mixed Effect Models

The purpose of this study is to investigate the impact of observations below the limit of quantification (BQL) occurring in three distinctly different ways and assess the best method for prevention of bias in parameter estimates and for illustrating model fit using visual predictive checks (VPCs). Three typical ways in which BQL can occur in a model was investigated with simulations from three different models and different levels of the limit of quantification (LOQ). Model A was used to represent a case with BQL observations in an absorption phase of a PK model whereas model B represented a case with BQL observations in the elimination phase. The third model, C, an indirect response model illustrated a case where the variable of interest in some cases decreases below the LOQ before returning towards baseline. Different approaches for handling of BQL data were compared with estimation of the full dataset for 100 simulated datasets following models A, B, and C. An improved standard for VPCs was suggested to better evaluate simulation properties both for data above and below LOQ. Omission of BQL data was associated with substantial bias in parameter estimates for all tested models even for seemingly small amounts of censored data. Best performance was seen when the likelihood of being below LOQ was incorporated into the model. In the tested examples this method generated overall unbiased parameter estimates. Results following substitution of BQL observations with LOQ/2 were in some cases shown to introduce bias and were always suboptimal to the best method. The new standard VPCs was found to identify model misfit more clearly than VPCs of data above LOQ only.

[1]  Mats O. Karlsson,et al.  Impact of Omission or Replacement of Data Below the Limit of Quantification on Parameter Estimates in a Two-Compartment Model , 2002, Pharmaceutical Research.

[2]  Mats O. Karlsson,et al.  Population Pharmacokinetics of Rifampin in Pulmonary Tuberculosis Patients, Including a Semimechanistic Model To Describe Variable Absorption , 2008, Antimicrobial Agents and Chemotherapy.

[3]  Stuart L. Beal,et al.  Ways to Fit a PK Model with Some Data Below the Quantification Limit , 2001, Journal of Pharmacokinetics and Pharmacodynamics.

[4]  E. Niclas Jonsson,et al.  Erratum to "PsN-Toolkit - A collection of computer intensive statistical methods for non-linear mixed effect modeling using NONMEM" [Comput. Methods Prog. Biomedicine 79 (2005) 241-257] , 2005, Comput. Methods Programs Biomed..

[5]  Mats O. Karlsson,et al.  Likelihood based approaches to handling data below the quantification limit using NONMEM VI , 2008, Journal of Pharmacokinetics and Pharmacodynamics.

[6]  E. Niclas Jonsson,et al.  Perl-speaks-NONMEM (PsN) - a Perl module for NONMEM related programming , 2004, Comput. Methods Programs Biomed..

[7]  Courtney V. Fletcher,et al.  Impact of censoring data below an arbitrary quantification limit on structural model misspecification , 2008, Journal of Pharmacokinetics and Pharmacodynamics.

[8]  Varun Garg,et al.  Comparison of four basic models of indirect pharmacodynamic responses , 1993, Journal of Pharmacokinetics and Biopharmaceutics.

[9]  E. Niclas Jonsson,et al.  PsN-Toolkit - A collection of computer intensive statistical methods for non-linear mixed effect modeling using NONMEM , 2005, Comput. Methods Programs Biomed..

[10]  Steven G. Woolfrey,et al.  Analysis of Toxicokinetic Data Using NONMEM: Impact of Quantification Limit and Replacement Strategies for Censored Data , 2001, Journal of Pharmacokinetics and Pharmacodynamics.

[11]  Nicholas H. G. Holford,et al.  The Visual Predictive Check Superiority to Standard Diagnostic (Rorschach) Plots , 2005 .

[12]  E N Jonsson,et al.  Xpose--an S-PLUS based population pharmacokinetic/pharmacodynamic model building aid for NONMEM. , 1999, Computer methods and programs in biomedicine.