Clinical Chemistry
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Clinical Chemistry 52: 1606-1608, 2006; 10.1373/clinchem.2006.069476
This Article
Right arrow Extract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Holmes, D. T.
Right arrow Articles by Buhr, K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Holmes, D. T.
Right arrow Articles by Buhr, K.
Related Collections
Right arrow Evidence Based Laboratory Medicine and Test Utilization
(Clinical Chemistry. 2006;52:1606-1608.)
© 2006 American Association for Clinical Chemistry, Inc.


Letters to the Editor

Mathematical Modeling: Assumptions Affect Results

Daniel T. Holmes1,a and Kevin Buhr2

1 St. Paul’s Hospital, Department of Laboratory Medicine, and the Department of Pathology, University of British Columbia, Vancouver, BC, Canada
2 Department of Mathematics, University of British Columbia, Vancouver, BC, Canada

aAddress correspondence to this author at: St. Paul’s Hospital Department of Laboratory Medicine and the University of British Columbia Department of Pathology. 1081 Burrard Street, Vancouver, BC, Canada, V6Z 1Y6. Fax 604-806h8815; e-mail dtholmes{at}interchange.ubc.ca.


To the Editor:

Apple et al. (1) recently presented an analysis of the effect of imprecision (CV) on false positive cardiac troponin I (cTnI) results. After examining their methodology, we have identified problems with the model and arguments.

Key mathematical details are omitted, and some assumptions appear unrealistic. The authors model the plasma cTnI concentrations in a healthy reference population by an exponential distribution and the analytical error as gaussian using 2 assumed CV profiles: a baseline profile and the same scaled by a factor of 0.40. We observe 2 difficulties.

First, the baseline profile is incompletely specified, and the points provided (CVs of 37.5%, 25.0%, and 9.4% at true cTnIs of 0.05, 0.07, and 0.14 µg/L) imply SDs that decrease with increasing cTnI, which does not reflect reality (2)(3). Further, the "SD of the blank" (cTnI = 0 µg/L) is entirely unspecified. To contrast, we constructed empirical profiles using mean concentration data and CVs from a recent paper comparing 15 cTnI assays (4). Fig. 1a compares the 5 cTnI assays whose data overlap the concentration range 0.05–0.14 µg/L to those used by Apple et al. (1). For the experimentally derived profiles, the SD function has both a positive slope and a "y-intercept", features assumed in theoretical models (5). Unfortunately the authors’ profile biases the results in favor of their conclusions.


Figure 1
View larger version (19K):
[in this window]
[in a new window]
 
Figure 1. (A), Experimentally determined precision profiles for several troponin assays, reconstructed from Panteghini et al. (4), are plotted as standard deviation vs mean concentration and contrasted with hypothetical precision profiles of Apple et al. (1).

(B), Hypothetical continuous probability density functions (PDFs) for a single determination of true or measured cTnI as generated using the assumptions of Apple et al. (1). As a consequence of continuous PDFs, the probability of making a measurement between 2 concentrations of cTnI is given by

Formula
where z1 is the lower cTnI concentration and z2 is the higher cTnI concentration. Equation for the true value PDF is given by p(x) = {lambda}exp(–{lambda}x), where {lambda}{approx}88.6 L/µg, and x is the concentration of cTnI in µg/L. The 2 measured value PDFs were generated using the convolution product (7) of the exponential distribution function of true values with a distribution of gaussian form whose SD is given by {varsigma}(x)=–0.0625x+0.021875 ("Baseline"), or the same profile scaled by exactly 0.4 ("Scaled"— i.e., the precision profiles described by Apple et al.).

Second, simulation with gaussian measurement error gives negative concentrations. If infrequent, such values might be dropped or taken as zero, but our attempts to duplicate the authors’ analysis give many negative values. Specifically, noting that the 3 supplied profile points are perfectly collinear when plotted as SD vs mean (in itself unrealistic), and assuming this collinearity holds everywhere, the resulting CV profile gives ~160 000 (32%) negative measured values in a run of 500 000 simulated values (Fig. 1bUp ). Fig. 1Up of the technical brief by Apple et al. does not lend clarity as it apparently shows true, not measured, values. How did the authors handle negative simulated values? Did they discard them, possibly biasing their assumed reference population toward higher cTnI concentrations?

Yet, if we agree that a suitable modification of their model might lead to results such as those presented, albeit not numerically identical, we fail to see how they lead to the conclusions drawn. Specifically, the authors first intend to show that, if 2 assays are identical except that one CV profile is 0.4 times that of the other, for 3 sequential measurements on 1000 healthy individuals, the more precise assay will result in 3–5 fewer patients exceeding the 99th percentile in at least one measurement. On this basis, the authors suggest that the analysis is relevant to clinical practice, as when patients present with chest pain. Second, they conclude that only the 99th percentile cutoff (and not the concentration at CV = 10%) should be used.

With respect to the first point, data collected at our institution indicate that emergency room patients have a very different cTnI distribution than healthy individuals, with a larger fraction of patients near the 99th percentile cutoff. Attempts to assess the error rates of serial measurements by use of a model based on healthy people will therefore have little bearing on real clinical practice. To appropriately compare more and less precise assays, with or without serial measurements, the model must simulate cTnI distributions appropriate to an emergency room setting; additionally, it must consider both false positives and false negatives.

With respect to the second point, we feel the results simply do not address it. For the 2 hypothetical assays, the CV = 10% concentrations of 0.07 µg/L and 0.14 µg/L exceed the respective 99th percentile cutoffs of 0.063 µg/L and 0.07 µg/L, but no results are presented for the CV = 10% cutoffs. How can conclusions be drawn about the relative merits of 2 alternate cutoffs when only 1 is considered?

To use the 99th percentile irrespective of assay imprecision is an attractive suggestion because this would be simpler and less arbitrary than a 10% CV cutoff (which corresponds to a different cTnI concentration on every platform). Although we applaud the authors for attempting to evaluate this problem in a systematic manner, this conclusion, however attractive, is not supported by the evidence provided; and the issue is certainly not closed, as was asserted with reference to the results of this technical brief (6). Mathematical models are only as good as the assumptions that underpin them. If the assumptions are incorrect or not explicitly stated, then the results will be misleading and difficult to asses.


References

  1. Apple FS, Parvin CA, Buechler KF, Christenson RH, Wu AH, Jaffe AS. Validation of the 99th percentile cutoff independent of assay imprecision (CV) for cardiac troponin monitoring for ruling out myocardial infarction. Clin Chem 2005;51:2198-2200.[Free Full Text]
  2. Linnet K, Boyd JC. Selection and analytical evaluation of methods - with statistical techniques. Burtis CA Ashwood ER Bruns DE eds. Tietz textbook of clinical chemistry and molecular diagnostics 4th ed. 2006:353-407 Elsevier Saunders St. Louis. .
  3. Raggatt PR. Duplicates or singletons? An analysis of the need for replication in immunoassay and a computer program to calculate the distribution of outliers, error rate, and the precision profile from assay duplicates. Ann Clin Biochem 1989;26:26-37.
  4. Panteghini M, Pagani F, Yeo KT, Apple FS, Christenson RH, Dati F, et al. Evaluation of imprecision for cardiac troponin assays at low-range concentrations. Clin Chem 2004;50:327-332.[Abstract/Free Full Text]
  5. Sadler WA, Smith MH, Legge HM. A method for direct estimation of imprecision profiles, with reference to immunoassay data. Clin Chem 1988;34:1058-1061.[Abstract/Free Full Text]
  6. Kupchak P, Wu AHB, Ghani F, Newby LK, Ohman EM, Christenson RH. Influence of imprecision on ROC curve analysis for cardiac markers. Clin Chem 2006;52:752-753.[Abstract/Free Full Text]
  7. Eric W. Weisstein. "Convolution". From MathWorld–A Wolfram Web Resource. http://mathworld.wolfram.com/Convolution.html (Accessed February 2006).




This Article
Right arrow Extract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Holmes, D. T.
Right arrow Articles by Buhr, K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Holmes, D. T.
Right arrow Articles by Buhr, K.
Related Collections
Right arrow Evidence Based Laboratory Medicine and Test Utilization


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS