Quasi-Monte Carlo Methods for Applications in Statistics
|
|
- Doreen Wells
- 5 years ago
- Views:
Transcription
1 Quasi-Monte Carlo Methods for Applications in Statistics Weights for QMC in Statistics Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
2 Quasi-Monte Carlo Methods for Applications in Statistics Joint work with Ian Sloan and Frances Kuo Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
3 Quasi-Monte Carlo Methods for Applications in Statistics Outline: Background on QMC A class of problems in Statistics Approach and preliminary results Future research Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
4 Integrals over Euclidean space Many problems involve evaluation of integrals given by I d (f, ρ) = f (x)ρ(x) dx, R d where ρ(x) is a probability density. Hence ρ(x) 0 for any x R d and R d ρ(x) dx = 1. Here we assume that the density is expressed as the product of univariate densities, hence d ρ(x) = φ(x j ). j=1 Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
5 Transforming the integral Integrals over unbounded regions are usually mapped to the unit cube. Then, a QMC method (for instance a shifted lattice rule) can be used to approximate the integral. In the 1-dimensional case, we can use the following transform: u = Φ(x) = x φ(t) dt, x R. The inverse mapping will be Φ 1 : (0, 1) R, Φ 1 (u) = x. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
6 Mapping to the unit cube In the d-dimensional case, if x = (x 1, x 2,..., x d ) R d, then Φ(x) = (Φ(x 1 ), Φ(x 2 ),..., Φ(x d )). In the same manner, the inverse mapping will also be applied component-wise. The integral becomes I d (f, ρ) = f (Φ 1 (u)) du = [0,1] d g(u) du = I d (g), [0,1] d where g = f Φ 1 (applied component-wise). These integrals can be approximated by quadrature rules of the form Q n,d (g) = 1 n 1 g(w k ) = 1 n 1 f (t k ), n n k=0 k=0 where t k = Φ 1 (w k ) R d, for all 0 k n 1. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
7 QMC error The typical error bound is expressed as I d (g) Q n,d (g) e worst n,d,γ (P) g (0,1) d, where P is the set of quadrature points. Note that we have g (0,1) d = f R d. Provided that g <, the vast majority of research papers was focused on obtaining the optimal convergence rate for the worst-case error. It has been proved that good lattice rules do exist and that there are efficient algorithms (such as the CBC construction) to construct such lattice rules under various settings. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
8 Weights A widely used setting is provided by weighted reproducing kernel Hilbert spaces (RKHS) and L 2 -type discrepancies. Typical examples are weighted Sobolev and weighted Korobov spaces of functions. The concepts of weights is related to the fact that variables or group of variables have different importance. Accordingly, we associate a weight γ u, for each subset of coordinates u D := {1, 2,..., d}. It is known that under suitable conditions on the weights, the weighted discrepancy or the worst-case error is independent of the dimension. A typical example of convergence rate on the worst-case error is O(n 1+δ ) with the involved constant independent of the dimension. It is also known how to construct lattice rules that achieve such a convergence rate. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
9 Our approach Instead of analysing only the worst-case error, we look at the error bound itself, namely we analyse the product between the worst-case error and the norm of the integrand. This idea is not new and there are however some earlier works by Dick, Sloan, Wang and Woźniakowski (2004) and Larcher, Leobacher and Scheicher (2003). More recent contributions coming from the analysis of the error bound is due to Kuo, Schwab and Sloan who obtained the optimal weights for a class of QMC methods suitable for PDE. We follow a similar idea, but for a class of important problems in statistics. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
10 A class of problems in Statistics We are interested in a class of highly structured generalised response models in Statistics, in particular on the so-called generalised linear mixed model (GLMM). The Poisson-type likelihood of such a model can be expressed by L(β, Σ) = R d d exp(y j (w j + β) e w j +β exp ( 1 2 wt Σ 1 w ) dw. y j (2π) d det(σ) j=1 However a naïve approach consisting of diagonalising the covariance matrix and then mapping to the unit cube gives poor results!! See Kuo, Dunsmuir, Sloan, Wand and Womersley (2008) for details. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
11 Recentering and rescaling Good results however have been obtained after an appropriate linear transformation consisting of recentering and rescaling. Denote F (w) = d j=1 ( y j (w j + β) e w j +β ) 1 2 wt Σ 1 w. We use a change of variables of the form w A x + w, where w satisfies F (w ) = 0 and A A T = ( 2 F (w )) 1. That is, w is the stationary point, while the matrix A A T describes the convexity of F around the stationary point. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
12 The transformed integral The integral we want to evaluate becomes d exp F (w) dw = det(a ) F (A x + w ) R d R d where φ(x) is a univariate probability density. j=1 1 φ(x j ) } {{ } f (x) d φ(x j ) dx, This integral can then be mapped to the unit cube as explained earlier. j=1 Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
13 Our function We want to evaluate the norm of functions f given by d ( ) f (x) := exp y j ((A x) j + β + wj ) e (A x) j +β+wj exp j=1 [ 1 ] 2 (A x + w ) T Σ 1 (A x + w ) d h(x j ), where h(x) := 1/φ(x). The norm of such functions is defined by f := ( ) 2 γ 1 u 1/2 u f ((x u, 0)) dx u. x u u D R u j=1 Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
14 How to find a bound on the norm of f? The task of finding bounds is also difficult. However we follow three main steps: 1 We find the partial mixed first derivatives of f. 2 We bound the derivatives. 3 We bound the norm. Following these steps, we obtained a generic bound on the norm f and a specific bound in the case of an assumed logistic density. Similar ideas can be used to deduce corresponding bounds for other densities (future research). Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
15 The partial derivatives of f After some tedious calculations, by applying Faà di Bruno formula and the formula for differentiating products of functions, we obtain: u f (x) = exp(g(x)) d e (A x+w ) j +β aji + T w (x) x u v u π Π(v) w π j=1 i w d d h(x i ) h (x i ), where i=1 i u v i=1 i u v Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
16 The partial derivatives of f g(x) := d j=1 ( ) y j ((A x) j + β + wj ) e (A x) j +β+wj and 1 2 (A x + w ) T Σ 1 (A x + w ), d j=1 y ja ji (A ) i Σ 1 (A x + w ), if T w (x) := (A ) T i Σ 1 (A ) k, if w = {i, k}, 0, if w 3. w = {i}, Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
17 Expressing the bounds on the derivatives Theorem The partial derivatives of f can be bounded by ( u f (x) d (A ) T i τ exp J(x) + Σ 1 A x ) x u α i i=1 d B v v!ω v α i h(x i ) v u 2 i v i=1 i u v d i=1 i u v h (x i ), where B m denotes the Bell number of order m, that is the number of partitions of a set with cardinality m. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
18 Expressing the bounds on the derivatives J(x) := d j=1 y j ((A x) j + β + w j ) 1 2 (A x + w ) T Σ 1 (A x + w ), α i := max j=1,...,d a ji, d τ := exp 1 d y j aji + (A ) T i Σ 1 w, α i i=1 ω := 1 + Σ 1 1, j=1 where by 1, we mean the entrywise norm of the matrix. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
19 Bounding the norm Theorem The norm of f can be bounded by f 2 τ 2 e r u D R u v u k u γ 1 u (h(0)) 2(d u ) exp( λ min x 2 k + µ kx k + θ k x k ) B v v!ω v 2 i v α i h(x i ) h (x i ) i v i u v 2 dx u. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
20 Bounding the norm where λ min is the smallest eigenvalue of the matrix A T Σ 1 A (which is positive definite), and r := 2 d j=1 y j (β + w j ) w T Σ 1 w, µ k := 2(w T Σ 1 A ) k + 2 θ k := 2 d i=1 d y j ajk, j=1 ((A ) T i Σ 1 A ) k α i. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
21 Logistic density The logistic univariate density is defined by This leads to φ(x) = e x/ν ν(1 + e x/ν ) 2. h(x) = ν(e x/ν e x/ν ) 4νe x /ν, and h (x) = e x/ν e x/ν e x /ν. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
22 Logistic density We then have Denote also h(x i ) h (x i ) (4ν) v e x i /ν. i v i u v i u θ k := θ k + 2 ν. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
23 Bounding the norm for the logistic density Theorem f 2 τ 2 e ( ) u π r γ 1 u (4ν) 2(d u ) B u 2 λ ( u!)2 min u D ) 2 (1 + 4νω 1/2 α k k u ( ( ) ( )) ( θ k µ k ) 2 ( θ k + µ k ) 2 exp + exp. 4λ min 4λ min k u Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
24 On the possible weights We remark that there are two distinctive parts in the RHS of the bound. One part is depending only on the cardinality of the set u and the second one is depending on the coordinates belonging to u. Consequently, the optimal weights might be a hybrid between product weights and order dependent weights. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
25 POD weights POD weights who were first introduced by Kuo, Schwab and Sloan (2011) in connection with the use of QMC methods for a class of elliptic PDE. Such weights are defined by γ u := Γ u γ j, where Γ u relates only to the cardinality of u (the order dependent part), while γ j is a weight associated with the individual coordinate j. It remains to figure out the expression of the optimal weights for our problem. j u Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
26 Worst-case error Here, it was the norm of the integrand that received full attention. We already know that there exist shifted lattice rules such that the worst-case error satisfies en,d,γ sh 1 (z) (ϕ(n)) 2λ γ λ u(ϱ(λ)) u u D 1 2λ, where n is the number of points of the lattice, z is the generating vector, λ ( 1 2, 1] and ϱ is a function depending only on λ. These lattice rules can be constructed via the CBC algorithm. Contributions by Kuo, Nichols, Schwab, Sloan, Wasilkowski, Waterhouse. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
27 Error bound The error bound is given by (ϕ(n)) 1 2λ 1 2λ γ λ u(ϱ(λ)) u u D u D The weights that minimise such error bounds are γ u = ( β u ( ϱ(λ)) u k u γ 1 u b k ) 2 1+λ. β u k u b k 1 2. Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
28 Future research 1 Finding the optimal weights for other statistical models. 2 Develop similar analysis for other important applications, such as finance (work in progress with some preliminary results). Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
29 Thank you! Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February / 24
Quasi-Monte Carlo integration over the Euclidean space and applications
Quasi-Monte Carlo integration over the Euclidean space and applications f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with James Nichols (UNSW) Journal of Complexity 30
More informationHot new directions for QMC research in step with applications
Hot new directions for QMC research in step with applications Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia Frances Kuo @ UNSW Australia p.1 Outline Motivating example
More informationNew Multilevel Algorithms Based on Quasi-Monte Carlo Point Sets
New Multilevel Based on Quasi-Monte Carlo Point Sets Michael Gnewuch Institut für Informatik Christian-Albrechts-Universität Kiel 1 February 13, 2012 Based on Joint Work with Jan Baldeaux (UTS Sydney)
More informationTutorial on quasi-monte Carlo methods
Tutorial on quasi-monte Carlo methods Josef Dick School of Mathematics and Statistics, UNSW, Sydney, Australia josef.dick@unsw.edu.au Comparison: MCMC, MC, QMC Roughly speaking: Markov chain Monte Carlo
More informationApplication of QMC methods to PDEs with random coefficients
Application of QMC methods to PDEs with random coefficients a survey of analysis and implementation Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with Ivan Graham
More informationWhat s new in high-dimensional integration? designing for applications
What s new in high-dimensional integration? designing for applications Ian H. Sloan i.sloan@unsw.edu.au The University of New South Wales UNSW Australia ANZIAM NSW/ACT, November, 2015 The theme High dimensional
More informationLattice rules and sequences for non periodic smooth integrands
Lattice rules and sequences for non periodic smooth integrands dirk.nuyens@cs.kuleuven.be Department of Computer Science KU Leuven, Belgium Joint work with Josef Dick (UNSW) and Friedrich Pillichshammer
More informationQuadrature using sparse grids on products of spheres
Quadrature using sparse grids on products of spheres Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at 3rd Workshop on High-Dimensional Approximation UNSW,
More informationFAST QMC MATRIX-VECTOR MULTIPLICATION
FAST QMC MATRIX-VECTOR MULTIPLICATION JOSEF DICK, FRANCES Y. KUO, QUOC T. LE GIA, CHRISTOPH SCHWAB Abstract. Quasi-Monte Carlo (QMC) rules 1/N N 1 n=0 f(y na) can be used to approximate integrals of the
More informationOptimal Randomized Algorithms for Integration on Function Spaces with underlying ANOVA decomposition
Optimal Randomized on Function Spaces with underlying ANOVA decomposition Michael Gnewuch 1 University of Kaiserslautern, Germany October 16, 2013 Based on Joint Work with Jan Baldeaux (UTS Sydney) & Josef
More informationPart III. Quasi Monte Carlo methods 146/349
Part III Quasi Monte Carlo methods 46/349 Outline Quasi Monte Carlo methods 47/349 Quasi Monte Carlo methods Let y = (y,...,y N ) be a vector of independent uniform random variables in Γ = [0,] N, u(y)
More informationarxiv: v3 [math.na] 18 Sep 2016
Infinite-dimensional integration and the multivariate decomposition method F. Y. Kuo, D. Nuyens, L. Plaskota, I. H. Sloan, and G. W. Wasilkowski 9 September 2016 arxiv:1501.05445v3 [math.na] 18 Sep 2016
More informationAPPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION. I.H. Sloan, X. Wang and H.
APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION I.H. Sloan, X. Wang and H. Wozniakowski Published in Journal of Complexity, Volume 20, Number 1,
More informationQuadrature using sparse grids on products of spheres
Quadrature using sparse grids on products of spheres Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ANZIAM NSW/ACT Annual Meeting, Batemans Bay. Joint
More informationOn the Behavior of the Weighted Star Discrepancy Bounds for Shifted Lattice Rules
On the Behavior of the Weighted Star Discrepancy Bounds for Shifted Lattice Rules Vasile Sinescu and Pierre L Ecuyer Abstract We examine the question of constructing shifted lattice rules of rank one with
More informationFast evaluation of mixed derivatives and calculation of optimal weights for integration. Hernan Leovey
Fast evaluation of mixed derivatives and calculation of optimal weights for integration Humboldt Universität zu Berlin 02.14.2012 MCQMC2012 Tenth International Conference on Monte Carlo and Quasi Monte
More informationQMC methods in quantitative finance. and perspectives
, tradition and perspectives Johannes Kepler University Linz (JKU) WU Research Seminar What is the content of the talk Valuation of financial derivatives in stochastic market models using (QMC-)simulation
More informationAdapting quasi-monte Carlo methods to simulation problems in weighted Korobov spaces
Adapting quasi-monte Carlo methods to simulation problems in weighted Korobov spaces Christian Irrgeher joint work with G. Leobacher RICAM Special Semester Workshop 1 Uniform distribution and quasi-monte
More informationApplicability of Quasi-Monte Carlo for lattice systems
Applicability of Quasi-Monte Carlo for lattice systems Andreas Ammon 1,2, Tobias Hartung 1,2,Karl Jansen 2, Hernan Leovey 3, Andreas Griewank 3, Michael Müller-Preussker 1 1 Humboldt-University Berlin,
More informationANOVA decomposition of convex piecewise linear functions
ANOVA decomposition of convex piecewise linear functions W. Römisch Abstract Piecewise linear convex functions arise as integrands in stochastic programs. They are Lipschitz continuous on their domain,
More informationExponential Convergence and Tractability of Multivariate Integration for Korobov Spaces
Exponential Convergence and Tractability of Multivariate Integration for Korobov Spaces Josef Dick, Gerhard Larcher, Friedrich Pillichshammer and Henryk Woźniakowski March 12, 2010 Abstract In this paper
More informationA Belgian view on lattice rules
A Belgian view on lattice rules Ronald Cools Dept. of Computer Science, KU Leuven Linz, Austria, October 14 18, 2013 Atomium Brussels built in 1958 height 103m figure = 2e coin 5 10 6 in circulation Body
More informationIntegration of permutation-invariant. functions
permutation-invariant Markus Weimar Philipps-University Marburg Joint work with Dirk Nuyens and Gowri Suryanarayana (KU Leuven, Belgium) MCQMC2014, Leuven April 06-11, 2014 un Outline Permutation-invariant
More informationProgress in high-dimensional numerical integration and its application to stochastic optimization
Progress in high-dimensional numerical integration and its application to stochastic optimization W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Short course,
More informationOn Infinite-Dimensional Integration in Weighted Hilbert Spaces
On Infinite-Dimensional Integration in Weighted Hilbert Spaces Sebastian Mayer Universität Bonn Joint work with M. Gnewuch (UNSW Sydney) and K. Ritter (TU Kaiserslautern). HDA 2013, Canberra, Australia.
More informationLow Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed?
Low Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed? Xiaoqun Wang 1,2 and Ian H. Sloan 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084,
More informationSparse grid quadrature on products of spheres
Sparse grid quadrature on products of spheres Markus Hegland a, Paul Leopardi a, a Centre for Mathematics and its Applications, Australian National University. Abstract This paper examines sparse grid
More informationQMC designs and covering of spheres by spherical caps
QMC designs and covering of spheres by spherical caps Ian H. Sloan University of New South Wales, Sydney, Australia Joint work with Johann Brauchart, Josef Dick, Ed Saff (Vanderbilt), Rob Womersley and
More informationIan Sloan and Lattice Rules
www.oeaw.ac.at Ian Sloan and Lattice Rules P. Kritzer, H. Niederreiter, F. Pillichshammer RICAM-Report 2016-34 www.ricam.oeaw.ac.at Ian Sloan and Lattice Rules Peter Kritzer, Harald Niederreiter, Friedrich
More informationInformation geometry for bivariate distribution control
Information geometry for bivariate distribution control C.T.J.Dodson + Hong Wang Mathematics + Control Systems Centre, University of Manchester Institute of Science and Technology Optimal control of stochastic
More informationDigital lattice rules
Digital lattice rules Multivariate integration and discrepancy estimates A thesis presented to The University of New South Wales in fulfilment of the thesis requirement for the degree of Doctor of Philosophy
More informationQuasi-Monte Carlo for Highly Structured Generalised Response Models
Quasi-Monte Carlo for Highly Structured Generalised Response Models BY F.Y. KUO, W.T.M. DUNSMUIR, I.H. SLOAN, M.P. WAND AND R.S. WOMERSLEY School of Mathematics and Statistics, University of New South
More informationEffective dimension of some weighted pre-sobolev spaces with dominating mixed partial derivatives
Effective dimension of some weighted pre-sobolev spaces with dominating mixed partial derivatives Art B. Owen Stanford University August 2018 Abstract This paper considers two notions of effective dimension
More information2 FRED J. HICKERNELL the sample mean of the y (i) : (2) ^ 1 N The mean square error of this estimate may be written as a sum of two parts, a bias term
GOODNESS OF FIT STATISTICS, DISCREPANCIES AND ROBUST DESIGNS FRED J. HICKERNELL Abstract. The Cramer{Von Mises goodness-of-t statistic, also known as the L 2 -star discrepancy, is the optimality criterion
More informationLow-discrepancy point sets lifted to the unit sphere
Low-discrepancy point sets lifted to the unit sphere Johann S. Brauchart School of Mathematics and Statistics, UNSW MCQMC 2012 (UNSW) February 17 based on joint work with Ed Saff (Vanderbilt) and Josef
More informationKernel Principal Component Analysis
Kernel Principal Component Analysis Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 Outlines Overview Introduction Linear Algebra Probability Linear Regression
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationA MONTE CARLO ALGORITHM FOR WEIGHTED INTEGRATION OVER R d
MATHEMATICS OF COMPUTATION Volume 73, Number 246, Pages 813 825 S 25-5718(3)1564-3 Article electronically published on August 19, 23 A MONTE CARLO ALGORITHM FOR WEIGHTED INTEGRATION OVER R d PIOTR GAJDA,
More informationConvergence Rates of Kernel Quadrature Rules
Convergence Rates of Kernel Quadrature Rules Francis Bach INRIA - Ecole Normale Supérieure, Paris, France ÉCOLE NORMALE SUPÉRIEURE NIPS workshop on probabilistic integration - Dec. 2015 Outline Introduction
More informationScenario generation in stochastic programming with application to energy systems
Scenario generation in stochastic programming with application to energy systems W. Römisch Humboldt-University Berlin Institute of Mathematics www.math.hu-berlin.de/~romisch SESO 2017 International Thematic
More informationTractability of Multivariate Problems
Erich Novak University of Jena Chemnitz, Summer School 2010 1 Plan for the talk Example: approximation of C -functions What is tractability? Tractability by smoothness? Tractability by sparsity, finite
More informationAPPLIED MATHEMATICS REPORT AMR04/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR. J. Dick and F. Pillichshammer
APPLIED MATHEMATICS REPORT AMR4/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR J. Dick and F. Pillichshammer January, 24 Diaphony, discrepancy, spectral test and worst-case error Josef Dick
More informationIntegral representations in models with long memory
Integral representations in models with long memory Georgiy Shevchenko, Yuliya Mishura, Esko Valkeila, Lauri Viitasaari, Taras Shalaiko Taras Shevchenko National University of Kyiv 29 September 215, Ulm
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 12, 2007 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 9, 2011 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationReproducing Kernel Hilbert Spaces Class 03, 15 February 2006 Andrea Caponnetto
Reproducing Kernel Hilbert Spaces 9.520 Class 03, 15 February 2006 Andrea Caponnetto About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationError Analysis for Quasi-Monte Carlo Methods
Error Analysis for Quasi-Monte Carlo Methods Fred J. Hickernell Department of Applied Mathematics, Illinois Institute of Technology hickernell@iit.edu mypages.iit.edu/~hickernell Thanks to the Guaranteed
More informationInfinite-Dimensional Integration on Weighted Hilbert Spaces
Infinite-Dimensional Integration on Weighted Hilbert Spaces Michael Gnewuch Department of Computer Science, Columbia University, 1214 Amsterdam Avenue, New Yor, NY 10027, USA Abstract We study the numerical
More informationOn the existence of higher order polynomial. lattices with large figure of merit ϱ α.
On the existence of higher order polynomial lattices with large figure of merit Josef Dick a, Peter Kritzer b, Friedrich Pillichshammer c, Wolfgang Ch. Schmid b, a School of Mathematics, University of
More informationSTRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE CARLO ALGORITHMS
MATHEMATICS OF COMPUTATION Volume 72, Number 242, Pages 823 838 S 0025-5718(02)01440-0 Article electronically published on June 25, 2002 STRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE
More informationSupport points. May 26, Abstract
Support points Simon Mak, V. Roshan Joseph May 26, 2018 arxiv:1609.01811v6 [math.st] 8 Sep 2017 Abstract This paper introduces a new way to compact a continuous probability distribution F into a set of
More informationHigh dimensional integration of kinks and jumps smoothing by preintegration. Andreas Griewank, Frances Y. Kuo, Hernan Leövey and Ian H.
High dimensional integration of kinks and jumps smoothing by preintegration Andreas Griewank, Frances Y. Kuo, Hernan Leövey and Ian H. Sloan August 2017 arxiv:1712.00920v1 [math.na] 4 Dec 2017 Abstract
More informationStochastic optimization, multivariate numerical integration and Quasi-Monte Carlo methods. W. Römisch
Stochastic optimization, multivariate numerical integration and Quasi-Monte Carlo methods W. Römisch Humboldt-University Berlin Institute of Mathematics www.math.hu-berlin.de/~romisch Chemnitzer Mathematisches
More informationON TRACTABILITY OF WEIGHTED INTEGRATION OVER BOUNDED AND UNBOUNDED REGIONS IN R s
MATHEMATICS OF COMPTATION Volume 73, Number 248, Pages 1885 1901 S 0025-5718(04)01624-2 Article electronically published on January 5, 2004 ON TRACTABILITY OF WEIGHTED INTEGRATION OVER BONDED AND NBONDED
More informationLifting the Curse of Dimensionality
Lifting the Curse of Dimensionality Frances Y. Kuo and Ian H. Sloan Introduction Richard Bellman [1] coined the phrase the curse of dimensionality to describe the extraordinarily rapid growth in the difficulty
More informationDifferential Stein operators for multivariate continuous distributions and applications
Differential Stein operators for multivariate continuous distributions and applications Gesine Reinert A French/American Collaborative Colloquium on Concentration Inequalities, High Dimensional Statistics
More informationConstruction Algorithms for Higher Order Polynomial Lattice Rules
Construction Algorithms for Higher Order Polynomial Lattice Rules Jan Baldeaux, Josef Dick, Julia Greslehner and Friedrich Pillichshammer Dedicated to Gerhard Larcher on the occasion of his 50th birthday
More informationERROR ESTIMATES FOR THE ANOVA METHOD WITH POLYNOMIAL CHAOS INTERPOLATION: TENSOR PRODUCT FUNCTIONS
ERROR ESTIMATES FOR THE ANOVA METHOD WITH POLYNOMIAL CHAOS INTERPOLATION: TENSOR PRODUCT FUNCTIONS ZHONGQIANG ZHANG, MINSEOK CHOI AND GEORGE EM KARNIADAKIS Abstract. We focus on the analysis of variance
More informationApproximation of High-Dimensional Numerical Problems Algorithms, Analysis and Applications
Approximation of High-Dimensional Numerical Problems Algorithms, Analysis and Applications Christiane Lemieux (University of Waterloo), Ian Sloan (University of New South Wales), Henryk Woźniakowski (University
More informationA Multilevel Method for the Global Optimization of Potential Energy Functions
A Multilevel Method for the Global Optimization of Potential Energy Functions Nadine Pastuszka Born 7th June 1989 in Brühl, Germany 29th April 2014 Master s Thesis Mathematics Advisor: Prof. Dr. Michael
More informationNumerical Methods II
Numerical Methods II Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute MC Lecture 13 p. 1 Quasi-Monte Carlo As in lecture 6, quasi-monte Carlo methods offer much greater
More informationCIS 520: Machine Learning Oct 09, Kernel Methods
CIS 520: Machine Learning Oct 09, 207 Kernel Methods Lecturer: Shivani Agarwal Disclaimer: These notes are designed to be a supplement to the lecture They may or may not cover all the material discussed
More informationNEW FUNCTIONAL INEQUALITIES
1 / 29 NEW FUNCTIONAL INEQUALITIES VIA STEIN S METHOD Giovanni Peccati (Luxembourg University) IMA, Minneapolis: April 28, 2015 2 / 29 INTRODUCTION Based on two joint works: (1) Nourdin, Peccati and Swan
More informationUniform Distribution and Quasi-Monte Carlo Methods
Workshop on Uniform Distribution and Quasi-Monte Carlo Methods October 14-18, 2013 as part of the Radon Special Semester 2013 on Applications of Algebra and Number Theory Metric number theory, lacunary
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 11, 2009 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationDuality of multiparameter Hardy spaces H p on spaces of homogeneous type
Duality of multiparameter Hardy spaces H p on spaces of homogeneous type Yongsheng Han, Ji Li, and Guozhen Lu Department of Mathematics Vanderbilt University Nashville, TN Internet Analysis Seminar 2012
More informationLecture 4 February 2
4-1 EECS 281B / STAT 241B: Advanced Topics in Statistical Learning Spring 29 Lecture 4 February 2 Lecturer: Martin Wainwright Scribe: Luqman Hodgkinson Note: These lecture notes are still rough, and have
More informationarxiv: v2 [math.na] 15 Aug 2018
Noname manuscript No. (will be inserted by the editor) Lattice rules in non-periodic subspaces of Sobolev spaces Takashi Goda Kosuke Suzuki Takehito Yoshiki arxiv:72.2572v2 [math.na] 5 Aug 28 Received:
More information13 Notes on Markov Chain Monte Carlo
13 Notes on Markov Chain Monte Carlo Markov Chain Monte Carlo is a big, and currently very rapidly developing, subject in statistical computation. Many complex and multivariate types of random data, useful
More informationApproximation of Functions Using Digital Nets
Approximation of Functions Using Digital Nets Josef Dick 1, Peter Kritzer 2, and Frances Y. Kuo 3 1 UNSW Asia, 1 Kay Siang Road, Singapore 248922. E-mail: j.dick@unswasia.edu.sg 2 Fachbereich Mathematik,
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationWeighted Geometric Discrepancies and Numerical Integration on Reproducing Kernel Hilbert Spaces
Weighted Geometric Discrepancies and Numerical Integration on Reproducing Kernel Hilbert Spaces ichael Gnewuch Institut für Informatik, Christian-Albrechts-Universität Kiel, Christian-Albrechts-Platz 4,
More informationThe ANOVA decomposition of a non-smooth function of infinitely many variables can have every term smooth
Wegelerstraße 6 53115 Bonn Germany phone +49 228 73-3427 fax +49 228 73-7527 www.ins.uni-bonn.de M. Griebel, F.Y. Kuo and I.H. Sloan The ANOVA decomposition of a non-smooth function of infinitely many
More informationCan we do statistical inference in a non-asymptotic way? 1
Can we do statistical inference in a non-asymptotic way? 1 Guang Cheng 2 Statistics@Purdue www.science.purdue.edu/bigdata/ ONR Review Meeting@Duke Oct 11, 2017 1 Acknowledge NSF, ONR and Simons Foundation.
More informationPolynomial interpolation on the sphere, reproducing kernels and random matrices
Polynomial interpolation on the sphere, reproducing kernels and random matrices Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at MASCOS Workshop on Stochastics
More informationKernel methods and the exponential family
Kernel methods and the exponential family Stéphane Canu 1 and Alex J. Smola 2 1- PSI - FRE CNRS 2645 INSA de Rouen, France St Etienne du Rouvray, France Stephane.Canu@insa-rouen.fr 2- Statistical Machine
More informationarxiv: v7 [math.st] 9 Sep 2018
Submitted to the Annals of Statistics SUPPORT POINTS By Simon Mak and V. Roshan Joseph Georgia Institute of Technology arxiv:1609.01811v7 [math.st] 9 Sep 2018 Abstract This paper introduces a new way to
More informationSparse Approximation via Penalty Decomposition Methods
Sparse Approximation via Penalty Decomposition Methods Zhaosong Lu Yong Zhang February 19, 2012 Abstract In this paper we consider sparse approximation problems, that is, general l 0 minimization problems
More informationCarleson Measures for Besov-Sobolev Spaces and Non-Homogeneous Harmonic Analysis
Carleson Measures for Besov-Sobolev Spaces and Non-Homogeneous Harmonic Analysis Brett D. Wick Georgia Institute of Technology School of Mathematics & Humboldt Fellow Institut für Mathematik Universität
More informationLecture 10: Support Vector Machine and Large Margin Classifier
Lecture 10: Support Vector Machine and Large Margin Classifier Applied Multivariate Analysis Math 570, Fall 2014 Xingye Qiao Department of Mathematical Sciences Binghamton University E-mail: qiao@math.binghamton.edu
More informationIntegration by Substitution
November 22, 2013 Introduction 7x 2 cos(3x 3 )dx =? 2xe x2 +5 dx =? Chain rule The chain rule: d dx (f (g(x))) = f (g(x)) g (x). Use the chain rule to find f (x) and then write the corresponding anti-differentiation
More information6.2 Fubini s Theorem. (µ ν)(c) = f C (x) dµ(x). (6.2) Proof. Note that (X Y, A B, µ ν) must be σ-finite as well, so that.
6.2 Fubini s Theorem Theorem 6.2.1. (Fubini s theorem - first form) Let (, A, µ) and (, B, ν) be complete σ-finite measure spaces. Let C = A B. Then for each µ ν- measurable set C C the section x C is
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationRANDOM FIELDS AND GEOMETRY. Robert Adler and Jonathan Taylor
RANDOM FIELDS AND GEOMETRY from the book of the same name by Robert Adler and Jonathan Taylor IE&M, Technion, Israel, Statistics, Stanford, US. ie.technion.ac.il/adler.phtml www-stat.stanford.edu/ jtaylor
More informationConvergence guarantees for kernel-based quadrature rules in misspecified settings
Convergence guarantees for kernel-based quadrature rules in misspecified settings Motonobu Kanagawa, Bharath K Sriperumbudur, Kenji Fukumizu The Institute of Statistical Mathematics, Tokyo 19-856, Japan
More informationDifferential Equations Preliminary Examination
Differential Equations Preliminary Examination Department of Mathematics University of Utah Salt Lake City, Utah 84112 August 2007 Instructions This examination consists of two parts, called Part A and
More informationA PIECEWISE CONSTANT ALGORITHM FOR WEIGHTED L 1 APPROXIMATION OVER BOUNDED OR UNBOUNDED REGIONS IN R s
A PIECEWISE CONSTANT ALGORITHM FOR WEIGHTED L 1 APPROXIMATION OVER BONDED OR NBONDED REGIONS IN R s FRED J. HICKERNELL, IAN H. SLOAN, AND GRZEGORZ W. WASILKOWSKI Abstract. sing Smolyak s construction [5],
More informationarxiv: v2 [math.na] 15 Oct 2015
Optimal quasi-monte Carlo rules on order 2 digital nets for the numerical integration of multivariate periodic functions Aicke Hinrichs a, Lev Markhasin b, Jens Oettershagen c, Tino Ullrich c arxiv:1501.01800v2
More informationBayesian Models for Regularization in Optimization
Bayesian Models for Regularization in Optimization Aleksandr Aravkin, UBC Bradley Bell, UW Alessandro Chiuso, Padova Michael Friedlander, UBC Gianluigi Pilloneto, Padova Jim Burke, UW MOPTA, Lehigh University,
More informationConditional Sampling for Max Stable Random Fields
Conditional Sampling for Max Stable Random Fields Yizao Wang Department of Statistics, the University of Michigan April 30th, 0 th GSPC at Duke University, Durham, North Carolina Joint work with Stilian
More informationKernels A Machine Learning Overview
Kernels A Machine Learning Overview S.V.N. Vishy Vishwanathan vishy@axiom.anu.edu.au National ICT of Australia and Australian National University Thanks to Alex Smola, Stéphane Canu, Mike Jordan and Peter
More informationFOCM Workshop B5 Information Based Complexity
FOCM 2014 - Workshop B5 Information Based Complexity B5 - December 15, 14:30 15:00 The ANOVA decomposition of a non-smooth function of an infinite number of variables Ian Sloan University of New South
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationCMS winter meeting 2008, Ottawa. The heat kernel on connected sums
CMS winter meeting 2008, Ottawa The heat kernel on connected sums Laurent Saloff-Coste (Cornell University) December 8 2008 p(t, x, y) = The heat kernel ) 1 x y 2 exp ( (4πt) n/2 4t The heat kernel is:
More informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationA Practical Guide to Quasi-Monte Carlo Methods
A Practical Guide to Quasi-Monte Carlo Methods Frances Y. Kuo and Dirk Nuyens Abstract These notes are prepared for the short course on High-dimensional Integration: the Quasi-Monte Carlo Way, to be held
More informationLiberating the Dimension for Function Approximation and Integration
Liberating the Dimension for Function Approximation and Integration G W Wasilkowski Abstract We discuss recent results on the complexity and tractability of problems dealing with -variate functions Such
More informationMax stable Processes & Random Fields: Representations, Models, and Prediction
Max stable Processes & Random Fields: Representations, Models, and Prediction Stilian Stoev University of Michigan, Ann Arbor March 2, 2011 Based on joint works with Yizao Wang and Murad S. Taqqu. 1 Preliminaries
More informationSmall ball probabilities and metric entropy
Small ball probabilities and metric entropy Frank Aurzada, TU Berlin Sydney, February 2012 MCQMC Outline 1 Small ball probabilities vs. metric entropy 2 Connection to other questions 3 Recent results for
More information