A stochastic frontier analysis for enhanced treatment quality of high-dose-rate brachytherapy plans

The purpose of the present study is to develop patient specific unbiased quality control (QC) models for high dose rate (HDR) brachytherapy plans. The proposed models are based on the stochastic frontier analysis formalism, a method of economic modeling. They act as a QC tool by predicting before the treatment planning process starts, the dosimetric coverage achievable for a HDR brachytherapy prostate plan. The geometric parameters considered in developing the models were: patient clinical target volume (CTV), organs at risk (OAR) volume, the bidirectional Hausdorff distance between CTV and OARs, and a fourth parameter measuring the catheters degree of non-parallelism within the target volume. Dosimetry parameters of interest are V100 for the CTV, V75 (bladder, rectum) and D10 (urethra). Results show that the built models can provide valuable information on the personalization of the optimization process based on the patient geometric parameters. The impact on the quality plan due to the planner's experience variability and judgment can be reduced by using those models, since the planner will attempt to achieve dosimetric parameters predicted by the models. Furthermore, the models provide information on the better trade-off between the target volume coverage and OARs sparing that can be achieved, regardless of the planner's experience; the latter being achieved by moving each plan at least around their respective frontier for V100, V75 and D10. The shortfall of the dosimetric parameters values computed by the treatment planning system (TPS) from those predicted by the models for a proportion of plans in the dataset reveals that optimized plans from a TPS, even clinically acceptable, are not necessarily the best that could be achieved. These represent 83% of plans in the training set for the target volume coverage (V100), ∼50% for the bladder (V75) and  ∼72% for the urethra (D10).

[1]  Alper Atamtürk,et al.  IPIP: A new approach to inverse planning for HDR brachytherapy by directly optimizing dosimetric indices. , 2010, Medical physics.

[2]  I. Rosenberg,et al.  Radiation Oncology Physics: A Handbook for Teachers and Students , 2008, British Journal of Cancer.

[3]  Yang Xiang,et al.  Generalized Simulated Annealing for Global Optimization: The GenSA Package , 2013, R J..

[4]  Jimm Grimm,et al.  Anatomy-based inverse planning simulated annealing optimization in high-dose-rate prostate brachytherapy: significant dosimetric advantage over other optimization techniques. , 2008, International journal of radiation oncology, biology, physics.

[5]  W. Greene A Gamma-distributed stochastic frontier model , 1990 .

[6]  Chun-qiang Zhu,et al.  Application of Hausdorff distance in image matching , 2014, 2014 IEEE Workshop on Electronics, Computer and Applications.

[7]  Etienne Lessard Development and clinical introduction of an inverse planning dose optimization by simulated annealing (IPSA) for high dose rate brachytherapy , 2004 .

[8]  D. Aigner,et al.  P. Schmidt, 1977,?Formulation and estimation of stochastic frontier production function models,? , 1977 .

[9]  Jean Pouliot,et al.  Inverse planning for interstitial gynecologic template brachytherapy: truly anatomy-based planning. , 2002, International journal of radiation oncology, biology, physics.

[10]  E. Pissaloux,et al.  A guided image matching approach using Hausdorff distance with interesting points detection , 1994, Proceedings of 1st International Conference on Image Processing.

[11]  R. Stevenson Likelihood functions for generalized stochastic frontier estimation , 1980 .

[12]  C. Lovell,et al.  Estimating technical and allocative inefficiency relative to stochastic production and cost frontiers , 1979 .

[13]  Daniel P. Huttenlocher,et al.  Comparing Images Using the Hausdorff Distance , 1993, IEEE Trans. Pattern Anal. Mach. Intell..

[14]  Milan Sonka,et al.  3D Slicer as an image computing platform for the Quantitative Imaging Network. , 2012, Magnetic resonance imaging.

[15]  W. Meeusen,et al.  Efficiency Estimation from Cobb-Douglas Production Functions with Composed Error , 1977 .

[16]  J Pouliot,et al.  Inverse planning anatomy-based dose optimization for HDR-brachytherapy of the prostate using fast simulated annealing algorithm and dedicated objective function. , 2001, Medical physics.

[17]  Jean Pouliot,et al.  Comparison of inverse planning simulated annealing and geometrical optimization for prostate high-dose-rate brachytherapy. , 2004, Brachytherapy.

[18]  P. Halina,et al.  A comparison of anatomy-based inverse planning with simulated annealing and graphical optimization for high-dose-rate prostate brachytherapy. , 2008, Brachytherapy.

[19]  S. S. Wilks The Large-Sample Distribution of the Likelihood Ratio for Testing Composite Hypotheses , 1938 .

[20]  C. Lovell,et al.  On the estimation of technical inefficiency in the stochastic frontier production function model , 1982 .

[21]  F. Khan The physics of radiation therapy , 1985 .