University of Surrey

Test tubes in the lab Research in the ATI Dance Research

A comparison of Bayesian and non-linear regression methods for robust estimation of pharmacokinetics in DCE-MRI and how it affects cancer diagnosis

Dikaios, N, Atkinson, D, Tudisca, C, Purpura, P, Forster, M, Ahmed, H, Beale, T, Emberton, E and Punwani, S (2017) A comparison of Bayesian and non-linear regression methods for robust estimation of pharmacokinetics in DCE-MRI and how it affects cancer diagnosis Computerized Medical Imaging and Graphics, 56. pp. 1-10.

cmig.pdf - Accepted version Manuscript

Download (1MB) | Preview


The aim of this work is to compare Bayesian Inference for nonlinear models with commonly used traditional non-linear regression (NR) algorithms for estimating tracer kinetics in Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI). The algorithms are compared in terms of accuracy, and reproducibility under different initialization settings. Further it is investigated how a more robust estimation of tracer kinetics affects cancer diagnosis. The derived tracer kinetics from the Bayesian algorithm were validated against traditional NR algorithms (i.e. Levenberg-Marquardt, simplex) in terms of accuracy on a digital DCE phantom and in terms of goodness-of-fit (Kolmogorov-Smirnov test) on ROI-based concentration time courses from two different patient cohorts. The first cohort consisted of 76 men, 20 of whom had significant peripheral zone prostate cancer (any cancer-core-length (CCL) with Gleason > 3 + 3 or any-grade with CCL > = 4 mm) following transperineal template prostate mapping biopsy. The second cohort consisted of 9 healthy volunteers and 24 patients with head and neck squamous cell carcinoma. The diagnostic ability of the derived tracer kinetics was assessed with receiver operating characteristic area under curve (ROC AUC) analysis. The Bayesian algorithm accurately recovered the ground-truth tracer kinetics for the digital DCE phantom consistently improving the Structural Similarity Index (SSIM) across the 50 different initializations compared to NR. For optimized initialization, Bayesian did not improve significantly the fitting accuracy on both patient cohorts, and it only significantly improved the ve ROC AUC on the HN population from ROC AUC = 0.56 for the simplex to ROC AUC = 0.76. For both cohorts, the values and the diagnostic ability of tracer kinetic parameters estimated with the Bayesian algorithm weren’t affected by their initialization. To conclude, the Bayesian algorithm led to a more accurate and reproducible quantification of tracer kinetic parameters in DCE-MRI, improving their ROC-AUC and decreasing their dependence on initialization settings.

Item Type: Article
Subjects : Computing
Divisions : Faculty of Engineering and Physical Sciences > Electronic Engineering
Authors :
Atkinson, D
Tudisca, C
Purpura, P
Forster, M
Ahmed, H
Beale, T
Emberton, E
Punwani, S
Date : 5 February 2017
DOI : 10.1016/j.compmedimag.2017.01.003
Copyright Disclaimer : Crown Copyright © 2017 Published by Elsevier Ltd. All rights reserved. This manuscript version is made available under the CC-BY-NC-ND 4.0 license
Uncontrolled Keywords : DCE analysis; Bayesian inference for nonlinear model; Prostate cancer; Head and neck
Depositing User : Symplectic Elements
Date Deposited : 17 May 2017 13:58
Last Modified : 16 Jan 2019 18:52

Actions (login required)

View Item View Item


Downloads per month over past year

Information about this web site

© The University of Surrey, Guildford, Surrey, GU2 7XH, United Kingdom.
+44 (0)1483 300800