Using the Stochastic Collocation Method for the Uncertainty Quantification of Drug Concentration Due to Depot Shape Variability

Numerical simulations entail modeling assumptions that impact outcomes. Therefore, characterizing, in a probabilistic sense, the relationship between the variability of model selection and the variability of outcomes is important. Under certain assumptions, the stochastic collocation method offers a computationally feasible alternative to traditional Monte Carlo approaches for assessing the impact of model and parameter variability. We propose a framework that combines component shape parameterization with the stochastic collocation method to study the effect of drug depot shape variability on the outcome of drug diffusion simulations in a porcine model. We use realistic geometries segmented from MR images and employ level-set techniques to create two alternative univariate shape parameterizations. We demonstrate that once the underlying stochastic process is characterized, quantification of the introduced variability is quite straightforward and provides an important step in the validation and verification process.

[1]  N. Zabaras,et al.  Stochastic inverse heat conduction using a spectral approach , 2004 .

[2]  G. W. Snedecor Statistical Methods , 1964 .

[3]  Dongbin Xiu,et al.  High-Order Collocation Methods for Differential Equations with Random Inputs , 2005, SIAM J. Sci. Comput..

[4]  Guido Gerig,et al.  User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability , 2006, NeuroImage.

[5]  David E. Breen,et al.  A level-set approach for the metamorphosis of solid models , 1999, SIGGRAPH '99.

[6]  J. Escher,et al.  The volume preserving mean curvature flow near spheres , 1998 .

[7]  D. Xiu,et al.  Modeling uncertainty in flow simulations via generalized polynomial chaos , 2003 .

[8]  D. Xiu,et al.  Modeling Uncertainty in Steady State Diffusion Problems via Generalized Polynomial Chaos , 2002 .

[9]  Hoon Kim,et al.  Monte Carlo Statistical Methods , 2000, Technometrics.

[10]  Ross T. Whitaker,et al.  A Level-Set Approach to 3D Reconstruction from Range Data , 1998, International Journal of Computer Vision.

[11]  Kenneth I. Joy,et al.  Shell maps , 2005, ACM Trans. Graph..

[12]  Alejandro F. Frangi,et al.  3D MRA coronary axis determination using a minimum cost path approach , 2002, Magnetic resonance in medicine.

[13]  Baskar Ganapathysubramanian,et al.  Modeling diffusion in random heterogeneous media: Data-driven models, stochastic collocation and the variational multiscale method , 2007, J. Comput. Phys..

[14]  L. Trefethen Spectral Methods in MATLAB , 2000 .

[15]  Klaus Gärtner,et al.  Meshing Piecewise Linear Complexes by Constrained Delaunay Tetrahedralizations , 2005, IMR.

[16]  D. Xiu Efficient collocational approach for parametric uncertainty analysis , 2007 .

[17]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[18]  Robert Michael Kirby,et al.  Application of Stochastic Finite Element Methods to Study the Sensitivity of ECG Forward Modeling to Organ Conductivity , 2008, IEEE Transactions on Biomedical Engineering.

[19]  Steve Oudot,et al.  Provably good sampling and meshing of surfaces , 2005, Graph. Model..

[20]  D. Blumenthal,et al.  Efficacy of local dipyridamole therapy in a porcine model of arteriovenous graft stenosis. , 2006, Kidney international.

[21]  Marko Subasic,et al.  Level Set Methods and Fast Marching Methods , 2003 .