Optimal Confidence for Monte Carlo Integration of Smooth Functions

Duration: 36 mins
Share this media item:
Embed this media item:


About this item
Image inherited from collection
Description: Kunsch, R
Thursday 21st February 2019 - 13:40 to 14:15
 
Created: 2019-02-21 16:18
Collection: Approximation, sampling and compression in data science
Publisher: Isaac Newton Institute
Copyright: Kunsch, R
Language: eng (English)
Distribution: World     (downloadable)
Explicit content: No
Aspect Ratio: 16:9
Screencast: No
Bumper: UCS Default
Trailer: UCS Default
 
Abstract: 1−δ∈(0,1) using function evaluations as information within randomized algorithms. Methods that achieve optimal rates in terms of the root mean square error (RMSE) are not always optimal in terms of error at confidence, usually we need some non-linearity in order to suppress outliers. Besides, there are numerical problems which can be solved in terms of error at confidence but no algorithm can guarantee a finite RMSE, see [1]. Hence, the new error criterion seems to be more general than the classical RMSE. The sharp order for multivariate functions from classical isotropic Sobolev spaces Wrp([0,1]d) can be achieved via control variates, as long as the space is embedded in the space of continuous functions C([0,1]d). It turns out that the integrability index p has an effect on the influence of the uncertainty δ to the complexity, with the limiting case p=1 where deterministic methods cannot be improved by randomization. In general, higher smoothness reduces the effort we need to take in order to increase the confidence level. Determining the complexity n(ε,δ) is much more challenging for mixed smoothness spaces Wrp([0,1]d). While optimal rates are known for the classical RMSE (as long as Wrp([0,1]d) is embedded in L2([0,1]d)), see [2], basic modifications of the corresponding algorithms fail to match the theoretical lower bounds for approximating the integral with prescribed confidence.

Joint work with Daniel Rudolf

[1] R.J. Kunsch, E. Novak, D. Rudolf. Solvable integration problems and optimal sample size selection. To appear in Journal of Complexity.
[2] M. Ullrich. A Monte Carlo method for integration of multivariate smooth functions. SIAM Journal on Numerical Analysis, 55(3):1188-1200, 2017.

Available Formats
Format Quality Bitrate Size
MPEG-4 Video 640x360    1.94 Mbits/sec 524.14 MB View Download
WebM 640x360    530.95 kbits/sec 140.06 MB View Download
iPod Video 480x270    522.08 kbits/sec 137.66 MB View Download
MP3 44100 Hz 249.8 kbits/sec 65.93 MB Listen Download
Auto * (Allows browser to choose a format it supports)