Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Validation of Tikhonov adaptively regularized gamma variate fitting with 24-h plasma clearance in cirrhotic patients with ascites</title>
<author><name sortKey="Wesolowski, Carl Adam" sort="Wesolowski, Carl Adam" uniqKey="Wesolowski C" first="Carl Adam" last="Wesolowski">Carl Adam Wesolowski</name>
<affiliation><nlm:aff id="Aff1">Radiology, Memorial University of Newfoundland, St. John’s, NL Canada</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff6">Nuclear Medicine, The General Hospital, HSC, 300 Prince Philip Drive, St. John’s, Newfoundland Canada A1B 3V6</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ling, Lin" sort="Ling, Lin" uniqKey="Ling L" first="Lin" last="Ling">Lin Ling</name>
<affiliation><nlm:aff id="Aff1">Radiology, Memorial University of Newfoundland, St. John’s, NL Canada</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Xirouchakis, Elias" sort="Xirouchakis, Elias" uniqKey="Xirouchakis E" first="Elias" last="Xirouchakis">Elias Xirouchakis</name>
<affiliation><nlm:aff id="Aff2">The Royal Free Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Burniston, Maria T" sort="Burniston, Maria T" uniqKey="Burniston M" first="Maria T." last="Burniston">Maria T. Burniston</name>
<affiliation><nlm:aff id="Aff3">Department of Medical Physics, Royal Free Hampstead NHS Trust, London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Puetter, Richard C" sort="Puetter, Richard C" uniqKey="Puetter R" first="Richard C." last="Puetter">Richard C. Puetter</name>
<affiliation><nlm:aff id="Aff4">Center for Astrophysics & Space Sciences, University of California, San Diego, La Jolla, CA United States</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Babyn, Paul S" sort="Babyn, Paul S" uniqKey="Babyn P" first="Paul S." last="Babyn">Paul S. Babyn</name>
<affiliation><nlm:aff id="Aff5">Radiology, University of Saskatchewan, Saskatoon, SK Canada</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Giamalis, Ioannis G" sort="Giamalis, Ioannis G" uniqKey="Giamalis I" first="Ioannis G." last="Giamalis">Ioannis G. Giamalis</name>
<affiliation><nlm:aff id="Aff2">The Royal Free Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Burroughs, Andrew K" sort="Burroughs, Andrew K" uniqKey="Burroughs A" first="Andrew K." last="Burroughs">Andrew K. Burroughs</name>
<affiliation><nlm:aff id="Aff2">The Royal Free Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">21881927</idno>
<idno type="pmc">3218286</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3218286</idno>
<idno type="RBID">PMC:3218286</idno>
<idno type="doi">10.1007/s00259-011-1887-9</idno>
<date when="2011">2011</date>
<idno type="wicri:Area/Pmc/Corpus">000D26</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000D26</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Validation of Tikhonov adaptively regularized gamma variate fitting with 24-h plasma clearance in cirrhotic patients with ascites</title>
<author><name sortKey="Wesolowski, Carl Adam" sort="Wesolowski, Carl Adam" uniqKey="Wesolowski C" first="Carl Adam" last="Wesolowski">Carl Adam Wesolowski</name>
<affiliation><nlm:aff id="Aff1">Radiology, Memorial University of Newfoundland, St. John’s, NL Canada</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff6">Nuclear Medicine, The General Hospital, HSC, 300 Prince Philip Drive, St. John’s, Newfoundland Canada A1B 3V6</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ling, Lin" sort="Ling, Lin" uniqKey="Ling L" first="Lin" last="Ling">Lin Ling</name>
<affiliation><nlm:aff id="Aff1">Radiology, Memorial University of Newfoundland, St. John’s, NL Canada</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Xirouchakis, Elias" sort="Xirouchakis, Elias" uniqKey="Xirouchakis E" first="Elias" last="Xirouchakis">Elias Xirouchakis</name>
<affiliation><nlm:aff id="Aff2">The Royal Free Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Burniston, Maria T" sort="Burniston, Maria T" uniqKey="Burniston M" first="Maria T." last="Burniston">Maria T. Burniston</name>
<affiliation><nlm:aff id="Aff3">Department of Medical Physics, Royal Free Hampstead NHS Trust, London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Puetter, Richard C" sort="Puetter, Richard C" uniqKey="Puetter R" first="Richard C." last="Puetter">Richard C. Puetter</name>
<affiliation><nlm:aff id="Aff4">Center for Astrophysics & Space Sciences, University of California, San Diego, La Jolla, CA United States</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Babyn, Paul S" sort="Babyn, Paul S" uniqKey="Babyn P" first="Paul S." last="Babyn">Paul S. Babyn</name>
<affiliation><nlm:aff id="Aff5">Radiology, University of Saskatchewan, Saskatoon, SK Canada</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Giamalis, Ioannis G" sort="Giamalis, Ioannis G" uniqKey="Giamalis I" first="Ioannis G." last="Giamalis">Ioannis G. Giamalis</name>
<affiliation><nlm:aff id="Aff2">The Royal Free Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Burroughs, Andrew K" sort="Burroughs, Andrew K" uniqKey="Burroughs A" first="Andrew K." last="Burroughs">Andrew K. Burroughs</name>
<affiliation><nlm:aff id="Aff2">The Royal Free Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">European Journal of Nuclear Medicine and Molecular Imaging</title>
<idno type="ISSN">1619-7070</idno>
<idno type="eISSN">1619-7089</idno>
<imprint><date when="2011">2011</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Purpose</title>
<p>The aim was to compare late-time extrapolation of plasma clearance (<italic>CL</italic>
) from Tikhonov adaptively regularized gamma variate fitting (Tk-GV) and from mono-exponential (E1) fitting.</p>
</sec>
<sec><title>Methods</title>
<p>Ten <sup>51</sup>
Cr-ethylenediaminetetraacetic acid bolus IV studies in adults—8 with ascites—assessed for liver transplantation, with 12–16 plasma samples drawn from 5-min to 24-h, were fit with Tk-GV and E1 models and <italic>CL</italic>
results were compared using Passing-Bablok fitting.</p>
</sec>
<sec><title>Results</title>
<p>The 24-h <italic>CL</italic>
(Tk-GV) values ranged from 11.4 to 79.7 ml/min. Linear regression of 4- versus 24-h <italic>CL</italic>
(Tk-GV) yielded no significant departure from a slope of 1, whereas the 4- versus 24-h <italic>CL</italic>
(E1) slope, 1.56, was significantly increased. For <italic>CL</italic>
(Tk-GV-24-h) versus <italic>CL</italic>
(E1-24-h), there was a biased slope and intercept (0.85, 5.97 ml/min). Moreover, the quality of fitting of 24-h data was significantly better for Tk-GV than for E1, as follows. For 10 logarithm of concentration curves, higher <italic>r</italic>
values were obtained for each Tk-GV fit (median 0.998) than for its corresponding E1 fit (median 0.965), with <italic>p</italic>
< 0.0001 (paired <italic>t</italic>
-test of <italic>z</italic>
-statistics from Fisher <italic>r-z</italic>
transformations). The E1 fit quality degraded with increasing <italic>V/W</italic>
[volume of distribution (l) per kg body weight, <italic>p</italic>
= 0.003]. However, Tk-GV fit quality versus <italic>V/W</italic>
was uncorrelated (<italic>p</italic>
= 0.8).</p>
</sec>
<sec><title>Conclusion</title>
<p><italic>CL</italic>
(E1) values were dependent on sample time and the quality of fit was poor and degraded with increasing ascites, consistent with current opinion that <italic>CL</italic>
(E1) is contraindicated in ascitic patients. <italic>CL</italic>
(Tk-GV) was relatively more accurate and the good quality of fit was unaffected by ascites. <italic>CL</italic>
(Tk-GV) was the preferred method for the accurate calculation of <italic>CL</italic>
and was useful despite liver failure and ascites.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Calvert, Ah" uniqKey="Calvert A">AH Calvert</name>
</author>
<author><name sortKey="Egorin, Mj" uniqKey="Egorin M">MJ Egorin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fleming, Js" uniqKey="Fleming J">JS Fleming</name>
</author>
<author><name sortKey="Zivanovic, Ma" uniqKey="Zivanovic M">MA Zivanovic</name>
</author>
<author><name sortKey="Blake, Gm" uniqKey="Blake G">GM Blake</name>
</author>
<author><name sortKey="Burniston, M" uniqKey="Burniston M">M Burniston</name>
</author>
<author><name sortKey="Cosgriff, Ps" uniqKey="Cosgriff P">PS Cosgriff</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Florijn, Kw" uniqKey="Florijn K">KW Florijn</name>
</author>
<author><name sortKey="Barendregt, Jnm" uniqKey="Barendregt J">JNM Barendregt</name>
</author>
<author><name sortKey="Lentjes, Egw" uniqKey="Lentjes E">EGW Lentjes</name>
</author>
<author><name sortKey="Dam, W" uniqKey="Dam W">W Dam</name>
</author>
<author><name sortKey="Prodjosudjadi, W" uniqKey="Prodjosudjadi W">W Prodjosudjadi</name>
</author>
<author><name sortKey="Saase, Jl" uniqKey="Saase J">JL Saase</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Henriksen, Jh" uniqKey="Henriksen J">JH Henriksen</name>
</author>
<author><name sortKey="Br Chner Mortensen, J" uniqKey="Br Chner Mortensen J">J Brøchner-Mortensen</name>
</author>
<author><name sortKey="Malchow M Ller, A" uniqKey="Malchow M Ller A">A Malchow-Møller</name>
</author>
<author><name sortKey="Schlichting, P" uniqKey="Schlichting P">P Schlichting</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Moore, Ae" uniqKey="Moore A">AE Moore</name>
</author>
<author><name sortKey="Park Holohan, Sj" uniqKey="Park Holohan S">SJ Park-Holohan</name>
</author>
<author><name sortKey="Blake, Gm" uniqKey="Blake G">GM Blake</name>
</author>
<author><name sortKey="Fogelman, I" uniqKey="Fogelman I">I Fogelman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nielsen, Ss" uniqKey="Nielsen S">SS Nielsen</name>
</author>
<author><name sortKey="Havsteen, H" uniqKey="Havsteen H">H Havsteen</name>
</author>
<author><name sortKey="Petersen, Lk" uniqKey="Petersen L">LK Petersen</name>
</author>
<author><name sortKey="Nielsen, Le" uniqKey="Nielsen L">LE Nielsen</name>
</author>
<author><name sortKey="Rehling, M" uniqKey="Rehling M">M Rehling</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Passing, H" uniqKey="Passing H">H Passing</name>
</author>
<author><name sortKey="Bablok, W" uniqKey="Bablok W">W Bablok</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Perkinson, As" uniqKey="Perkinson A">AS Perkinson</name>
</author>
<author><name sortKey="Evans, Cj" uniqKey="Evans C">CJ Evans</name>
</author>
<author><name sortKey="Burniston, Mt" uniqKey="Burniston M">MT Burniston</name>
</author>
<author><name sortKey="Smye, Sw" uniqKey="Smye S">SW Smye</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rehling, M" uniqKey="Rehling M">M Rehling</name>
</author>
<author><name sortKey="Stadeager, C" uniqKey="Stadeager C">C Stadeager</name>
</author>
<author><name sortKey="Henriksen, Jh" uniqKey="Henriksen J">JH Henriksen</name>
</author>
<author><name sortKey="Siemsen, O" uniqKey="Siemsen O">O Siemsen</name>
</author>
<author><name sortKey="Krintel, Jj" uniqKey="Krintel J">JJ Krintel</name>
</author>
<author><name sortKey="Malchow M Ller, A" uniqKey="Malchow M Ller A">A Malchow-Møller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Russell, Cd" uniqKey="Russell C">CD Russell</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Russell, Cd" uniqKey="Russell C">CD Russell</name>
</author>
<author><name sortKey="Taylor, At" uniqKey="Taylor A">AT Taylor</name>
</author>
<author><name sortKey="Dubovsky, Ev" uniqKey="Dubovsky E">EV Dubovsky</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wesolowski, Ca" uniqKey="Wesolowski C">CA Wesolowski</name>
</author>
<author><name sortKey="Hogendoorn, P" uniqKey="Hogendoorn P">P Hogendoorn</name>
</author>
<author><name sortKey="Vandierendonck, R" uniqKey="Vandierendonck R">R Vandierendonck</name>
</author>
<author><name sortKey="Driedger, Aa" uniqKey="Driedger A">AA Driedger</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wesolowski, Ca" uniqKey="Wesolowski C">CA Wesolowski</name>
</author>
<author><name sortKey="Puetter, Rc" uniqKey="Puetter R">RC Puetter</name>
</author>
<author><name sortKey="Ling, L" uniqKey="Ling L">L Ling</name>
</author>
<author><name sortKey="Babyn, Ps" uniqKey="Babyn P">PS Babyn</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xirouchakis, E" uniqKey="Xirouchakis E">E Xirouchakis</name>
</author>
<author><name sortKey="Marelli, L" uniqKey="Marelli L">L Marelli</name>
</author>
<author><name sortKey="Cholongitas, E" uniqKey="Cholongitas E">E Cholongitas</name>
</author>
<author><name sortKey="Manousou, P" uniqKey="Manousou P">P Manousou</name>
</author>
<author><name sortKey="Calvaruso, V" uniqKey="Calvaruso V">V Calvaruso</name>
</author>
<author><name sortKey="Pleguezuelo, M" uniqKey="Pleguezuelo M">M Pleguezuelo</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Eur J Nucl Med Mol Imaging</journal-id>
<journal-title-group><journal-title>European Journal of Nuclear Medicine and Molecular Imaging</journal-title>
</journal-title-group>
<issn pub-type="ppub">1619-7070</issn>
<issn pub-type="epub">1619-7089</issn>
<publisher><publisher-name>Springer-Verlag</publisher-name>
<publisher-loc>Berlin/Heidelberg</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">21881927</article-id>
<article-id pub-id-type="pmc">3218286</article-id>
<article-id pub-id-type="publisher-id">1887</article-id>
<article-id pub-id-type="doi">10.1007/s00259-011-1887-9</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Original Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Validation of Tikhonov adaptively regularized gamma variate fitting with 24-h plasma clearance in cirrhotic patients with ascites</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes"><name><surname>Wesolowski</surname>
<given-names>Carl Adam</given-names>
</name>
<address><phone>+1-709-7778999</phone>
<fax>+1-709-7778267</fax>
<email>carl.wesolowski@gmail.com</email>
</address>
<xref ref-type="aff" rid="Aff1">1</xref>
<xref ref-type="aff" rid="Aff6">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ling</surname>
<given-names>Lin</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Xirouchakis</surname>
<given-names>Elias</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Burniston</surname>
<given-names>Maria T.</given-names>
</name>
<xref ref-type="aff" rid="Aff3">3</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Puetter</surname>
<given-names>Richard C.</given-names>
</name>
<xref ref-type="aff" rid="Aff4">4</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Babyn</surname>
<given-names>Paul S.</given-names>
</name>
<xref ref-type="aff" rid="Aff5">5</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Giamalis</surname>
<given-names>Ioannis G.</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Burroughs</surname>
<given-names>Andrew K.</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<aff id="Aff1"><label>1</label>
Radiology, Memorial University of Newfoundland, St. John’s, NL Canada</aff>
<aff id="Aff2"><label>2</label>
The Royal Free Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK</aff>
<aff id="Aff3"><label>3</label>
Department of Medical Physics, Royal Free Hampstead NHS Trust, London, UK</aff>
<aff id="Aff4"><label>4</label>
Center for Astrophysics & Space Sciences, University of California, San Diego, La Jolla, CA United States</aff>
<aff id="Aff5"><label>5</label>
Radiology, University of Saskatchewan, Saskatoon, SK Canada</aff>
<aff id="Aff6"><label>6</label>
Nuclear Medicine, The General Hospital, HSC, 300 Prince Philip Drive, St. John’s, Newfoundland Canada A1B 3V6</aff>
</contrib-group>
<pub-date pub-type="epub"><day>1</day>
<month>9</month>
<year>2011</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>1</day>
<month>9</month>
<year>2011</year>
</pub-date>
<pub-date pub-type="ppub"><month>12</month>
<year>2011</year>
</pub-date>
<volume>38</volume>
<issue>12</issue>
<fpage>2247</fpage>
<lpage>2256</lpage>
<history><date date-type="received"><day>10</day>
<month>4</month>
<year>2011</year>
</date>
<date date-type="accepted"><day>12</day>
<month>7</month>
<year>2011</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2011</copyright-statement>
</permissions>
<abstract id="Abs1"><sec><title>Purpose</title>
<p>The aim was to compare late-time extrapolation of plasma clearance (<italic>CL</italic>
) from Tikhonov adaptively regularized gamma variate fitting (Tk-GV) and from mono-exponential (E1) fitting.</p>
</sec>
<sec><title>Methods</title>
<p>Ten <sup>51</sup>
Cr-ethylenediaminetetraacetic acid bolus IV studies in adults—8 with ascites—assessed for liver transplantation, with 12–16 plasma samples drawn from 5-min to 24-h, were fit with Tk-GV and E1 models and <italic>CL</italic>
results were compared using Passing-Bablok fitting.</p>
</sec>
<sec><title>Results</title>
<p>The 24-h <italic>CL</italic>
(Tk-GV) values ranged from 11.4 to 79.7 ml/min. Linear regression of 4- versus 24-h <italic>CL</italic>
(Tk-GV) yielded no significant departure from a slope of 1, whereas the 4- versus 24-h <italic>CL</italic>
(E1) slope, 1.56, was significantly increased. For <italic>CL</italic>
(Tk-GV-24-h) versus <italic>CL</italic>
(E1-24-h), there was a biased slope and intercept (0.85, 5.97 ml/min). Moreover, the quality of fitting of 24-h data was significantly better for Tk-GV than for E1, as follows. For 10 logarithm of concentration curves, higher <italic>r</italic>
values were obtained for each Tk-GV fit (median 0.998) than for its corresponding E1 fit (median 0.965), with <italic>p</italic>
< 0.0001 (paired <italic>t</italic>
-test of <italic>z</italic>
-statistics from Fisher <italic>r-z</italic>
transformations). The E1 fit quality degraded with increasing <italic>V/W</italic>
[volume of distribution (l) per kg body weight, <italic>p</italic>
= 0.003]. However, Tk-GV fit quality versus <italic>V/W</italic>
was uncorrelated (<italic>p</italic>
= 0.8).</p>
</sec>
<sec><title>Conclusion</title>
<p><italic>CL</italic>
(E1) values were dependent on sample time and the quality of fit was poor and degraded with increasing ascites, consistent with current opinion that <italic>CL</italic>
(E1) is contraindicated in ascitic patients. <italic>CL</italic>
(Tk-GV) was relatively more accurate and the good quality of fit was unaffected by ascites. <italic>CL</italic>
(Tk-GV) was the preferred method for the accurate calculation of <italic>CL</italic>
and was useful despite liver failure and ascites.</p>
</sec>
</abstract>
<kwd-group><title>Keywords</title>
<kwd>Renal function</kwd>
<kwd>Tikhonov regularization</kwd>
<kwd>Minimization of errors</kwd>
<kwd>Gamma variate</kwd>
<kwd>Ascites</kwd>
</kwd-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© Springer-Verlag 2011</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1"><title>Introduction</title>
<p>This paper presents a validation of a curve fit method for measuring plasma clearance (<italic>CL</italic>
, ml/min). The <italic>CL</italic>
method validated here uses <bold>T</bold>
i<bold>k</bold>
honov adaptively regularized <bold>g</bold>
amma <bold>v</bold>
ariate (Tk-GV) fitting. The word “adaptive” is used to imply minimization of relative error of <italic>CL</italic>
as opposed to an ordinary curve fit minimization (<xref rid="Sec6" ref-type="sec">Appendix</xref>
and [<xref ref-type="bibr" rid="CR13">13</xref>
]). In general, radiometric <italic>CL</italic>
is indicated:<list list-type="bullet"><list-item><p>To detect nephrotoxicity of chemotherapy drugs, especially in children</p>
</list-item>
<list-item><p>For many chemotherapeutic agents’ dose calculations</p>
</list-item>
<list-item><p>For detection of renal failure in patients in whom (a) serum creatinine results might be misleading, (b) missing a decline in renal function might be disastrous e.g. single kidney, renovascular disease or renal transplant and (c) a 24-h clearance measurement is difficult, e.g., in the elderly or those with learning difficulties</p>
</list-item>
<list-item><p>For assessment of potential live donors for kidney transplantation</p>
</list-item>
<list-item><p>For the evaluation and follow-up of chronic renal disease</p>
</list-item>
<list-item><p>For the evaluation of single kidney function in conjunction with relative renal function measurements from static or dynamic radionuclide imaging [<xref ref-type="bibr" rid="CR2">2</xref>
]</p>
</list-item>
</list>
</p>
<p>However, utilization of radiometric <italic>CL</italic>
methods is spotty with only occasional centres performing them, with needed but inconsistently applied <italic>CL</italic>
-value correction factors, without quality control of the individual <italic>CL</italic>
-values, and in certain cases, current methods are relatively contraindicated. For example, liver failure, ascites, oedema and low <italic>CL</italic>
are identified in the 2004 British Nuclear Medicine Society Guidelines as having problematically inaccurate <italic>CL</italic>
-values [<xref ref-type="bibr" rid="CR2">2</xref>
]. On the other hand, Xirouchakis et al. [<xref ref-type="bibr" rid="CR14">14</xref>
] outline the clinical importance of <italic>CL</italic>
measurements in liver failure: “Renal dysfunction is a well-established predictor of mortality in both acute liver failure and cirrhosis, particularly following the development of complications, such as sepsis. Inulin clearance and other direct methods using injected exogenous radiolabelled substances [<sup>51</sup>
Cr-EDTA (ethylenediamine tetra-acetic acid), <sup>125</sup>
I-Iothalamate, <sup>99m</sup>
Tc-DTPA (diethylenetriamine penta-acetic acid)] are the most accurate to assess renal function.” Thus, the need for <italic>CL</italic>
measurements in liver failure is clear, but current radiometric methods are inadequate for this purpose [<xref ref-type="bibr" rid="CR2">2</xref>
, <xref ref-type="bibr" rid="CR4">4</xref>
, <xref ref-type="bibr" rid="CR6">6</xref>
, <xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR9">9</xref>
].</p>
<p>In this study, Tk-GV <italic>CL</italic>
values will be obtained for ten liver failure patients most of whom (8/10) have ascites. Investigated here is whether the Tk-GV models overcome the need for <italic>CL</italic>
-value correction factors and maintain good quality even in the presence of ascites. In prior published work, the Tk-GV method has been applied to 4-h data using <sup>99m</sup>
Tc-DTPA and <sup>169</sup>
Yb-DTPA [<xref ref-type="bibr" rid="CR13">13</xref>
]. Tk-GV results provide quality control for each individual <italic>CL</italic>
-value. The current paper is a first application of Tk-GV to <sup>51</sup>
Cr-EDTA data, and to sampling data extended to 24-h. There is no a priori guarantee that Tk-GV will either usefully fit an untried pharmaceutical or 24-h data so testing is needed. The Tk-GV results are contrasted with those from mono-exponential (E1) fitting. E1 models fit to this new data are, from prior work, fairly certain to be problematic [<xref ref-type="bibr" rid="CR13">13</xref>
], especially in ascites [<xref ref-type="bibr" rid="CR2">2</xref>
, <xref ref-type="bibr" rid="CR4">4</xref>
, <xref ref-type="bibr" rid="CR6">6</xref>
, <xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR9">9</xref>
]. Thus, how poorly E1 performs for our 24-h data is also of interest.</p>
<p>Intravenous (IV) bolus-injected inert markers clear into the interstitium as well as in the urine. Consequently, <italic>CL</italic>
is slightly faster than glomerular filtration rate (<italic>GFR</italic>
) [<xref ref-type="bibr" rid="CR3">3</xref>
, <xref ref-type="bibr" rid="CR5">5</xref>
, <xref ref-type="bibr" rid="CR13">13</xref>
], e.g. <italic>CL</italic>
≈ 1.076 <italic>GFR</italic>
[<xref ref-type="bibr" rid="CR5">5</xref>
]. After IV bolus injection, <italic>CL</italic>
estimates from multiple plasma samples are calculated from the area under the plasma marker disappearance curve (<italic>AUC</italic>
) from time equals zero to infinity, i.e. <italic>CL</italic>
= dose/<italic>AUC</italic>
. <italic>AUC</italic>
has units of concentration min. and marker (e.g. drug) systemic exposure time is plasma clearance (<italic>CL</italic>
) related [<xref ref-type="bibr" rid="CR1">1</xref>
, <xref ref-type="bibr" rid="CR13">13</xref>
]. Even the most minimalistic technique for estimating <italic>AUC</italic>
, e.g. numerical integration [<xref ref-type="bibr" rid="CR5">5</xref>
], requires the use of a curve fit to extrapolate concentrations from the time of last sample to infinite time. In this paper, the E1 and Tk-GV curve fit methods for finding <italic>AUC</italic>
are tested to determine the accuracy of <italic>CL</italic>
from curve fitting of 4- versus 24-h data for ten <sup>51</sup>
Cr-EDTA concentration curves. A key issue explored is how well clearances determined only from 4-h data predict the results determined from 24-h data. When the 4-h <italic>CL</italic>
curve does a good job of predicting the 24-h concentrations, the model extrapolates the late-time behaviour well, and the <italic>AUC</italic>
values give accurate clearances. Conversely, if the 4-h <italic>CL</italic>
curve fits do not predict the 24-h concentrations, then the ability of the fitting procedure to predict accurately the late-time behaviour is in question and the clearance values predicted from the method are suspect.</p>
<p>On a lnE1 versus <italic>t</italic>
plot, lnE1 models are, of course, straight lines. However, let <italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) be the observed concentrations. Then, on a ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus <italic>t</italic>
plot, the data are curved; thus E1 models do not fit data well. This same trend, i.e. curvature, is also seen in a plot of ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus lnE1, and the Pearson <italic>r</italic>
-value is less than 1. Indeed, from ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus lnE1, the best R<sup>2</sup>
-value that can be achieved by a linear fit (i.e. an E1 model of any kind, fit to any subset of the data) is the Pearson <italic>r</italic>
-value squared, i.e. <italic>r</italic>
<sup>2</sup>
, as this is the limitation of the data’s linearity. To achieve a better fit and to gain better predictability, one needs to use curves that follow the curvature of the ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus <italic>t</italic>
plot. As we shall see, such a curve is the Tk-GV model. If such a model successfully follows the data’s curvature, the observed log concentrations, ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
), will be much better correlated with the model predictions. In other words, in a plot of ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus ln<italic>C</italic>
<sub>Tk-GV</sub>
(<italic>t</italic>
), the relationship should be almost perfectly linear, unlike the ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus ln<italic>C</italic>
<sub>E1</sub>
(<italic>t</italic>
) plots that show curvature and thus imperfect correlation.</p>
</sec>
<sec id="Sec2" sec-type="materials|methods"><title>Materials and methods</title>
<p>The patients are part of a research project comparing methods of assessing renal function in cirrhosis approved by the Royal Free Hospital Research Ethics Committee. Thirteen patients with cirrhosis under assessment for liver transplantation had <sup>51</sup>
Cr-EDTA bolus IV studies. Of these, three patients did not have plasma samples drawn at various times, e.g. at 12- or 24-h, and were discarded. Of the remaining ten patients, four had hepatitis C-related cirrhosis (one also with hepatocellular carcinoma), three had alcoholic cirrhosis, two had primary biliary cirrhosis, and one had non-alcoholic, fatty liver disease-associated cirrhosis. Eight were men and two women: eight were white, one black and one South Asian Indian. The median weight was 79.5 kg (range 51.3–105.8), median serum bilirubin was 66 mmol/l (range 13–390), median serum albumin was 35 g/l (range 26–40) and median INR was 1.7 (range 1.5–2.4). Hepatic encephalopathy was present in eight and ascites in eight. The patients had 12–16 plasma samples drawn from 5-min to 24-h. Their concentration curves were fit with Tk-GV and E1 models and <italic>CL</italic>
results were compared using Passing-Bablok fitting and paired (correlated) samples <italic>t</italic>
-tests. A minimum of two samples determines an E1 solution, and fits to data from 5-min to 4-, 12- and 24-h are examined here. Also examined is a recommended clinical method of obtaining E1 functions from plasma samples at 2- and at 4-h [<xref ref-type="bibr" rid="CR2">2</xref>
]. Exponential models have uniform concentrations (i.e. are well mixed, aka, instant mixing) in one or more partial volumes. For Tk-GV fitting, a minimum of four samples are needed. The GV concentration, <italic>C</italic>
(<italic>t</italic>
), model used for Tk-GV fitting is <inline-formula id="IEq9"><alternatives><tex-math id="M1">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\text{GV}} = C(t) = K{t^{{\alpha - 1}}}{e^{{ - \beta t}}} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq9.gif"></inline-graphic>
</alternatives>
</inline-formula>
, where <italic>C</italic>
(<italic>t</italic>
) is the marker concentration, <italic>α</italic>
is the dimensionless volume of distribution scaling factor with <italic>V = αCL/β</italic>
, <italic>β</italic>
is the renal elimination rate constant, <italic>K</italic>
is a constant, and the mean residence time <italic>MRT = α/β</italic>
. The Tk-GV method uses least relative error of <italic>CL</italic>
as a minimization target (see the <xref rid="Sec6" ref-type="sec">Appendix</xref>
) and is more robust than GV fitting with least-squares regression. The fit results for well-posed GV <italic>CL</italic>
-values (0 < <italic>α</italic>
≤ 1, <italic>β</italic>
> 0, and fits from 5 or 10 min to at least 3-h of data) are extensively described elsewhere [<xref ref-type="bibr" rid="CR13">13</xref>
]. In specific, the 0 < <italic>α</italic>
≤ 1 constraint is <italic>only</italic>
observed with Tk-GV fitting and does not occur with least-squares GV fitting. A Tk-GV model of an inert marker is a poorly mixed system with a concentration gradient such that the physical volume of the system is factor of <italic>α</italic>
times smaller than the <italic>virtual</italic>
volume (where the latter is <italic>CL</italic>
/<italic>β</italic>
). Moreover, the dynamic equilibrium concentration gradient is established by renal function and disappears, i.e. <italic>α</italic>
→ 1, only when renal function approaches zero [<xref ref-type="bibr" rid="CR13">13</xref>
].</p>
<p>Passing-Bablok fitting is used here for comparison of methods to provide unbiased linear regression slope and intercept [<xref ref-type="bibr" rid="CR7">7</xref>
]. The Passing-Bablok method includes the cusum test for linearity [<xref ref-type="bibr" rid="CR7">7</xref>
], which is useful for comparison of methods because it detects nonlinear relationships. Cusum <italic>p</italic>
> 0.1 does not suggest nonlinearity, 0.05 < <italic>p</italic>
< 0.1 suggests nonlinearity, and <italic>p</italic>
< 0.05 is highly suggestive of nonlinearity. Spearman rank correlation is used once for extracting a correlation in a slightly nonlinear relationship. Paired (i.e. two, correlated) samples <italic>t</italic>
-testing is used here for comparison of methods using all the data. Two-sample <italic>t</italic>
-testing is also applied to (paired) <italic>z</italic>
-statistics from Fisher <italic>r-z</italic>
transformation for detecting improved correlation. This transformation changes the otherwise non-normally distributed correlation coefficients (<italic>r</italic>
) to more normally distributed <italic>z</italic>
-statistics, using the inverse hyperbolic tangent function, <italic>z</italic>
= atanh (<italic>r</italic>
). <italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus ln(<italic>t</italic>
) plots are used to explain modelling effects. However, for calculating model fit quality correlations, we use the logarithm of both the observed concentrations and the predicted concentrations from the E1 and Tk-GV models. To test how linear the ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) versus ln<italic>C</italic>
<sub>Model</sub>
(<italic>t</italic>
) relationship is, one simply calculates its Pearson <italic>r</italic>
value. (The correlation coefficient of determination, R<sup>2</sup>
, is then the Pearson <italic>r</italic>
-value squared, i.e. <italic>r</italic>
<sup>2</sup>
.) The test for ascites looks at the trend of fit quality, <italic>r</italic>
, versus <italic>V/W</italic>
, where <italic>V</italic>
is the volume of distribution (ml) of the Tk-GV model, and <italic>W</italic>
is the patient mass in kilograms. A model <italic>r</italic>
-value that degrades with increasing ascites (as indexed to <italic>V/W</italic>
) suggests a poorly behaved model. If the model <italic>r</italic>
-value is unchanged by increasing ascites, the model is relatively bulletproof.</p>
</sec>
<sec id="Sec3"><title>Results</title>
<p>Figure <xref rid="Fig1" ref-type="fig">1</xref>
shows plots of 24-h concentrations of <sup>51</sup>
Cr-EDTA in counts per ml min, <italic>C</italic>
(<italic>t</italic>
), versus the logarithm of time, <italic>t</italic>
, in min. Figure <xref rid="Fig1" ref-type="fig">1</xref>
plots show the highest and lowest <italic>CL</italic>
(Tk-GV) cases (79.7 and 11.4 ml/min) in this series. These plots suggest an improvement in fit quality for 24-h data fits using Tk-GV functions compared to the E1 functions. <italic>C</italic>
(<italic>t</italic>
) versus ln(<italic>t</italic>
) plots roughly linearize the early-time <italic>C</italic>
(<italic>t</italic>
). Figure <xref rid="Fig1" ref-type="fig">1</xref>
demonstrates that the exponential fit assumption of a constant concentration at time zero causes underestimation of the early concentrations of actual data. Let us define late-time <italic>C</italic>
(<italic>t</italic>
) as being after the concentration curve becomes nonlinear in ln(<italic>t</italic>
), i.e. asymptotic to the <italic>C</italic>
(<italic>t</italic>
) = 0 axis. Figure <xref rid="Fig1" ref-type="fig">1</xref>
a demonstrates that when an exponential fit includes the late-time concentrations, late-time concentrations are accurately fit. Figure <xref rid="Fig1" ref-type="fig">1</xref>
b illustrates that when an exponential fit does not include late-time concentrations, the late-time concentrations are underestimated.<fig id="Fig1"><label>Fig. 1</label>
<caption><p><bold>a</bold>
Plot of 24-h concentrations of <sup>51</sup>
Cr-EDTA in counts per ml min, <italic>C</italic>
(<italic>t</italic>
), per logarithm of time, <italic>t</italic>
, in min for patient 3’s data. For patient 3’s 24-h data, R<sup>2</sup>
(Tk-GV) = 0.999 > R<sup>2</sup>
= 0.977 from a one exponential term model (E1). The <italic>CL</italic>
(Tk-GV) for patient 3 is the highest <italic>CL</italic>
in this series (79.7 ml/min). <bold>b</bold>
Plot of 24-h concentrations, <italic>C</italic>
(<italic>t</italic>
), per logarithm of time, <italic>t</italic>
, for patient 10’s data, having the lowest <italic>CL</italic>
(Tk-GV) in this series (11.4 ml/min). Right-hand areas contribute more to <italic>AUC</italic>
than those on the left do, for this plot type. For example in <bold>b</bold>
, the <italic>AUC</italic>
for the Tk-GV curve from 1,440 min to infinite time is 47.8% of the total area</p>
</caption>
<graphic xlink:href="259_2011_1887_Fig1_HTML" id="MO1"></graphic>
</fig>
</p>
<p>For several reasons, the correlations (<italic>r</italic>
) used for comparison of methods are calculated from the observed versus predicted logarithms of concentration values. These <italic>r</italic>
-values (Table <xref rid="Tab1" ref-type="table">1</xref>
) are only indirectly related to <italic>r</italic>
-values calculated without taking logarithms (e.g. the R<sup>2</sup>
in Fig. <xref rid="Fig1" ref-type="fig">1</xref>
). Comparing methods, the improvement in correlation to logarithm concentration of ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) and ln<italic>C</italic>
<sub><italic>Model</italic>
</sub>
(<italic>t</italic>
) was highly significant (<italic>p</italic>
< 0.0001) for Tk-GV functions compared to the E1 functions using paired samples <italic>t</italic>
-testing of the <italic>z</italic>
-values from <italic>z</italic>
= atanh(<italic>r</italic>
). Moreover, the <italic>r</italic>
-values are better for Tk-GV than for E1 fits in each of the ten cases (Table <xref rid="Tab1" ref-type="table">1</xref>
). With a median 24-h Tk-GV fit of <italic>r</italic>
= 0.998, there is little unexplained variance in the logarithm concentration data models (1−<italic>r</italic>
<sup>2</sup>
= 0.4%). This gives us confidence that extrapolation from 1,440 min to infinite time and from 5 min to zero time is accurate. Next, one uses Tk-GV fits to estimate the percentage of area under the curve (<italic>AUC</italic>
) of concentration versus time not directly under the samples. Not surprisingly, there is a lesser percentage of extrapolated <italic>AUC</italic>
for higher <italic>CL</italic>
(e.g. Table <xref rid="Tab1" ref-type="table">1</xref>
, case 3, 5.4%), and a higher percentage for lower <italic>CL</italic>
(case 10, 48.3%), and it takes longer than 24-h to clear <sup>51</sup>
Cr-EDTA in renal failure. However, the goodness of fit for Tk-GV allows us to have confidence that 24 h of <italic>AUC</italic>
data is long enough to calculate accurate <italic>CL</italic>
-values. An indication of <italic>CL</italic>
(Tk-GV) precision is given as standard deviation (SD) of <italic>CL</italic>
, which values range from 0.4 to 3.5 ml/min in Table <xref rid="Tab1" ref-type="table">1</xref>
. The Tk-GV method produces these estimates during processing. These estimates correspond to 95% confidence intervals for errors of 24-h <italic>CL</italic>
-values that range from ±0.8 to ±6.9 ml/min.<table-wrap id="Tab1"><label>Table 1</label>
<caption><p>From 24-h fitting to logarithms of concentrations for Tk-GV, <italic>CL</italic>
and <italic>V</italic>
(l) ± 1 standard deviation (SD) error, <italic>V</italic>
as percent of <italic>W</italic>
(weight-kg), and percentage of <italic>AUC</italic>
from extrapolation from left tails (0 to 5 min) and right tails (1,440 min to infinite time). Also shown are correlation coefficients, <italic>r</italic>
, for Tk-GV and E1 fitting</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Case</th>
<th><italic>CL</italic>
±1 SD (ml/min)</th>
<th><italic>V</italic>
±1 SD (l)</th>
<th><italic>V/W</italic>
(%)</th>
<th><italic>AUC</italic>
of tails (%)</th>
<th>Tk-GV <italic>r</italic>
</th>
<th>E1 <italic>r</italic>
all fits</th>
</tr>
</thead>
<tbody><tr><td>3</td>
<td>79.7 ± 0.4</td>
<td>29.7 ± 0.2</td>
<td>34.9</td>
<td>5.4</td>
<td>1.000</td>
<td>0.988</td>
</tr>
<tr><td>5</td>
<td>76.5 ± 1.6</td>
<td>38.9 ± 1.2</td>
<td>45.2</td>
<td>13.0</td>
<td>0.999</td>
<td>0.933</td>
</tr>
<tr><td>6</td>
<td>56.5 ± 0.9</td>
<td>21.3 ± 0.6</td>
<td>29.6</td>
<td>6.8</td>
<td>0.998</td>
<td>0.975</td>
</tr>
<tr><td>9</td>
<td>46.5 ± 0.8</td>
<td>18.6 ± 0.3</td>
<td>22.8</td>
<td>6.9</td>
<td>0.999</td>
<td>0.977</td>
</tr>
<tr><td>4</td>
<td>45.8 ± 0.8</td>
<td>21.9 ± 0.4</td>
<td>35.6</td>
<td>9.0</td>
<td>0.999</td>
<td>0.971</td>
</tr>
<tr><td>12</td>
<td>42.5 ± 2.1</td>
<td>26.9 ± 0.9</td>
<td>39.7</td>
<td>13.7</td>
<td>0.995</td>
<td>0.960</td>
</tr>
<tr><td>13</td>
<td>36.4 ± 3.1</td>
<td>41.9 ± 1.8</td>
<td>39.6</td>
<td>29.3</td>
<td>0.988</td>
<td>0.937</td>
</tr>
<tr><td>8</td>
<td>29.9 ± 3.5</td>
<td>60.8 ± 4.7</td>
<td>59.9</td>
<td>45.5</td>
<td>0.998</td>
<td>0.848</td>
</tr>
<tr><td>11</td>
<td>21.2 ± 0.8</td>
<td>17.7 ± 0.4</td>
<td>34.7</td>
<td>19.5</td>
<td>0.998</td>
<td>0.970</td>
</tr>
<tr><td>10</td>
<td>11.4 ± 0.9</td>
<td>23.4 ± 0.6</td>
<td>30.4</td>
<td>48.3</td>
<td>0.994</td>
<td>0.929</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
<p>As per the “<xref rid="Sec2" ref-type="sec">Materials and methods</xref>
” section, ln<italic>C</italic>
from E1 fit equations from any sample subsets have at best a quality of fit to ln<italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
) of the correlation coefficient evaluated at all the sample times (median as <italic>r</italic>
= 0.965, see Table <xref rid="Tab1" ref-type="table">1</xref>
). The correlation of <italic>r</italic>
-values from E1 24-h data fits with <italic>V/W</italic>
(Table <xref rid="Tab1" ref-type="table">1</xref>
) is significantly negative, <inline-formula id="IEq10"><alternatives><tex-math id="M2">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ r\left[ {r\left( {{\text{E}}1} \right),{ }V\,/\,W} \right] = - 0.823 $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq10.gif"></inline-graphic>
</alternatives>
</inline-formula>
, <italic>p</italic>
= 0.0034. However, the correlation of <italic>r</italic>
-values from 24-h Tk-GV fits to the data with the <italic>V/W</italic>
from the same Tk-GV model is uncorrelated, <inline-formula id="IEq11"><alternatives><tex-math id="M3">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ r\left[ {r\left( {{\text{Tk - GV}}} \right),{ }V\,/\,W} \right] = - 0.084 $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq11.gif"></inline-graphic>
</alternatives>
</inline-formula>
, <italic>p</italic>
= 0.8165.</p>
<p>Changing sample subsets can produce different <italic>CL</italic>
-values, especially for E1, as follows. With reference to Table <xref rid="Tab2" ref-type="table">2</xref>
and Fig. <xref rid="Fig2" ref-type="fig">2</xref>
a, the 4-h (vertical axis) versus 24-h (horizontal axis) <italic>CL</italic>
(Tk-GV) plot yielded no significant Passing-Bablok departure from a line slope of 1 or an intercept of 0. This means that 4-h Tk-GV clearances unbiasedly predict 24-h Tk-GV clearances. However, the R<sup>2</sup>
-value for this relationship is only 0.738, and the more powerful paired samples <italic>t</italic>
-test shows a (borderline) significant bias (<italic>p</italic>
= 0.04). This borderline significant mean difference is 8.6 ml/min greater <italic>CL</italic>
measured for the 4-h data than for the 24-h data. Much of the variance for predicting 24-h <italic>CL</italic>
from 4-h data is not related to the method of calculation but to the data themselves, i.e. from propagation of errors. For example, the furthest 4-h versus 24-h <italic>CL</italic>
(Tk-GV) outlier of Fig. <xref rid="Fig2" ref-type="fig">2</xref>
a, 64.4 versus 29.9 ml/min, is only a near outlier at 2.6 interquartile ranges (IQR) and would not usually require explanation. To decrease the 4- versus 24-h <italic>CL</italic>
(Tk-GV) difference from 11.5 to 3.9 ml/min (1 SD), one carries the fits out to 12-h (instead of 4-h) of data when comparing to the 24-h <italic>CL.</italic>
Note (Table <xref rid="Tab2" ref-type="table">2</xref>
) that for 12- versus 24-h, there is no significant departure from a line slope of 1 or an intercept of 0 and the R<sup>2</sup>
improves (0.968). However, the <italic>t</italic>
-test, 3.1 ml/min mean bias, although small, is still (borderline) significant (<italic>p</italic>
= 0.03).<table-wrap id="Tab2"><label>Table 2</label>
<caption><p>Tk-GV and E1 line fit tests for slope, <italic>m</italic>
=1, and intercept <italic>b</italic>
= 0 and test for significance of different mean values of <italic>CL</italic>
</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Passing-Bablok fit lines for</th>
<th>Parameter values</th>
<th></th>
<th colspan="2">95% CI</th>
<th>CI</th>
<th colspan="2">Paired <italic>t</italic>
test (Δ in ml/min)</th>
<th><italic>R</italic>
<sup>2</sup>
</th>
</tr>
<tr><th>(h) = 4, 12, 24, 2&4</th>
<th>(<italic>b</italic>
in ml/min)</th>
<th></th>
<th colspan="2">Low to high</th>
<th>Significance</th>
<th>Probability</th>
<th>Significance (Δ)</th>
<th></th>
</tr>
</thead>
<tbody><tr><td rowspan="2">Tk-GV(4) = <italic>m</italic>
•Tk-GV(24) + <italic>b</italic>
</td>
<td><italic>m</italic>
0.971</td>
<td></td>
<td colspan="2">0.568 to 1.50</td>
<td>NS</td>
<td>0.0416</td>
<td>B (8.6)</td>
<td>0.738</td>
</tr>
<tr><td><italic>b</italic>
</td>
<td>6.7</td>
<td colspan="2">−15.4 to 28.5</td>
<td>NS</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td rowspan="2">Tk-GV(12) = <italic>m</italic>
•Tk-GV(24) + <italic>b</italic>
</td>
<td><italic>m</italic>
</td>
<td>0.986</td>
<td colspan="2">0.795 to 1.15</td>
<td>NS</td>
<td>0.0337</td>
<td>B (3.1)</td>
<td>0.968</td>
</tr>
<tr><td><italic>b</italic>
</td>
<td>3.3</td>
<td colspan="2">−5.0 to 12.8</td>
<td>NS</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td rowspan="2">E1(4) = <italic>m</italic>
•E1(24) + <italic>b</italic>
</td>
<td><italic>m</italic>
</td>
<td>1.56</td>
<td colspan="2">1.16 to 9.77</td>
<td>S</td>
<td>0.0001</td>
<td>S (32.1)</td>
<td>0.738</td>
</tr>
<tr><td><italic>b</italic>
</td>
<td>5.2</td>
<td colspan="2">−347 to 17.4</td>
<td>NS</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td rowspan="2">E1(12) = <italic>m</italic>
•E1(24) + <italic>b</italic>
</td>
<td><italic>m</italic>
</td>
<td>1.11</td>
<td colspan="2">0.98 to 1.45</td>
<td>NS</td>
<td>0.0002</td>
<td>S (9.1)</td>
<td>0.951</td>
</tr>
<tr><td><italic>b</italic>
</td>
<td>3.0</td>
<td colspan="2">−10.7 to 7.63</td>
<td>NS</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td rowspan="2">E1(4) = <italic>m</italic>
•Tk-GV(4) + <italic>b</italic>
<sup>†</sup>
</td>
<td><italic>m</italic>
</td>
<td>1.31</td>
<td colspan="2">1.04 to 1.62</td>
<td>B</td>
<td><0.0001</td>
<td>S (23.9)</td>
<td>0.928<sup>a</sup>
</td>
</tr>
<tr><td><italic>b</italic>
</td>
<td>6.5</td>
<td colspan="2">−8.6 to 18.9</td>
<td>NS</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td rowspan="2">E1(24) = <italic>m</italic>
•Tk-GV(24) + <italic>b</italic>
</td>
<td><italic>m</italic>
</td>
<td>0.850</td>
<td colspan="2">0.417 to 0.940</td>
<td>S</td>
<td>0.8386</td>
<td>NS (0.4)</td>
<td>0.948</td>
</tr>
<tr><td><italic>b</italic>
</td>
<td>6.0</td>
<td colspan="2">1.0 to 26.6</td>
<td>B</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td rowspan="2">E1(2&4) = <italic>m</italic>
•Tk-GV(24) + <italic>b</italic>
</td>
<td><italic>m</italic>
</td>
<td>1.01</td>
<td colspan="2">0.626 to 1.59</td>
<td>NS</td>
<td>0.0008</td>
<td>S (20.7)</td>
<td>0.740</td>
</tr>
<tr><td><italic>b</italic>
</td>
<td>16.4</td>
<td colspan="2">−7.0 to 36.0</td>
<td>NS</td>
<td></td>
<td></td>
<td></td>
</tr>
</tbody>
</table>
<table-wrap-foot><p><italic>CI</italic>
confidence interval, <italic>R</italic>
<sup>2</sup>
is coefficient of determination, <italic>S</italic>
significant (<italic>p</italic>
< 0.01), <italic>NS</italic>
not significant (<italic>p</italic>
> 0.05), <italic>B</italic>
borderline significant (0.01 < <italic>p</italic>
< 0.05)</p>
<p><sup>a</sup>
Nonlinear, 0.05 < <italic>p</italic>
< 0.1 from cusum. However, <italic>R</italic>
<sup>2</sup>
(0.928) and the Spearman rank coefficient squared (0.929) are similar</p>
</table-wrap-foot>
</table-wrap>
<fig id="Fig2"><label>Fig. 2</label>
<caption><p><bold>a</bold>
The Tk-GV plasma clearance method, <italic>CL</italic>
(Tk-GV), was applied to 4- and 24-h data and plotted. When this is fit with a Passing-Bablok regression line, there is no significant departure from unity slope and no significant <italic>y</italic>
-intercept. <bold>b</bold>
The 4-h data <italic>CL</italic>
(E1), plotted against 24-h <italic>CL</italic>
(E1), has a slope significantly >1, i.e. the 24-h <italic>CL</italic>
(E1) values are overestimated by 4-h <italic>CL</italic>
(E1) calculations. <bold>c</bold>
<italic>CL</italic>
(Tk-GV) for 12- and 24-h data and no significant departure from a slope of 1 and an intercept of zero. <bold>d</bold>
Even for 24-h data, <inline-formula id="IEq12"><alternatives><tex-math id="M4">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ CL\left( {{\text{Tk - GV}}} \right)\, \ne \,CL\left( {\text{E1}} \right) $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq12.gif"></inline-graphic>
</alternatives>
</inline-formula>
, with both the slope <inline-formula id="IEq13"><alternatives><tex-math id="M5">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \ne 1 $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq13.gif"></inline-graphic>
</alternatives>
</inline-formula>
and the intercept <inline-formula id="IEq14"><alternatives><tex-math id="M6">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \ne 0 $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq14.gif"></inline-graphic>
</alternatives>
</inline-formula>
of the Passing-Bablok linear fit</p>
</caption>
<graphic xlink:href="259_2011_1887_Fig2_HTML" id="MO3"></graphic>
</fig>
</p>
<p>Significantly poorer results were obtained from the 4-h versus 24-h E1 fitting (see Table <xref rid="Tab2" ref-type="table">2</xref>
). Figure <xref rid="Fig2" ref-type="fig">2</xref>
b shows these results plotted with a regression slope of 1.56, i.e. the 4-h <italic>CL</italic>
-values are significantly larger than the 24-h <italic>CL</italic>
-values with a very significant (<italic>p</italic>
= 0.0001) mean difference of 32.1 ml/min. However, the R<sup>2</sup>
value for this plot is 0.738, whereas the R<sup>2</sup>
value (or Spearman rank sum squared, <italic>rs</italic>
<sup>2</sup>
) for 4-h <italic>CL</italic>
(Tk-GV) versus <italic>CL</italic>
(E1) data would be significantly better, i.e. 0.928 (<italic>rs</italic>
<sup>2</sup>
= 0.929). That the R<sup>2</sup>
values of 4- versus 24-h <italic>CL</italic>
(Tk-GV) and 4- versus 24-h <italic>CL</italic>
(E1) are equal to 3 decimal places (i.e. 0.738), and 4-h <italic>CL</italic>
(Tk-GV) and <italic>CL</italic>
(E1) agree better with each other (0.928) than either 4-h <italic>CL</italic>
agrees with its corresponding 24-h result (i.e. 0.738) appears to confirm error propagation as explaining much of the variance for predicting 24-h <italic>CL</italic>
from 4-h data.</p>
<p>When the sample times included in the fit are increased to 12 h, then <italic>CL</italic>
(E1-12-h) > <italic>CL</italic>
(E1-24-h) in all cases. The paired <italic>t</italic>
-test for this is highly significant (<italic>p</italic>
= 0.0002, mean difference 9.1 ml/min), with an improved R<sup>2</sup>
(0.951), (see Table <xref rid="Tab2" ref-type="table">2</xref>
). From the 24-h data (Fig. <xref rid="Fig1" ref-type="fig">1</xref>
d), the plot of <italic>CL</italic>
(Tk-GV) versus <italic>CL</italic>
(E1) has a significantly biased slope (0.85) and borderline significant intercept (5.97 ml/min, 0.01 < <italic>p</italic>
< 0.05). Note that from Table <xref rid="Tab2" ref-type="table">2</xref>
, the <italic>t</italic>
-test <italic>p</italic>
= 0.8386 mean difference of 0.4 ml/min is not significant. In other words, on the average, <italic>CL</italic>
(Tk-GV-24-h) values did not differ from <italic>CL</italic>
(E1-24-h) values, but, for low <italic>CL</italic>
-values, <italic>CL</italic>
(E1) > <italic>CL</italic>
(Tk-GV), and for high <italic>CL</italic>
-values, <italic>CL</italic>
(E1) < <italic>CL</italic>
(Tk-GV).</p>
<p>The slope intercept (E1) 2- and 4-h plasma sample method results follow. The Passing-Bablok regression slope of <italic>CL</italic>
(E1-2- & 4-h) versus <italic>CL</italic>
(Tk-GV-24-h) has no bias detected for slope (1.01) and intercept (16.4 ml/min). However, the intercept is larger than for <italic>CL</italic>
(Tk-GV-4-h) versus <italic>CL</italic>
(Tk-GV-24-h) (6.7 ml/min). Indeed, paired <italic>t</italic>
testing (Table <xref rid="Tab2" ref-type="table">2</xref>
) shows the mean bias (20.7 ml/min) to be highly significant (<italic>p</italic>
= 0.0008) for <italic>CL</italic>
(E1-2- & 4-h) versus <italic>CL</italic>
(Tk-GV-24-h). This compares to the more moderate mean bias (8.6 ml/min) of borderline significance (<italic>p</italic>
= 0.0337) for the <italic>CL</italic>
(Tk-GV-4-h) versus <italic>CL</italic>
(Tk-GV-24-h) plot. The Passing-Bablok regression of <italic>CL</italic>
(E1-24-h) (vertical axis) versus <italic>CL</italic>
(E1-2- & 4-h) (horizontal axis) has a significantly lesser slope (<italic>p</italic>
< 0.01) than 1, i.e. 0.65. However, both <italic>CL</italic>
(E1-24-h) and <italic>CL</italic>
(E1-2- & 4-h) have <italic>r</italic>
-values for fitting the logarithms of concentrations that are very significantly poorer (<italic>z</italic>
-statistic <italic>t</italic>
-test <italic>p</italic>
< 0.0001) than <italic>CL</italic>
(Tk-GV-24-h).</p>
</sec>
<sec id="Sec4"><title>Discussion</title>
<p>The Tk-GV models have a less obvious dependence on sampling-time effects than E1 fits. They also have less bias for extrapolation and fit the data very significantly better (<italic>p</italic>
< 0.0001). Of some concern is that the predictions of 24-h <italic>CL</italic>
from 4-h <italic>CL</italic>
data were imprecise (R<sup>2</sup>
= 0.738). This lower than optimal coefficient of determination is partly due to a restricted range of <italic>CL</italic>
-values from examining low renal function patients only, may partly be due to biological fluctuation of <italic>CL</italic>
in time and is due to the propagation of errors of measured concentrations, which in turn is due to errors in sample measurements, e.g. see [<xref ref-type="bibr" rid="CR2">2</xref>
]. In early prospective clinical experience, <italic>CL</italic>
(Tk-GV) has been implemented for four-sample 4-h <sup>99m</sup>
Tc-DTPA, with some notable features being that the ease of use, stability of the solutions and that the error calculations help to improve sampling technique and inspire confidence in the results. Some of the sources for error are easily addressed. For example, extravasation in the injection site is reduced by applying the flush-bolus technique and checked for by imaging of the injection site [<xref ref-type="bibr" rid="CR12">12</xref>
]. Whatever the causes for error, use of Tk-GV has the advantage of increased accuracy for extrapolated <italic>CL</italic>
(defined as lack of bias) while having precision (defined as R<sup>2</sup>
for <italic>CL</italic>
extrapolation) as good as or better than that of E1 fitting.</p>
<p>Fractal structures have been postulated as causing a gamma variate (GV) shape for inert substance <italic>CL</italic>
-curves [<xref ref-type="bibr" rid="CR13">13</xref>
]. However, it is only when a GV is found using a <italic>CL</italic>
relative-error minimizing technique that GV (1) accurately predicts the late-time concentration and (2) consistently forms a model with a dynamic concentration gradient established by renal function. The results here suggest that the 24-h <italic>CL</italic>
(Tk-GV) values are a standard against which to judge clearance, because of the very good quality of all of the curve fits (median <italic>r</italic>
= 0.998) (see Fig. <xref rid="Fig1" ref-type="fig">1</xref>
and Table <xref rid="Tab1" ref-type="table">1</xref>
). Moreover, prediction of 24-h <italic>CL</italic>
(Tk-GV) from 4- and 12-h <italic>CL</italic>
(Tk-GV) was good, suggesting that 24-h Tk-GV fits should extrapolate concentration correctly at even later times. In particular for six statistical tests, <italic>CL</italic>
(Tk-GV) extrapolation was unbiased for four tests and borderline biased (i.e. 0.01 < <italic>p</italic>
< 0.05) for two, whereas <italic>CL</italic>
(E1) was significantly biased in three of six tests. The 4-h <italic>CL</italic>
(Tk-GV) versus <italic>CL</italic>
(E1) (Table <xref rid="Tab2" ref-type="table">2</xref>
) is somewhat nonlinear (cusum probability, 0.05 < <italic>p</italic>
< 0.1 [<xref ref-type="bibr" rid="CR7">7</xref>
]) and the clinically used slope intercept method with 2- and 4-h plasma samples differs from 24-h <italic>CL</italic>
(Tk-GV) by a constant (mean difference 20.7 ml/min, <italic>p</italic>
= 0.0008), so there is no “quick-fix” linear transformation of <italic>CL</italic>
(E1) values to estimate <italic>CL</italic>
(Tk-GV). This is discordant with the literature suggesting that <italic>CL</italic>
(E1-2- & 4-h) may be corrected by multiplying by 0.87, i.e. the Chantler correction factor, and does not suggest subtracting a constant [<xref ref-type="bibr" rid="CR2">2</xref>
]. However, it is easy enough to see how such an impression could arise as the Passing-Bablok regression of <italic>CL</italic>
(E1-24-h) (vertical axis) versus <italic>CL</italic>
(E1-2- & 4-h) (horizontal axis) only has a significant slope (0.65, <italic>p</italic>
< 0.01 of slope =1) and an insignificant intercept. However, both <italic>CL</italic>
(E1-24-h) and <italic>CL</italic>
(E1-2- & 4-h) have <italic>r</italic>
-values for fitting the logarithms of concentrations that are very significantly poorer (<italic>z</italic>
-statistic <italic>t</italic>
-test <italic>p</italic>
< 0.0001) than <italic>CL</italic>
(Tk-GV-24-h). Thus, subtraction of a 16.4 (or 20.7) ml/min constant from the <italic>CL</italic>
(E1-2- & 4-h) values may be the more accurate correction, at least for our patients with ascites [and if one accepts that <italic>CL</italic>
(Tk-GV-24-h) is a gold standard].</p>
<p>A final caveat is that eight of ten patients had ascites. Ascites has been associated with increased <italic>CL</italic>
for compartmental models by 16–20 ml/min [<xref ref-type="bibr" rid="CR4">4</xref>
, <xref ref-type="bibr" rid="CR6">6</xref>
, <xref ref-type="bibr" rid="CR9">9</xref>
] and is considered by some to relatively contraindicate <italic>CL</italic>
testing due to gross inaccuracy [<xref ref-type="bibr" rid="CR2">2</xref>
]. However, ascites seems to have had little effect on the <italic>CL</italic>
(Tk-GV)-values here. This is significant and it needs to be emphasized that the Tk-GV method is valid in patients with low <italic>GFR</italic>
and with ascites. This comes about as the Tk-GV volumes of distribution (<italic>V</italic>
) (and <italic>V/W</italic>
, Table <xref rid="Tab1" ref-type="table">1</xref>
) are larger than usual (e.g. see [<xref ref-type="bibr" rid="CR2">2</xref>
]). However, that is an expected result for ascites (also see [<xref ref-type="bibr" rid="CR8">8</xref>
]), as follows. For a bolus clearance model, <italic>MRT=V/CL</italic>
, where <italic>MRT</italic>
is the mean residence time. If in that equation <italic>CL</italic>
is held constant, then for a longer residence time, <italic>V</italic>
must be larger, and for a good model, that is an expected outcome of ascites. <italic>V/W</italic>
is our measure of ascites from the 24-h Tk-GV model. The correlation of <italic>r</italic>
-values from E1 24-h data fits with <italic>V/W</italic>
(Table <xref rid="Tab1" ref-type="table">1</xref>
) is significantly negative, <inline-formula id="IEq15"><alternatives><tex-math id="M7">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ r\left[ {r\left( {{\text{E}}1} \right),{ }V\,/\,W} \right] = - 0.823 $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq15.gif"></inline-graphic>
</alternatives>
</inline-formula>
, <italic>p</italic>
=0.0034. However, the correlation of <italic>r</italic>
-values from 24-h Tk-GV fits to the data with <italic>V/W</italic>
from the same Tk-GV model is insignificant, <inline-formula id="IEq16"><alternatives><tex-math id="M8">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ r\left[ {r\left( {{\text{Tk - GV}}} \right),{ }V/W} \right] = - 0.084 $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq16.gif"></inline-graphic>
</alternatives>
</inline-formula>
, <italic>p</italic>
= 0.8165. For <inline-formula id="IEq17"><alternatives><tex-math id="M9">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ r\left[ {r\left( {{\text{E}}1} \right),{ }V/W} \right] $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq17.gif"></inline-graphic>
</alternatives>
</inline-formula>
, we have not caused a correlation by reusing parameters, because the only sharing between models is from the characteristics of the data themselves. For <inline-formula id="IEq18"><alternatives><tex-math id="M10">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ r\left[ {r\left( {{\text{Tk - GV}}} \right),{ }V/W} \right] $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq18.gif"></inline-graphic>
</alternatives>
</inline-formula>
, we have reused parameters, and despite this, there is no correlation. We conclude that using <italic>V/W</italic>
(from 24-h Tk-GV) as an index of the severity of ascites, E1 fits degrade as ascites increases, and Tk-GV fit quality is not measurably affected. This finding confirms the observation that use of E1 (or E1 extrapolation) is problematic in ascites [<xref ref-type="bibr" rid="CR4">4</xref>
, <xref ref-type="bibr" rid="CR6">6</xref>
, <xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR9">9</xref>
] and presents a solution to this problem using Tk-GV. To inspect this in detail, let us examine the GV rate equation,<disp-formula id="Equ1"><label>1</label>
<alternatives><tex-math id="M11">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \frac{1}{{C(t)}}\frac{{dC(t)}}{{dt}} = \frac{{\alpha - 1}}{t} - \beta ..., $$\end{document}</tex-math>
<graphic xlink:href="259_2011_1887_Article_Equ1.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
where (<italic>α</italic>
−1)/<italic>t</italic>
is definite negative, is associated with tissue clearance, and gets slower and slower in time (it is proportional to 1/<italic>t</italic>
). <italic>CL</italic>
><italic>GFR</italic>
implies <italic>α</italic>
≤1, where <italic>α</italic>
>1 is a common [<xref ref-type="bibr" rid="CR8">8</xref>
], but ill-posed solution for <italic>CL</italic>
for GV fitting [<xref ref-type="bibr" rid="CR13">13</xref>
]. Here, using Tk-GV fitting, the entire range of <italic>α</italic>
for our 24-h data is from 0.62 to 0.88, i.e. <italic>α</italic>
is well behaved. So <italic>α</italic>
does not change by much to accompany an increase of fluid volume. Here, a low <italic>α</italic>
=0.67 coincides with the most inflated volumes of distribution as a percentage of patient mass 59.9% (case 8, Table <xref rid="Tab1" ref-type="table">1</xref>
).</p>
<p>Normally slow or pathologically slow clearance areas in the body include ascites or abdominal fluid, some adenomas and neoplasias, some chemotherapy, pleural or pericardial fluid, lymphatic circulation and lymphoedema, cerebrospinal fluid circulation, synovial fluid, circulation in aqueous humour, abscess, cellulitis, fibrous dysplasia, etc. The entire spectrum of slow tissue marker leakage rates is included in GV rate Eq. <xref rid="Equ1" ref-type="">1</xref>
because for all <italic>t</italic>
>(1−<italic>α</italic>
)/<italic>β</italic>
, the rates of losses of marker into tissue are slower than renal loss. That is, tissue marker leakage rates converging to zero (i.e. to an infinite half-life) are an integral part of the GV model. By contrast, an infinite sum of exponential terms would be needed to approximate the effect of the temporal spectrum of the single, gradually decreasing rate term of the GV model.</p>
<p>Why only compare Tk-GV to E1 fitting? Why not then test E2, E3, E4 or E<italic>n</italic>
(i.e. sums of exponential terms) in general? The reason is that a robust model is strongly desirable, and E(<italic>n</italic>
> 1) is not stable for regression analysis unless the number of samples is substantially greater than <italic>n</italic>
[<xref ref-type="bibr" rid="CR10">10</xref>
, <xref ref-type="bibr" rid="CR13">13</xref>
]. Although methods exist for stabilizing E2 models (e.g. see [<xref ref-type="bibr" rid="CR11">11</xref>
]), in general for regression, E<italic>n</italic>
is not robust, indeed ordinary least-squares GV function fitting is also not robust. However, GV fitting can be made robust by regularization (as in Tk-GV [<xref ref-type="bibr" rid="CR13">13</xref>
]). It is at best tedious and at worst impossible to obtain four (or more) stable parameter estimates with an unstable E2 (or higher) regression model. By contrast, a single attempt suffices to obtain the globally converged Tk-GV model’s three parameters, with superior quality curve fits, better extrapolation to late-time concentrations from fits to early concentration data and using fewer samples [<xref ref-type="bibr" rid="CR13">13</xref>
]. For practical measurements, one strongly desires either a simple robust model (e.g. E1 with correction factors) or another robust model (e.g. Tk-GV with no need of correction factors). However, as shown in the results above (Table <xref rid="Tab2" ref-type="table">2</xref>
) the magnitude and the type of the correction factors needed for E1 <italic>CL</italic>
-values are a function of the times used to obtain the samples. Even if it were clear how to correct E1 <italic>CL</italic>
-values (and it is not), any such timing correction formulae would have empirical coefficients. For example, these formulae would not apply to both small and large mammals without allometric rescaling. Moreover, E1 fit quality degrades in ascites, rendering the utility of the model questionable. In summary, making “corrections” to an inherently flawed exponential term model is not optimal. The better approach using Tk-GV better fits the observed real-life data, without corrections needed.</p>
<p>Is 24-h a long enough data collection time to determine <italic>CL</italic>
in renal failure (<italic>CL</italic>
< 15 ml/min)? The results suggest that if we use E1 functions for fitting, the answer would be no. It is common to plot logarithm concentration versus time, which would linearize ln<italic>C</italic>
(<italic>t</italic>
)≈−<italic>at</italic>
+ <italic>b</italic>
, i.e. an E1 function. In fact, note in Fig. <xref rid="Fig1" ref-type="fig">1</xref>
that <italic>C</italic>
(<italic>t</italic>
) at early times is approximately linear in ln<italic>t</italic>
, i.e. <italic>C</italic>
(<italic>t</italic>
)≈−<italic>a</italic>
ln<italic>t</italic>
+ <italic>b</italic>
. This latter plot type helps visually confirm the correlation results, in that E1 functions have the wrong shape and do not follow the concentration curves well. Adding more exponential terms does not solve the curve shape problem. Let us define late-time <italic>C</italic>
(<italic>t</italic>
) as being after the concentration curve becomes nonlinear in ln(<italic>t</italic>
), i.e. asymptotic to the <italic>C</italic>
(<italic>t</italic>
)=0 axis. Figure <xref rid="Fig1" ref-type="fig">1</xref>
b demonstrates that when an exponential fit does not include late-time concentrations, the late-time concentrations are underestimated. This explains the overestimation of <italic>CL</italic>
(from underestimation of <italic>AUC</italic>
) for (all choices of) E1 early-time extrapolations to later-time fits in Table <xref rid="Tab2" ref-type="table">2</xref>
. However, with the better-fitting Tk-GV functions, the extrapolation errors are not so severe, as suggested by extrapolation testing of early data. Using Tk-GV for the lowest <italic>CL</italic>
value (11.4 ml/min, case 10, Fig. <xref rid="Fig1" ref-type="fig">1</xref>
b), one predicts that it would take 6 days of data collection for the concentration to go below 1% of the concentration measured in the 5-min sample. Adding a “terminal” exponential to the 1-day data collection would underestimate the concentration at day 6, as exponentials go to zero too quickly to mimic late concentration. Adding another early exponential term to the Fig. <xref rid="Fig1" ref-type="fig">1</xref>
E1 fits would not remove the exponential fit assumption of a constant concentration at time zero causing underestimation of the early concentrations of actual data. On the other hand, finding the <italic>AUC</italic>
is not especially challenging when the data collected include all the concentrations that contribute to them, and any arbitrary interpolative curve fitting, or numerical integration [<xref ref-type="bibr" rid="CR5">5</xref>
], will suffice for that purpose. Especially as relates to early-time only data collections, proper extrapolation is not as arbitrary as interpolation and must be validated by extrapolation testing, done here for Tk-GV. Thus, if a 4- or 12-h test is desired for finding <italic>CL</italic>
, Tk-GV fitting appears to provide an accurate (i.e. relatively unbiased) method of predicting <italic>CL</italic>
for 24-h, even in the presence of ascites.</p>
</sec>
<sec id="Sec5"><title>Conclusion</title>
<p><italic>CL</italic>
(E1) values were dependent on sample time and the quality of fit was poor and degraded with increasing ascites, consistent with current opinion that <italic>CL</italic>
(E1) is contraindicated in ascitic patients. <italic>CL</italic>
(Tk-GV) was relatively more accurate and the good quality of fit was unaffected by ascites. <italic>CL</italic>
(Tk-GV) was the preferred method for the accurate calculation of <italic>CL</italic>
and was useful despite liver failure and ascites.</p>
</sec>
</body>
<back><ack><p><bold>Conflicts of interest</bold>
None.</p>
<p><bold>Open Access</bold>
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.</p>
</ack>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Calvert</surname>
<given-names>AH</given-names>
</name>
<name><surname>Egorin</surname>
<given-names>MJ</given-names>
</name>
</person-group>
<article-title>Carboplatin dosing formulae: gender bias and the use of creatinine-based methodologies</article-title>
<source>Eur J Cancer</source>
<year>2002</year>
<volume>38</volume>
<fpage>11</fpage>
<lpage>16</lpage>
<pub-id pub-id-type="doi">10.1016/S0959-8049(01)00340-9</pub-id>
<pub-id pub-id-type="pmid">11750834</pub-id>
</mixed-citation>
</ref>
<ref id="CR2"><label>2.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fleming</surname>
<given-names>JS</given-names>
</name>
<name><surname>Zivanovic</surname>
<given-names>MA</given-names>
</name>
<name><surname>Blake</surname>
<given-names>GM</given-names>
</name>
<name><surname>Burniston</surname>
<given-names>M</given-names>
</name>
<name><surname>Cosgriff</surname>
<given-names>PS</given-names>
</name>
<collab>British Nuclear Medicine Society</collab>
</person-group>
<article-title>Guidelines for the measurement of glomerular filtration rate using plasma sampling</article-title>
<source>Nucl Med Commun</source>
<year>2004</year>
<volume>25</volume>
<fpage>759</fpage>
<lpage>769</lpage>
<pub-id pub-id-type="doi">10.1097/01.mnm.0000136715.71820.4a</pub-id>
<pub-id pub-id-type="pmid">15266169</pub-id>
</mixed-citation>
</ref>
<ref id="CR3"><label>3.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Florijn</surname>
<given-names>KW</given-names>
</name>
<name><surname>Barendregt</surname>
<given-names>JNM</given-names>
</name>
<name><surname>Lentjes</surname>
<given-names>EGW</given-names>
</name>
<name><surname>Dam</surname>
<given-names>W</given-names>
</name>
<name><surname>Prodjosudjadi</surname>
<given-names>W</given-names>
</name>
<name><surname>Saase</surname>
<given-names>JL</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Glomerular filtration rate measurement by “single-shot” injection of inulin</article-title>
<source>Kidney Int</source>
<year>1994</year>
<volume>46</volume>
<fpage>252</fpage>
<lpage>259</lpage>
<pub-id pub-id-type="doi">10.1038/ki.1994.267</pub-id>
<pub-id pub-id-type="pmid">7933844</pub-id>
</mixed-citation>
</ref>
<ref id="CR4"><label>4.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Henriksen</surname>
<given-names>JH</given-names>
</name>
<name><surname>Brøchner-Mortensen</surname>
<given-names>J</given-names>
</name>
<name><surname>Malchow-Møller</surname>
<given-names>A</given-names>
</name>
<name><surname>Schlichting</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Over-estimation of glomerular filtration rate by single injection [51Cr]EDTA plasma clearance determination in patients with ascites</article-title>
<source>Scand J Clin Lab Invest</source>
<year>1980</year>
<volume>40</volume>
<fpage>279</fpage>
<lpage>284</lpage>
<pub-id pub-id-type="doi">10.3109/00365518009095579</pub-id>
<pub-id pub-id-type="pmid">6777855</pub-id>
</mixed-citation>
</ref>
<ref id="CR5"><label>5.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moore</surname>
<given-names>AE</given-names>
</name>
<name><surname>Park-Holohan</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Blake</surname>
<given-names>GM</given-names>
</name>
<name><surname>Fogelman</surname>
<given-names>I</given-names>
</name>
</person-group>
<article-title>Conventional measurements of GFR using 51Cr-EDTA overestimate true renal clearance by 10 percent</article-title>
<source>Eur J Nucl Med Mol Imaging</source>
<year>2003</year>
<volume>30</volume>
<fpage>4</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.1007/s00259-002-1007-y</pub-id>
<pub-id pub-id-type="pmid">12483403</pub-id>
</mixed-citation>
</ref>
<ref id="CR6"><label>6.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nielsen</surname>
<given-names>SS</given-names>
</name>
<name><surname>Havsteen</surname>
<given-names>H</given-names>
</name>
<name><surname>Petersen</surname>
<given-names>LK</given-names>
</name>
<name><surname>Nielsen</surname>
<given-names>LE</given-names>
</name>
<name><surname>Rehling</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Patients with ovarian cancer have elevated (51)Cr-EDTA plasma clearance early post-operatively</article-title>
<source>Nucl Med Commun</source>
<year>2002</year>
<volume>23</volume>
<fpage>917</fpage>
<lpage>920</lpage>
<pub-id pub-id-type="doi">10.1097/00006231-200209000-00016</pub-id>
<pub-id pub-id-type="pmid">12195097</pub-id>
</mixed-citation>
</ref>
<ref id="CR7"><label>7.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Passing</surname>
<given-names>H</given-names>
</name>
<name><surname>Bablok</surname>
<given-names>W</given-names>
</name>
</person-group>
<article-title>A new biometrical procedure for testing the equality of measurements from two different analytical methods. Application of linear regression procedures for method comparison studies in clinical chemistry, Part I</article-title>
<source>J Clin Chem Clin Biochem</source>
<year>1983</year>
<volume>21</volume>
<fpage>709</fpage>
<lpage>720</lpage>
<pub-id pub-id-type="pmid">6655447</pub-id>
</mixed-citation>
</ref>
<ref id="CR8"><label>8.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Perkinson</surname>
<given-names>AS</given-names>
</name>
<name><surname>Evans</surname>
<given-names>CJ</given-names>
</name>
<name><surname>Burniston</surname>
<given-names>MT</given-names>
</name>
<name><surname>Smye</surname>
<given-names>SW</given-names>
</name>
</person-group>
<article-title>The effect of improved modelling of plasma clearance in paediatric patients with expanded body spaces on estimation of the glomerular filtration rate</article-title>
<source>Physiol Meas</source>
<year>2010</year>
<volume>31</volume>
<fpage>183</fpage>
<lpage>192</lpage>
<pub-id pub-id-type="doi">10.1088/0967-3334/31/2/005</pub-id>
<pub-id pub-id-type="pmid">20016117</pub-id>
</mixed-citation>
</ref>
<ref id="CR9"><label>9.</label>
<mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Rehling</surname>
<given-names>M</given-names>
</name>
<name><surname>Stadeager</surname>
<given-names>C</given-names>
</name>
<name><surname>Henriksen</surname>
<given-names>JH</given-names>
</name>
<name><surname>Siemsen</surname>
<given-names>O</given-names>
</name>
<name><surname>Krintel</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Malchow-Møller</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<person-group person-group-type="editor"><name><surname>Thomsen</surname>
<given-names>HS</given-names>
</name>
<name><surname>Nally</surname>
<given-names>JP</given-names>
<suffix>Jr</suffix>
</name>
<name><surname>Britton</surname>
<given-names>K</given-names>
</name>
<name><surname>Frùkiñr</surname>
<given-names>J</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Measurement of glomerular filtration rate in patients with ascites</article-title>
<source>Radionuclides in nephrourology</source>
<year>1998</year>
<publisher-loc>Copenhagen</publisher-loc>
<publisher-name>FADL</publisher-name>
<fpage>114</fpage>
<lpage>118</lpage>
</mixed-citation>
</ref>
<ref id="CR10"><label>10.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Russell</surname>
<given-names>CD</given-names>
</name>
</person-group>
<article-title>Optimum sample times for single-injection, multisample renal clearance methods</article-title>
<source>J Nucl Med</source>
<year>1993</year>
<volume>34</volume>
<fpage>1761</fpage>
<lpage>1765</lpage>
<pub-id pub-id-type="pmid">8410295</pub-id>
</mixed-citation>
</ref>
<ref id="CR11"><label>11.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Russell</surname>
<given-names>CD</given-names>
</name>
<name><surname>Taylor</surname>
<given-names>AT</given-names>
</name>
<name><surname>Dubovsky</surname>
<given-names>EV</given-names>
</name>
</person-group>
<article-title>A Bayesian regression model for plasma clearance</article-title>
<source>J Nucl Med</source>
<year>2002</year>
<volume>43</volume>
<fpage>762</fpage>
<lpage>766</lpage>
<pub-id pub-id-type="pmid">12050320</pub-id>
</mixed-citation>
</ref>
<ref id="CR12"><label>12.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wesolowski</surname>
<given-names>CA</given-names>
</name>
<name><surname>Hogendoorn</surname>
<given-names>P</given-names>
</name>
<name><surname>Vandierendonck</surname>
<given-names>R</given-names>
</name>
<name><surname>Driedger</surname>
<given-names>AA</given-names>
</name>
</person-group>
<article-title>Bolus injections of measured amounts of radioactivity</article-title>
<source>J Nucl Med Technol</source>
<year>1988</year>
<volume>16</volume>
<fpage>1</fpage>
<lpage>4</lpage>
</mixed-citation>
</ref>
<ref id="CR13"><label>13.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wesolowski</surname>
<given-names>CA</given-names>
</name>
<name><surname>Puetter</surname>
<given-names>RC</given-names>
</name>
<name><surname>Ling</surname>
<given-names>L</given-names>
</name>
<name><surname>Babyn</surname>
<given-names>PS</given-names>
</name>
</person-group>
<article-title>Tikhonov adaptively regularized gamma variate fitting to assess plasma clearance of inert renal markers</article-title>
<source>J Pharmacokinet Pharmacodyn</source>
<year>2010</year>
<volume>37</volume>
<fpage>435</fpage>
<lpage>474</lpage>
<pub-id pub-id-type="doi">10.1007/s10928-010-9167-z</pub-id>
<pub-id pub-id-type="pmid">20865304</pub-id>
</mixed-citation>
</ref>
<ref id="CR14"><label>14.</label>
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xirouchakis</surname>
<given-names>E</given-names>
</name>
<name><surname>Marelli</surname>
<given-names>L</given-names>
</name>
<name><surname>Cholongitas</surname>
<given-names>E</given-names>
</name>
<name><surname>Manousou</surname>
<given-names>P</given-names>
</name>
<name><surname>Calvaruso</surname>
<given-names>V</given-names>
</name>
<name><surname>Pleguezuelo</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Comparison of cystatin C and creatinine-based glomerular filtration rate formulas with 51Cr-EDTA clearance in patients with cirrhosis</article-title>
<source>Clin J Am Soc Nephrol</source>
<year>2011</year>
<volume>6</volume>
<fpage>84</fpage>
<lpage>92</lpage>
<pub-id pub-id-type="doi">10.2215/CJN.03400410</pub-id>
<pub-id pub-id-type="pmid">20829419</pub-id>
</mixed-citation>
</ref>
</ref-list>
<app-group><app id="App1"><sec id="Sec6"><title>Appendix</title>
<p>Tikhonov regularization</p>
<p>For an overdetermined system of linear equations, <bold>Ax = b</bold>
, the Tikhonov regularization (Tk) of this problem introduces the penalty function Γ<bold>x</bold>
and seeks to find a solution that minimizes <inline-formula id="IEq19"><alternatives><tex-math id="M12">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\left\| {Ax - b} \right\|^2} + {\left\| {\Gamma x} \right\|^2} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq19.gif"></inline-graphic>
</alternatives>
</inline-formula>
. This latter is the square of a norm of the residuals, <inline-formula id="IEq20"><alternatives><tex-math id="M13">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\left\| {Ax - b} \right\|^2} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq20.gif"></inline-graphic>
</alternatives>
</inline-formula>
, plus the square of a norm of the product of the Tikhonov matrix, Γ, with the <bold>x</bold>
fit parameters (unknowns). The more general Γ<sup>T</sup>
Γ regularizing term is often, as it is here, replaced by <italic>λ</italic>
<bold>I</bold>
, where <bold>I</bold>
is the identity matrix, and <italic>λ</italic>
is a Lagrange (i.e. constraint) multiplier, also commonly called the shrinkage, Tikhonov or damping factor.</p>
<p>A constraint on ln<italic>K</italic>
</p>
<p>A most common constraint for regression is to require the fit function to pass through the data mean point [aka the centroid, (<inline-formula id="IEq21"><alternatives><tex-math id="M14">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \bar{x},\bar{y} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq21.gif"></inline-graphic>
</alternatives>
</inline-formula>
)]. Because the logarithm of concentrations is the more homoscedastic quantity, it is common to fit the logarithms of marker concentrations rather than the concentrations themselves. Thus, for the Tk-GV method, the GV function is written <inline-formula id="IEq22"><alternatives><tex-math id="M15">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \ln C = \ln K + \left( {\alpha - 1} \right)\ln t - \beta t $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq22.gif"></inline-graphic>
</alternatives>
</inline-formula>
, where the constant term <inline-formula id="IEq23"><alternatives><tex-math id="M16">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \ln K $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq23.gif"></inline-graphic>
</alternatives>
</inline-formula>
need not be independent, but can be determined from the other fit parameters using a mean value constraint. Taking averages over the data<disp-formula id="Equ2"><label>2</label>
<alternatives><tex-math id="M17">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \ln K = \overline {\ln C(t)} - \left( {\alpha - 1} \right)\overline {\ln t} + \beta \bar{t} = \bar{b} - \left( {\alpha - 1} \right){\bar{a}_1} + \beta {\bar{a}_2}..., $$\end{document}</tex-math>
<graphic xlink:href="259_2011_1887_Article_Equ2.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
such that <inline-formula id="IEq24"><alternatives><tex-math id="M18">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \bar{b} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq24.gif"></inline-graphic>
</alternatives>
</inline-formula>
, <inline-formula id="IEq25"><alternatives><tex-math id="M19">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\bar{a}_1} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq25.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq26"><alternatives><tex-math id="M20">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\bar{a}_2} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq26.gif"></inline-graphic>
</alternatives>
</inline-formula>
are data constants, where <inline-formula id="IEq27"><alternatives><tex-math id="M21">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \bar{b} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq27.gif"></inline-graphic>
</alternatives>
</inline-formula>
is the mean value of the logarithms of the concentrations, <inline-formula id="IEq28"><alternatives><tex-math id="M22">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\bar{a}_1} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq28.gif"></inline-graphic>
</alternatives>
</inline-formula>
is the mean of the logarithms of the sample times and <inline-formula id="IEq29"><alternatives><tex-math id="M23">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\bar{a}_2} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq29.gif"></inline-graphic>
</alternatives>
</inline-formula>
is the mean of the sample times. Then, Eq. <xref rid="Equ2" ref-type="">2</xref>
is used to remove <italic>K</italic>
from the formula for <italic>CL</italic>
, and an expression is derived for the errors in <italic>CL</italic>
with only <italic>α</italic>
and <italic>β</italic>
as independent parameters, as follows<disp-formula id="Equ3"><label>3</label>
<alternatives><tex-math id="M24">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ CL = \frac{D}{{AUC}} = \frac{{D{\beta^{\alpha }}}}{{K\Gamma \left( \alpha \right)}} = \frac{{D{\beta^{\alpha }}}}{{\exp \left[ {\bar{b} - \left( {\alpha - 1} \right){{\bar{a}}_1} + \beta {{\bar{a}}_2}} \right]\Gamma \left( \alpha \right)}}..., $$\end{document}</tex-math>
<graphic xlink:href="259_2011_1887_Article_Equ3.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
where the gamma function, <inline-formula id="IEq30"><alternatives><tex-math id="M25">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \Gamma (x) = \int_0^{\infty } {{t^{{x - 1}}}{e^{{ - t}}}dt} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq30.gif"></inline-graphic>
</alternatives>
</inline-formula>
.</p>
<p>Error propagation and the error adaptive equation</p>
<p>One applies error propagation (Taylor series expansion) to the right-hand <italic>CL</italic>
expression of Eq. <xref rid="Equ3" ref-type="">3</xref>
with respect to <italic>α</italic>
and <italic>β</italic>
yielding<disp-formula id="Equ4"><label>4</label>
<alternatives><tex-math id="M26">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\left( {\frac{{{s_{{CL}}}}}{{CL}}} \right)^2} = s_{\alpha }^2{\left( {{{\bar{a}}_1} + \ln \beta - \Psi \left( \alpha \right)} \right)^2} + s_{\beta }^2{\left( {\frac{\alpha }{\beta } - {{\bar{a}}_2}} \right)^2} + 2{s_{{\alpha \beta }}}\left( {{{\bar{a}}_1} + \ln \beta - \Psi \left( \alpha \right)} \right)\left( {\frac{\alpha }{\beta } - {{\bar{a}}_2}} \right)... $$\end{document}</tex-math>
<graphic xlink:href="259_2011_1887_Article_Equ4.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
where Ψ(<italic>α</italic>
) is the digamma function of <italic>α</italic>
and <inline-formula id="IEq31"><alternatives><tex-math id="M27">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \Psi \left( \alpha \right) = d\left[ {\ln \Gamma \left( \alpha \right)} \right]/d\alpha = \Gamma \prime \left( \alpha \right)/\Gamma \left( \alpha \right) $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq31.gif"></inline-graphic>
</alternatives>
</inline-formula>
, the subscripted <italic>s</italic>
variables are the standard deviations of the subscripted quantities and (<italic>s</italic>
<sub><italic>CL</italic>
</sub>
/<italic>CL</italic>
)<sup>2</sup>
is the squared coefficient of variation CV<sup>2</sup>
of <italic>CL.</italic>
Minimizing the right-hand side of Eq. <xref rid="Equ4" ref-type="">4</xref>
as a function of the shrinkage, <italic>λ</italic>
, selects a <italic>λ</italic>
value that produces the <italic>CL</italic>
value with the smallest relative error achievable. Also, minimizing the relative error in <italic>CL</italic>
is indispensable for making reliable measures of <italic>CL</italic>
when <italic>CL</italic>
is small.</p>
<p>The squared coefficient of variation, CV<sup>2</sup>
, of the individual volumes of distribution, i.e. (<italic>s</italic>
<sub><italic>V</italic>
</sub>
/<italic>V</italic>
)<sup>2</sup>
, is calculated from the SDs produced from minimizing Eq. <xref rid="Equ4" ref-type="">4</xref>
substituted into<disp-formula id="Equ5"><label>5</label>
<alternatives><tex-math id="M28">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\left( {\frac{{{s_V}}}{V}} \right)^2} = s_{\alpha }^2{\left( {\frac{1}{\alpha } + {{\bar{a}}_1} + \ln \beta - \Psi \left( \alpha \right)} \right)^2} + s_{\beta }^2{\left( {\frac{\alpha }{\beta } - \frac{1}{\beta } - {{\bar{a}}_2}} \right)^2} + 2{s_{{\alpha \beta }}}\left( {\frac{1}{\alpha } + {{\bar{a}}_1} + \ln \beta - \Psi \left( \alpha \right)} \right)\left( {\frac{\alpha }{\beta } - \frac{1}{\beta } - {{\bar{a}}_2}} \right)... $$\end{document}</tex-math>
<graphic xlink:href="259_2011_1887_Article_Equ5.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
Equation <xref rid="Equ5" ref-type="">5</xref>
is from application of the propagation of error formula to the substitution of Eq. <xref rid="Equ3" ref-type="">3</xref>
into <italic>V</italic>
=(<italic>α/β</italic>
)<italic>CL</italic>
.</p>
</sec>
</app>
</app-group>
<glossary><title>Abbreviations</title>
<def-list><def-list><def-item><term>α</term>
<def><p>Shape parameter of GV, volume scale factor, and plasma leak term Eq. <xref rid="Equ1" ref-type="">1</xref>
</p>
</def>
</def-item>
<def-item><term><italic>AUC</italic>
</term>
<def><p>Area under the curve of concentration in time, i.e. <inline-formula id="IEq1"><alternatives><tex-math id="M29">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \int_0^{\infty } {{C_{{obs}}}(t)dt} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq1.gif"></inline-graphic>
</alternatives>
</inline-formula>
</p>
</def>
</def-item>
<def-item><term><italic>β</italic>
</term>
<def><p>Rate constant for GV (per min)</p>
</def>
</def-item>
<def-item><term><italic>CL</italic>
</term>
<def><p>Plasma clearance, <italic>D</italic>
/<italic>AUC</italic>
</p>
</def>
</def-item>
<def-item><term><italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
)</term>
<def><p>Observed concentration, in Bq/ml, often scaled as dose percent per ml</p>
</def>
</def-item>
<def-item><term><italic>C</italic>
(<italic>t</italic>
)</term>
<def><p>Continuous function for estimation of <italic>C</italic>
<sub><italic>obs</italic>
</sub>
(<italic>t</italic>
)</p>
</def>
</def-item>
<def-item><term>CV</term>
<def><p>Coefficient of variation, SD/mean</p>
</def>
</def-item>
<def-item><term><italic>D</italic>
</term>
<def><p>The dosage in Bq or per cent</p>
</def>
</def-item>
<def-item><term>E1, E2, E<italic>n</italic>
</term>
<def><p>Sum of exponential terms model with 1, 2 or <italic>n</italic>
terms, also, 1, 2, or <italic>n</italic>
-compartment model</p>
</def>
</def-item>
<def-item><term>GV</term>
<def><p>Gamma variate, <inline-formula id="IEq2"><alternatives><tex-math id="M30">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ C(t) = K{t^{{\alpha - 1}}}{e^{{ - \beta t}}} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq2.gif"></inline-graphic>
</alternatives>
</inline-formula>
</p>
</def>
</def-item>
<def-item><term>Γ</term>
<def><p>Tikhonov matrix [do not confuse with Γ(<italic>x</italic>
)]</p>
</def>
</def-item>
<def-item><term>Γ(<italic>x</italic>
)</term>
<def><p>The gamma function, <inline-formula id="IEq3"><alternatives><tex-math id="M31">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \Gamma (x) = \int_0^{\infty } {{t^{{x - 1}}}{e^{{ - t}}}dt} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq3.gif"></inline-graphic>
</alternatives>
</inline-formula>
, is a generalization of the factorial, <inline-formula id="IEq4"><alternatives><tex-math id="M32">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \Gamma (n) = \left( {n - 1} \right)! $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq4.gif"></inline-graphic>
</alternatives>
</inline-formula>
</p>
</def>
</def-item>
<def-item><term><italic>K</italic>
</term>
<def><p>Constant of proportionality for GV</p>
</def>
</def-item>
<def-item><term><italic>MRT</italic>
</term>
<def><p>Mean residence time</p>
</def>
</def-item>
<def-item><term>Ψ(<italic>α</italic>
)</term>
<def><p>(Psi) digamma function of <italic>α</italic>
is <inline-formula id="IEq5"><alternatives><tex-math id="M33">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \Psi \left( \alpha \right) = d\left[ {\ln \Gamma \left( \alpha \right)} \right]/d\alpha = \Gamma \prime \left( \alpha \right)/\Gamma \left( \alpha \right) $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq5.gif"></inline-graphic>
</alternatives>
</inline-formula>
</p>
</def>
</def-item>
<def-item><term><inline-formula id="IEq6"><alternatives><tex-math id="M34">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {s_{\alpha }},{s_{\beta }},{s_{{CL}}} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq6.gif"></inline-graphic>
</alternatives>
</inline-formula>
</term>
<def><p>SDs of <inline-formula id="IEq7"><alternatives><tex-math id="M35">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \alpha, \beta, CL $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq7.gif"></inline-graphic>
</alternatives>
</inline-formula>
for Tk-GV Eq. <xref rid="Equ4" ref-type="">4</xref>
</p>
</def>
</def-item>
<def-item><term><inline-formula id="IEq8"><alternatives><tex-math id="M36">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {s_V} $$\end{document}</tex-math>
<inline-graphic xlink:href="259_2011_1887_Article_IEq8.gif"></inline-graphic>
</alternatives>
</inline-formula>
</term>
<def><p>SD of <italic>V</italic>
from Tk-GV Eq. <xref rid="Equ5" ref-type="">5</xref>
</p>
</def>
</def-item>
<def-item><term>SD</term>
<def><p>Standard deviation</p>
</def>
</def-item>
<def-item><term>Tk-GV</term>
<def><p>Tikhonov regularization of GV fit of CV of <italic>CL</italic>
</p>
</def>
</def-item>
<def-item><term><italic>V</italic>
</term>
<def><p>Volume of distribution of marker (ml or l)</p>
</def>
</def-item>
</def-list>
</def-list>
</glossary>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Sante/explor/LymphedemaV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000D260 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000D260 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Wicri/Sante |area= LymphedemaV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.31. | ![]() |