Non-linear mixed-effects modelling and population-based model selection for 131I kinetics in benign thyroid disease
Year: 2025
Type: article
Abstract: Purpose This study aimed to determine a mathematical model for accurately calculating time-integrated activities (TIAs) of target tissue in 131I therapy for benign thyroid disease using the population-based model selection and non-linear mixed-effects (PBMS-NLME) method. Methods Biokinetic data of 131I in target tissue were collected from seventy-three patients at 2, 6, 24, 48, and 96 (N = 53) or 120 (N = 20) h after oral capsule administration with 1 MBq 131I. Based on the Akaike weight, the best sum-of-exponential function (SOEF) describing the biokinetic data was selected using PBMS-NLME modelling. Nine SOEF with three to six parameters (including the function from the European Association of Nuclear Medicine Standard Operational Procedure (EANM SOP)) were used. The fittings were repeated 1000 times with different starting values of the SOE parameters to find the optimal fit. Akaike weight was used to identify the performance of the best model from PBMS-NLME and the EANM SOP SOE function with individual fitting. Results Based on the PBMS-NLME analysis, the SOEF $$\frac{{\lambda_{1} }}{{\lambda_{2} + \lambda_{1} - \lambda_{3} }}\left( {e^{{ - \left( {\lambda_{3} + \lambda_{phys} } \right)t}} - e^{{ - \left( {\lambda_{1} + \lambda_{2} + \lambda_{phys} } \right)t}} } \right) + a_{1} e^{{ - \left( {\lambda_{1} + \lambda_{2} + \lambda_{phys} } \right)t}}$$ λ 1 λ 2 + λ 1 - λ 3 e - λ 3 + λ phys t - e - λ 1 + λ 2 + λ phys t + a 1 e - λ 1 + λ 2 + λ phys t was selected as the function most supported by the data. The Akaike weight of the best function was approximately 100%. The best SOEF from the PBMS-NLME approach shows a better performance in describing the biokinetic data of 131I in the thyroid gland than the function from the EANM SOP with individual fitting, based on the Akaike weight. Conclusions The best mathematical model from the PBMS-NLME approach has one more free parameter than the EANM SOP function, which could lead to more accurate TIAs.
Source: EJNMMI Physics
Language: English
Subfield: Radiology, Nuclear Medicine and Imaging
Field: Medicine
Domain: Health Sciences
Open Access status: gold
Estimated APC paid: $2,490
OpenAlex ID: W4409277976