Abstract:Abstract Objective . To present an efficient uncertainty quantification method for range and set-up errors in Monte Carlo (MC) dose calculations. Further, we show that uncertainty induced by interplay and other dynamic influences may be approximated using suitable error correlation models. Approach . We introduce an importance (re-)weighting method in MC history scoring to concurrently construct estimates for error scenarios, the expected dose and its variance from a single set of MC simulated particle histories. The approach relies on a multivariate Gaussian input and uncertainty model, which assigns probabilities to the initial phase space sample, enabling the use of different correlation models. Through modification of the phase space parameterization, accuracy can be traded between that of the uncertainty or the nominal dose estimate. Main results . The method was implemented using the MC code TOPAS and validated for proton intensity-modulated particle therapy (IMPT) with reference scenario estimates. We achieve accurate results for set-up uncertainties ( γ 2 mm/2% ≥ 99.01% ( E [ d ]), γ 2 mm/2% ≥ 98.04% ( σ ( d ))) and expectedly lower but still sufficient agreement for range uncertainties, which are approximated with uncertainty over the energy distribution. Here pass rates of 99.39% ( E [ d ])/ 93.70% ( σ ( d )) (range errors) and 99.86% ( E [ d ])/ 96.64% ( σ ( d )) (range and set-up errors) can be achieved. Initial evaluations on a water phantom, a prostate and a liver case from the public CORT dataset show that the CPU time decreases by more than an order of magnitude. Significance . The high precision and conformity of IMPT comes at the cost of susceptibility to treatment uncertainties in particle range and patient set-up. Yet, dose uncertainty quantification and mitigation, which is usually based on sampled error scenarios, becomes challenging when computing the dose with computationally expensive but accurate MC simulations. As the results indicate, the proposed method could reduce computational effort while also facilitating the use of high-dimensional uncertainty models.

Importance sampling method for Feldman-Cousins confidence intervals

An importance sampling method for Feldman-Cousins confidence intervals

Importance Sampling of Interval Markov Chains

The Profiled Feldman-Cousins technique for confidence interval construction in the presence of nuisance parameters

A Unified Approach to the Classical Statistical Analysis of Small Signals

Undercoverage in high-statistics counting experiments with finite MC samples

Nonasymptotic Bounds for Suboptimal Importance Sampling

Estimating Ratios of Normalizing Constants Using Linked Importance Sampling

Statistical significance testing for mixed priors: a combined Bayesian and frequentist analysis

Uncertainty Quantification and Confidence Intervals for Naive Rare-Event Estimators

Efficient Importance Sampling for Rare Event Simulation with Applications

Efficient uncertainty quantification for Monte Carlo dose calculations using importance (re-)weighting

Importance Sampling for the Empirical Measure of Weakly Interacting Diffusions

Guaranteed Conservative Fixed Width Confidence Intervals Via Monte Carlo Sampling

Predicting the statistical error of analog particle tracing Monte Carlo

Efficient and Unbiased Sampling of Boltzmann Distributions via Consistency Models

Iterative Importance Sampling Algorithms for Parameter Estimation

On the importance of factorization for fast binned likelihood inference

On importance sampling and independent Metropolis-Hastings with an unbounded weight function

A probabilistic simulation method for sensitivity analysis of input epistemic uncertainties on failure probability

Repro Samples Method for Finite- and Large-Sample Inferences