Single-time-point dosimetry using model selection and nonlinear mixed-effects modelling: a proof of concept
Deni Hardiansyah, Ade Riana, Ambros J. Beer, Gerhard Glatting
Abstract
Abstract Purpose This project aims to develop and evaluate a method for accurately determining time-integrated activities (TIAs) in single-time-point (STP) dosimetry for molecular radiotherapy. It performs a model selection (MS) within the framework of the nonlinear mixed-effects (NLME) model (MS–NLME). Methods Biokinetic data of [ 111 In]In-DOTATATE activity in kidneys at T 1 = (2.9 ± 0.6) h, T 2 = (4.6 ± 0.4) h, T 3 = (22.8 ± 1.6) h, T 4 = (46.7 ± 1.7) h, and T 5 = (70.9 ± 1.0) h post injection were obtained from eight patients using planar imaging. Eleven functions were derived from various parameterisations of mono-, bi-, and tri-exponential functions. The functions’ fixed and random effects parameters were fitted simultaneously (in the NLME framework) to the biokinetic data of all patients. The Akaike weights were used to select the fit function most supported by the data. The relative deviations (RD) and the root-mean-square error (RMSE) of the calculated TIAs for the STP dosimetry at T 3 = (22.8 ± 1.6) h and T 4 = (46.7 ± 1.7) h p.i. were determined for all functions passing the goodness-of-fit test. Results The function $$f_{4d} \left( t \right) = A_{1} /\left\{ {\left( {\frac{1 - \alpha }{{\lambda_{1} + \lambda_{{{\text{phys}}}} }}} \right) - \left( {\frac{\alpha }{{\lambda_{2} + \lambda_{{{\text{phys}}}} }}} \right) - \left( {\frac{1 - 2\alpha }{{\lambda_{bc} + \lambda_{{{\text{phys}}}} }}} \right)} \right\} \cdot e^{{ - \lambda_{{{\text{phys}}}} t}} \cdot \left\{ {\left( {1 - \alpha } \right) \cdot e^{{ - \lambda_{1} t}} - \alpha \cdot e^{{ - \lambda_{2} t}} - \left( {1 - 2\alpha } \right) \cdot e^{{ - \lambda_{bc} t}} } \right\}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mi>f</mml:mi> <mml:mrow> <mml:mn>4</mml:mn> <mml:mi>d</mml:mi> </mml:mrow> </mml:msub> <mml:mfenced> <mml:mi>t</mml:mi> </mml:mfenced> <mml:mo>=</mml:mo> <mml:msub> <mml:mi>A</mml:mi> <mml:mn>1</mml:mn> </mml:msub> <mml:mo>/</mml:mo> <mml:mfenced> <mml:mrow> <mml:mfenced> <mml:mfrac> <mml:mrow> <mml:mn>1</mml:mn> <mml:mo>-</mml:mo> <mml:mi>α</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mi>λ</mml:mi> <mml:mn>1</mml:mn> </mml:msub> <mml:mo>+</mml:mo> <mml:msub> <mml:mi>λ</mml:mi> <mml:mtext>phys</mml:mtext> </mml:msub> </mml:mrow> </mml:mfrac> </mml:mfenced> <mml:mo>-</mml:mo> <mml:mfenced> <mml:mfrac> <mml:mi>α</mml:mi> <mml:mrow> <mml:msub> <mml:mi>λ</mml:mi> <mml:mn>2</mml:mn> </mml:msub> <mml:mo>+</mml:mo> <mml:msub> <mml:mi>λ</mml:mi> <mml:mtext>phys</mml:mtext> </mml:msub> </mml:mrow> </mml:mfrac> </mml:mfenced> <mml:mo>-</mml:mo> <mml:mfenced> <mml:mfrac> <mml:mrow> <mml:mn>1</mml:mn> <mml:mo>-</mml:mo> <mml:mn>2</mml:mn> <mml:mi>α</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mi>λ</mml:mi> <mml:mrow> <mml:mi>bc</mml:mi> </mml:mrow> </mml:msub> <mml:mo>+</mml:mo> <mml:msub> <mml:mi>λ</mml:mi> <mml:mtext>phys</mml:mtext> </mml:msub> </mml:mrow> </mml:mfrac> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo>·</mml:mo> <mml:msup> <mml:mi>e</mml:mi> <mml:mrow> <mml:mo>-</mml:mo> <mml:msub> <mml:mi>λ</mml:mi> <mml:mtext>phys</mml:mtext> </mml:msub> <mml:mi>t</mml:mi> </mml:mrow> </mml:msup> <mml:mo>·</mml:mo> <mml:mfenced> <mml:mrow> <mml:mfenced> <mml:mrow> <mml:mn>1</mml:mn> <mml:mo>-</mml:mo> <mml:mi>α</mml:mi> </mml:mrow> </mml:mfenced> <mml:mo>·</mml:mo> <mml:msup> <mml:mi>e</mml:mi> <mml:mrow> <mml:mo>-</mml:mo> <mml:msub> <mml:mi>λ</mml:mi> <mml:mn>1</mml:mn> </mml:msub> <mml:mi>t</mml:mi> </mml:mrow> </mml:msup> <mml:mo>-</mml:mo> <mml:mi>α</mml:mi> <mml:mo>·</mml:mo> <mml:msup> <mml:mi>e</mml:mi> <mml:mrow> <mml:mo>-</mml:mo> <mml:msub> <mml:mi>λ</mml:mi> <mml:mn>2</mml:mn> </mml:msub> <mml:mi>t</mml:mi> </mml:mrow> </mml:msup> <mml:mo>-</mml:mo> <mml:mfenced> <mml:mrow> <mml:mn>1</mml:mn> <mml:mo>-</mml:mo> <mml:mn>2</mml:mn> <mml:mi>α</mml:mi> </mml:mrow> </mml:mfenced> <mml:mo>·</mml:mo> <mml:msup> <mml:mi>e</mml:mi> <mml:mrow> <mml:mo>-</mml:mo> <mml:msub> <mml:mi>λ</mml:mi> <mml:mrow> <mml:mi>bc</mml:mi> </mml:mrow> </mml:msub> <mml:mi>t</mml:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> with four adjustable parameters and $$\lambda_{bc} = \frac{{{\text{ln}}\left( 2 \right)}}{{1\;{\text{ min}}}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mi>λ</mml:mi> <mml:mrow> <mml:mi>bc</mml:mi> </mml:mrow> </mml:msub> <mml:mo>=</mml:mo> <mml:mfrac> <mml:mrow> <mml:mtext>ln</mml:mtext> <mml:mfenced> <mml:mn>2</mml:mn> </mml:mfenced>