Hypothesis testing for Stochastic PDEs. Igor Cialenco
|
|
- Matthew McCarthy
- 6 years ago
- Views:
Transcription
1 Hypothesis testing for Stochastic PDEs Igor Cialenco Department of Applied Mathematics Illinois Institute of Technology Joint work with Liaosha Xu Research partially funded by NSF grants DMS Asymptotical Statistics of Stochastic Processes X Le Mans, March, 2015
2 Outline of the talk Overview of Parameter Estimation for SPDEs Hypothesis testing Error estimation Numerical results Ig. Cialenco, IIT March 17-20, 2015 Slide # 2
3 Bibliography I. Cialenco, L. Xu, Hypothesis testing for stochastic PDEs driven by additive noise, Stochastic Processes and their Applications, vol. 125, Issue 3, March 2015, pp I. Cialenco, L. Xu, A note on error estimation for hypothesis testing problems for some linear SPDEs, Stochastic Partial Differential Equations: Analysis and Computations, September 2014, vol. 2, No 3, pp Ig. Cialenco, IIT March 17-20, 2015 Slide # 3
4 The Setup Let (Ω, F, {F t } t 0, P) be a stochastic basis with usual assumptions. Consider heat equation driven by additive noise: du(t, x) θ U(t, x)dt = σ k=1 where x G, G is a bounded domain in R d, t [0, T ]; λ γ k h k(x)dw k (t), (0.1) θ > 0 is the parameter of interest; γ 0, σ R {0} are known. zero initial conditions and boundary values; {w k (t)} k N are independent Brownian motions; is the Laplace operator on G with zero boundary condition U = d j=1 2 U ; x 2 j {h k } are the eigenfunctions of in L 2 (G), and {ρ k } are the corresponding eigenvalues; λ k = ρ k k 1/d ; Ig. Cialenco, IIT March 17-20, 2015 Slide # 4
5 Existence of the solution Theorem (Existence and Uniqueness) If (γ s)/d > 1/2, then the SPDE (0.1) has a unique solution (weak in PDE sense, and strong in probability sense) U L 2 (Ω [0, T ]; H s+1 (G)) L 2 (Ω; C((0, T ); H s (G))). Ig. Cialenco, IIT March 17-20, 2015 Slide # 5
6 Construction of the test Part I: Hypothesis Testing for SPDEs Ig. Cialenco, IIT March 17-20, 2015 Slide # 6
7 The Problem: Simple Hypothesis Hypothesis Testing Construction of the test du(t, x) θ U(t, x)dt = σ k=1 Assume that θ can take only two values {, θ 1 }. Consider a simple hypothesis: H 0 θ =, H 1 θ = θ 1. For simplicity, assume θ 1 > and σ > 0. λ γ k h k(x)dw k (t). Ig. Cialenco, IIT March 17-20, 2015 Slide # 7
8 Construction of the Test Hypothesis Testing Construction of the test du(t, x) θ U(t, x)dt = σ k=1 λ γ k h k(x)dw k (t), U(0, x) = 0. The k-th Fourier coefficient u k (t) = U(t, x), h k (x) is given by du k = θλ 2 k u kdt + σλ γ k dw k(t), u k (0) = 0, u k (t) = σλ γ k t 0 e θλ2 k (t s) dw k, k 1. Let P N,T θ ( ) = P(U N T ) be the measure on C([0, T ]; RN ) generated by U N T (t) = (u 1,..., u N ) up to time T. Ig. Cialenco, IIT March 17-20, 2015 Slide # 8
9 Construction of the test Observable: First N Fourier coefficients, for all t [0, T ], U N T = ((u 1 (t),..., u N (t)), t [0, T ]). In this case, the Likelihood Ratio has the form L(, θ 1, UT N ) = PN,T θ 1 P N,T θ 1 ( 0 = exp ( (θ 1 )σ 2 T N λ 2+2γ k k=1 u k (t)du k (t) (θ 1 + )λ 2 k and the Maximum Likelihood Estimator is 0 T u 2 k (t)dt)), θ N T = N k=1 λ2+2γ k N k=1 λ4+2γ k T 0 u k(t)du k (t) T 0 u2 k (t)dt, N N, T > 0. Ig. Cialenco, IIT March 17-20, 2015 Slide # 9
10 Construction of the test Two Asymptotic Regimes: Large times, T ; Large number of Fourier modes, N. Theorem Assume that 2γ > d. Then, lim θ T N N = lim T θn T = θ, a.e. and where M = N k=1 λ2 k. (4/d + 2)θ ϖσ 2 T ). lim N 1/d+ 1 2 ( θ T N θ) = d N (0, N T ( θn T θ) = d N (0, 2θ/M), lim T Ig. Cialenco, IIT March 17-20, 2015 Slide # 10
11 Construction of the test Looking for rejection region R B(C([0, T ]; R N )). Type I error = P N,T (R); Type II error = 1 P N,T θ 1 Define the class of test (R), and power of the test = P N,T θ 1 (R) K α = {R B(C([0, T ]; R N )) P N,T (R) α}. with α (0, 1) being the significant level, fixed in what follows. Definition We say that a rejection region R K α is the most powerful in the class K α if P N,T θ 1 (R) P N,T θ 1 (R ), for all R K α. Ig. Cialenco, IIT March 17-20, 2015 Slide # 11
12 Neyman-Pearson Lemma Hypothesis Testing Most Powerful Test Theorem (C. and Xu, 14) Let c α be a real number such that P N,T (L(, θ 1, U N T ) c α ) = α. Then, R = {U N T L(, θ 1, U N T ) c α }, is the most powerful rejection region in the class K α. Ig. Cialenco, IIT March 17-20, 2015 Slide # 12
13 The Difficulty: c α has no explicit formula for finite T and N. Ig. Cialenco, IIT March 17-20, 2015 Slide # 13
14 The Difficulty: c α has no explicit formula for finite T and N. We take/suggest an Asymptotic Method (1) Fix N, let T ; (2) Fix T, let N. Inspired by large times Asymptotic Method for finite dimensional ergodic diffusion processes developed by Kutoyants 74, and 04. In this talk we focus on case (2), asymptotics in number of Fourier modes; For case (1) see [CX 14]. Ig. Cialenco, IIT March 17-20, 2015 Slide # 13
15 umber of Fourier modes N Define a new class K α = {(R N ) N N R N B(C([0, T ]; R N ), lim sup P N,T (R N ) α}, N where T is fixed, and α is the Asymptotic Significant Level. Ig. Cialenco, IIT March 17-20, 2015 Slide # 14
16 umber of Fourier modes N Define a new class K α = {(R N ) N N R N B(C([0, T ]; R N ), lim sup P N,T (R N ) α}, N where T is fixed, and α is the Asymptotic Significant Level. Goal: We want to find a rejection region ( R N ) N N such that lim N P N,T ( R N ) = α. Ig. Cialenco, IIT March 17-20, 2015 Slide # 14
17 Attempt: We still try Likelihood Ratio test. Then, what is c α? P N,T = P N,T where (L(, θ 1, U N T ) c α (N)) XN Y N X N = 8θ 3 0 ln c α (N) T M(θ 2 1 θ0 2) + 2θ0 T M(θ 1 ) 2(θ 1 + ) N θ0 (θ 1 ) σ 2 2T M(θ 1 + ) ( k=1 2θ0 T Y N = σ T M M = N k=1 λ 2 k. N λ 2+γ k k=1 0 λ 2+2γ k u 2 σ2 k (T ) N), 2 u k dw k, (θ 1 )N 8θ0 T M(θ 1 + ), Ig. Cialenco, IIT March 17-20, 2015 Slide # 15
18 For δ R, by splitting argument, we get P N,T + P N,T (L(, θ 1, UT N ) c α (N)) P N,T (X N δ) Y N 8θ 3 0 ln c α (N) T M(θ 2 1 θ0 2) + 2θ0 T M(θ 1 ) 2(θ 1 + ) (θ 1 )N 8θ0 T M(θ 1 + ) δ. Take δ = N 1/2d, and deduce that P N,T (X N δ) 0, N. Use the fact that Y N d N (0, 1) as N, and find Reasonable candidate for the threshold constant ĉ α (N) = exp (θ 1 ) 2 T M 4 + (θ 1 ) 2 N 8θ 2 0 T M(θ 2 1 θ0 2 ) q α 8θ 3 0. Ig. Cialenco, IIT March 17-20, 2015 Slide # 16
19 Theorem (C. and Xu) Suppose R N = {U N T L(, θ 1, U N T ) ĉ α (N)}, for all N, where ĉ α is given above. Then, the rejection region ( R N ) N N K α, and moreover lim N PN,T ( R N ) = α. Ig. Cialenco, IIT March 17-20, 2015 Slide # 17
20 The Next Question: Is ( R N ) N N (asymptotically) the most powerful test in K α? Ig. Cialenco, IIT March 17-20, 2015 Slide # 18
21 The Next Question: Is ( R N ) N N (asymptotically) the most powerful test in K α? Short answer: No! Ig. Cialenco, IIT March 17-20, 2015 Slide # 18
22 The Next Question: Is ( R N ) N N (asymptotically) the most powerful test in K α? Short answer: No! How does the power of this test P N,T θ 1 ( R N ) behave? Definition We say that a rejection region (R N ) K α is asymptotically the most powerful in the class K α if 1 P N,T θ lim inf 1 (R N ) N 1 P N,T θ 1 (RN ) 1, for all (R N) Kα. Similarly, we define asymptotically the most powerful rejection regions for a different given class of tests. Ig. Cialenco, IIT March 17-20, 2015 Slide # 18
23 Theorem (C. and Xu) (Criterion for Most Powerful Test) Consider the rejection region of the form R N = {U N T L(, θ 1, U N T ) c α(n)}, where c α(n) is a function of N such that, c α(n) > 0 for all N N, and lim N PN,T (R N) = α, lim N c α(n) 1 P N,T θ 1 (R N ) <. Then (R N ) is asymptotically the most powerful in K α. Is the test R N asymptotically the most powerful in the class K α? We need to control the power of the test P N,T θ 1 ( R N ), as N. Ig. Cialenco, IIT March 17-20, 2015 Slide # 19
24 Calculate the Cumulant Generating Function of the Log-Likelihood ratio m(ɛ) = E[exp(ɛ ln L(, θ 1, U N T ))] m(ɛ) = exp [ 1 2 N k=1 ln (cosh(γ k T ) p sinh(γ k T )) + Similar to Gapeev and Küchler [2008] Use Feynman-Kac Formula to derive a PDE Make some transforms and guess the solution ɛ(θ1 θ0) + θ1 MT ]. 2 Sharp Large Deviations (see Lin kov 1999, and Bercu and Roualt 2001 for T ) L T (ɛ) = T 1 ln E θ1 [exp (ɛ ln L(, θ 1, U N T ))] = L(ɛ) + T 1 H(ɛ) + T 1 R T (ɛ), L N (ɛ) = M 1 ln E θ1 [exp (ɛ ln L(, θ 1, U N T ))] = L(ɛ) + NM 1 H(ɛ) + M 1 RN (ɛ), N M = λ 2 k N 2/d+1. k=1 Ig. Cialenco, IIT March 17-20, 2015 Slide # 20
25 Corollary After a series of technical results, we get ĉ α (N)/ (1 P N,T θ 1 ( R N )) M N 1/d+1/2, and hence R N is NOT assymptoticaly most powerful in the class K α. Ig. Cialenco, IIT March 17-20, 2015 Slide # 21
26 Corollary After a series of technical results, we get ĉ α (N)/ (1 P N,T θ 1 ( R N )) M N 1/d+1/2, and hence R N is NOT assymptoticaly most powerful in the class K α. Solution: a new class of tests Similar remark is true for T : c α(t )/ (1 P N,T θ 1 (R T )) T, as T, where R T ={U N T L(, θ 1, U N T ) c α(t )}, for all T, c α(t ) = exp (θ1 θ0)2 MT θ2 1 θ0 2 MT q α Ig. Cialenco, IIT March 17-20, 2015 Slide # 21
27 Control of Type I error Hypothesis Testing Proposition (C. and Xu) For any x R, we have the following expansion, P N,T (I N x) = Φ(x) + Φ δ 1(x)M 1/2 + Φ δ 2(x)NM 1 + R δ N(x) (M 1 + NM 3/2 + N 2 M 2 ), where I N = 8θ 3 0 ln L(,θ 1,U N T ) T M(θ 2 1 θ 2 0 ) 2θ0 T M(θ 1 ) 2(θ 1 + ) + (θ 1 )N 8θ0 T M(θ 1 + ), and Φ( ) is the distribution of a standard Gaussian random variable, and Φ 1 ( ) and Φ 2 ( ) are bounded smooth functions and all their derivatives are bounded, and R N ( ) is uniformly bounded in x. Ig. Cialenco, IIT March 17-20, 2015 Slide # 22
28 Main result Hypothesis Testing Define ˆK α = {(R N ) lim sup (P N,T (R N ) α) M α 1 }. N Note that R N = {U N T IN q α }, so by the previous lemma we have where α 1 is an explicit constant. lim N (PN,T ( ˆR N ) α) M = ˆα 1, Ig. Cialenco, IIT March 17-20, 2015 Slide # 23
29 Main result Hypothesis Testing Define ˆK α = {(R N ) lim sup (P N,T (R N ) α) M α 1 }. N Note that R N = {U N T IN q α }, so by the previous lemma we have where α 1 is an explicit constant. Theorem (C. and Xu) lim N (PN,T ( ˆR N ) α) M = ˆα 1, The rejection region ( ˆR N ) is asymptotically the most powerful in ˆK α. Ig. Cialenco, IIT March 17-20, 2015 Slide # 23
30 Part II: Error Control and Numerical Results Ig. Cialenco, IIT March 17-20, 2015 Slide # 24
31 New Tests for Error Control, N Hypothesis Testing New Tests for Error Control Theorem (C. and Xu) Consider the test statistics of the form R 0 N = {U N T ln L(, θ 1, U N T ) ζ 0 M}, where ζ 0 is given by an explicit formula of the form (θ 1 ) 2 4 T + O(N 1/2 1/d ). If N N 0 (N 0 has explicit formula), then the Type I and Type II errors have the following bound estimates P N,T (RN 0 ) (1 + ϱ)α, 1 P N,T θ 1 (RN 0 ) (1 + ϱ) exp ( (θ 1 ) 2 16θ0 2 MT ), where ϱ denotes a given threshold of error tolerance. Ig. Cialenco, IIT March 17-20, 2015 Slide # 25
32 Table: Type I errors for various N N P N,T (RN 0 ) P N,T (RN ) Other parameters: α = 0.05, = 0.1, θ 1 = 0.2, T = 1, ϱ = 0.1, d = σ = 1, γ = 0
33 Table: T = Tb 1 given by theoretical results and Type I error for various α. α T 1 b P N,T (RT 0 ) Other parameters: = 0.1, θ 1 = 0.2, N = 3, ϱ = 0.1, d = σ = 1, γ = 0
34 Table: Type II errors for various T. T exp ( (θ 1 ) 2 16 MT ) P N,T θ 1 (R 0 T ) P N,T θ 1 (R T ) Other parameters: α = 0.05, = 0.1, θ 1 = 0.2, N = 3, ϱ = 0.1, d = σ = 1, γ = 0
35 Thank You! The end of the talk... but not of the story
A note on error estimation for hypothesis testing problems for some linear SPDEs
Stoch PDE: Anal Comp 214 2:48 431 DOI 1.17/s472-14-36-4 A note on error estimation for hypothesis testing problems for some linear SPDEs Igor Cialenco Liaosha Xu Received: 17 December 213 / Published online:
More informationIgor Cialenco. Department of Applied Mathematics Illinois Institute of Technology, USA joint with N.
Parameter Estimation for Stochastic Navier-Stokes Equations Igor Cialenco Department of Applied Mathematics Illinois Institute of Technology, USA igor@math.iit.edu joint with N. Glatt-Holtz (IU) Asymptotical
More informationBrownian Motion and Stochastic Calculus
ETHZ, Spring 17 D-MATH Prof Dr Martin Larsson Coordinator A Sepúlveda Brownian Motion and Stochastic Calculus Exercise sheet 6 Please hand in your solutions during exercise class or in your assistant s
More informationParameter Estimation for Stochastic Evolution Equations with Non-commuting Operators
Parameter Estimation for Stochastic Evolution Equations with Non-commuting Operators Sergey V. Lototsky and Boris L. Rosovskii In: V. Korolyuk, N. Portenko, and H. Syta (editors), Skorokhod s Ideas in
More informationEconomics 520. Lecture Note 19: Hypothesis Testing via the Neyman-Pearson Lemma CB 8.1,
Economics 520 Lecture Note 9: Hypothesis Testing via the Neyman-Pearson Lemma CB 8., 8.3.-8.3.3 Uniformly Most Powerful Tests and the Neyman-Pearson Lemma Let s return to the hypothesis testing problem
More informationBrownian survival and Lifshitz tail in perturbed lattice disorder
Brownian survival and Lifshitz tail in perturbed lattice disorder Ryoki Fukushima Kyoto niversity Random Processes and Systems February 16, 2009 6 B T 1. Model ) ({B t t 0, P x : standard Brownian motion
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationAsymptotic results for empirical measures of weighted sums of independent random variables
Asymptotic results for empirical measures of weighted sums of independent random variables B. Bercu and W. Bryc University Bordeaux 1, France Workshop on Limit Theorems, University Paris 1 Paris, January
More informationAsymptotic Properties of an Approximate Maximum Likelihood Estimator for Stochastic PDEs
Asymptotic Properties of an Approximate Maximum Likelihood Estimator for Stochastic PDEs M. Huebner S. Lototsky B.L. Rozovskii In: Yu. M. Kabanov, B. L. Rozovskii, and A. N. Shiryaev editors, Statistics
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationLarge deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationSTATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, = (2n 1)(2n 3) 3 1.
STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, 26 Problem Normal Moments (A) Use the Itô formula and Brownian scaling to check that the even moments of the normal distribution
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationMarkov Chain BSDEs and risk averse networks
Markov Chain BSDEs and risk averse networks Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) 2nd Young Researchers in BSDEs
More informationLecture 8: Information Theory and Statistics
Lecture 8: Information Theory and Statistics Part II: Hypothesis Testing and I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 23, 2015 1 / 50 I-Hsiang
More informationMore Empirical Process Theory
More Empirical Process heory 4.384 ime Series Analysis, Fall 2008 Recitation by Paul Schrimpf Supplementary to lectures given by Anna Mikusheva October 24, 2008 Recitation 8 More Empirical Process heory
More informationTopic 19 Extensions on the Likelihood Ratio
Topic 19 Extensions on the Likelihood Ratio Two-Sided Tests 1 / 12 Outline Overview Normal Observations Power Analysis 2 / 12 Overview The likelihood ratio test is a popular choice for composite hypothesis
More informationLecture 12 November 3
STATS 300A: Theory of Statistics Fall 2015 Lecture 12 November 3 Lecturer: Lester Mackey Scribe: Jae Hyuck Park, Christian Fong Warning: These notes may contain factual and/or typographic errors. 12.1
More informationCorrelation Detection and an Operational Interpretation of the Rényi Mutual Information
Correlation Detection and an Operational Interpretation of the Rényi Mutual Information Masahito Hayashi 1, Marco Tomamichel 2 1 Graduate School of Mathematics, Nagoya University, and Centre for Quantum
More informationOther properties of M M 1
Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected
More informationRegularity of Weak Solution to Parabolic Fractional p-laplacian
Regularity of Weak Solution to Parabolic Fractional p-laplacian Lan Tang at BCAM Seminar July 18th, 2012 Table of contents 1 1. Introduction 1.1. Background 1.2. Some Classical Results for Local Case 2
More informationAsymptotic analysis of the sieve estimator for a class of parabolic SPDEs
Asymptotic analysis of the sieve estimator for a class of parabolic SPDEs M. Huebner Michigan State University S. Lototsky Massachusetts Institute of Technology Published in Scandinavian Journal of Statistics,
More informationA stochastic particle system for the Burgers equation.
A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u
More informationExponential martingales: uniform integrability results and applications to point processes
Exponential martingales: uniform integrability results and applications to point processes Alexander Sokol Department of Mathematical Sciences, University of Copenhagen 26 September, 2012 1 / 39 Agenda
More informationA diamagnetic inequality for semigroup differences
A diamagnetic inequality for semigroup differences (Irvine, November 10 11, 2001) Barry Simon and 100DM 1 The integrated density of states (IDS) Schrödinger operator: H := H(V ) := 1 2 + V ω =: H(0, V
More informationJae Gil Choi and Young Seo Park
Kangweon-Kyungki Math. Jour. 11 (23), No. 1, pp. 17 3 TRANSLATION THEOREM ON FUNCTION SPACE Jae Gil Choi and Young Seo Park Abstract. In this paper, we use a generalized Brownian motion process to define
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationAsymptotic results for empirical measures of weighted sums of independent random variables
Asymptotic results for empirical measures of weighted sums of independent random variables B. Bercu and W. Bryc University Bordeaux 1, France Seminario di Probabilità e Statistica Matematica Sapienza Università
More informationNew ideas in the non-equilibrium statistical physics and the micro approach to transportation flows
New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,
More informationA PRACTICAL WAY FOR ESTIMATING TAIL DEPENDENCE FUNCTIONS
Statistica Sinica 20 2010, 365-378 A PRACTICAL WAY FOR ESTIMATING TAIL DEPENDENCE FUNCTIONS Liang Peng Georgia Institute of Technology Abstract: Estimating tail dependence functions is important for applications
More informationSPATIAL STRUCTURE IN LOW DIMENSIONS FOR DIFFUSION LIMITED TWO-PARTICLE REACTIONS
The Annals of Applied Probability 2001, Vol. 11, No. 1, 121 181 SPATIAL STRUCTURE IN LOW DIMENSIONS FOR DIFFUSION LIMITED TWO-PARTICLE REACTIONS By Maury Bramson 1 and Joel L. Lebowitz 2 University of
More informationGiovanni Migliorati. MATHICSE-CSQI, École Polytechnique Fédérale de Lausanne
Analysis of the stability and accuracy of multivariate polynomial approximation by discrete least squares with evaluations in random or low-discrepancy point sets Giovanni Migliorati MATHICSE-CSQI, École
More informationDirection: This test is worth 250 points and each problem worth points. DO ANY SIX
Term Test 3 December 5, 2003 Name Math 52 Student Number Direction: This test is worth 250 points and each problem worth 4 points DO ANY SIX PROBLEMS You are required to complete this test within 50 minutes
More informationFrom now on, we will represent a metric space with (X, d). Here are some examples: i=1 (x i y i ) p ) 1 p, p 1.
Chapter 1 Metric spaces 1.1 Metric and convergence We will begin with some basic concepts. Definition 1.1. (Metric space) Metric space is a set X, with a metric satisfying: 1. d(x, y) 0, d(x, y) = 0 x
More informationPARAMETER ESTIMATION FOR SPDES WITH MULTIPLICATIVE FRACTIONAL NOISE
PARAMETER ESTIMATION FOR SPDES WITH MULTIPLICATIVE FRACTIONAL NOISE IGOR CIALENCO Abstract. We study parameter estimation problem for diagonalizable stochastic partial differential equations driven by
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 informationTHE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL. Haseo Ki and Young One Kim
THE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL Haseo Ki Young One Kim Abstract. The zero-distribution of the Fourier integral Q(u)eP (u)+izu du, where P is a polynomial with leading
More informationAsymptotical distribution free test for parameter change in a diffusion model (joint work with Y. Nishiyama) Ilia Negri
Asymptotical distribution free test for parameter change in a diffusion model (joint work with Y. Nishiyama) Ilia Negri University of Bergamo (Italy) ilia.negri@unibg.it SAPS VIII, Le Mans 21-24 March,
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued
Introduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and
More informationOn a class of stochastic differential equations in a financial network model
1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University
More informationLAN property for ergodic jump-diffusion processes with discrete observations
LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationLower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions
Lower Tail Probabilities and Normal Comparison Inequalities In Memory of Wenbo V. Li s Contributions Qi-Man Shao The Chinese University of Hong Kong Lower Tail Probabilities and Normal Comparison Inequalities
More informationApproximating the gains from trade in large markets
Approximating the gains from trade in large markets Ellen Muir Supervisors: Peter Taylor & Simon Loertscher Joint work with Kostya Borovkov The University of Melbourne 7th April 2015 Outline Setup market
More informationVISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5.
VISCOSITY SOLUTIONS PETER HINTZ We follow Han and Lin, Elliptic Partial Differential Equations, 5. 1. Motivation Throughout, we will assume that Ω R n is a bounded and connected domain and that a ij C(Ω)
More informationBayesian Estimations for Diagonalizable Bilinear SPDEs
ayesian Estimations for Diagonalizable ilinear SPDEs Ziteng Cheng zcheng7@haw.iit.edu Igor Cialenco cialenco@iit.edu http://math.iit.edu/~igor uoting Gong rgong2@iit.edu http://mypages.iit.edu/~rgong2
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationAnnealed Brownian motion in a heavy tailed Poissonian potential
Annealed Brownian motion in a heavy tailed Poissonian potential Ryoki Fukushima Research Institute of Mathematical Sciences Stochastic Analysis and Applications, Okayama University, September 26, 2012
More informationProperties of an infinite dimensional EDS system : the Muller s ratchet
Properties of an infinite dimensional EDS system : the Muller s ratchet LATP June 5, 2011 A ratchet source : wikipedia Plan 1 Introduction : The model of Haigh 2 3 Hypothesis (Biological) : The population
More informationLecture 8: Information Theory and Statistics
Lecture 8: Information Theory and Statistics Part II: Hypothesis Testing and Estimation I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 22, 2015
More informationWiener Measure and Brownian Motion
Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca October 22nd, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationSequential Monte Carlo Samplers for Applications in High Dimensions
Sequential Monte Carlo Samplers for Applications in High Dimensions Alexandros Beskos National University of Singapore KAUST, 26th February 2014 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Alex
More informationFractional Cahn-Hilliard equation
Fractional Cahn-Hilliard equation Goro Akagi Mathematical Institute, Tohoku University Abstract In this note, we review a recent joint work with G. Schimperna and A. Segatti (University of Pavia, Italy)
More informationDiophantine approximation of beta expansion in parameter space
Diophantine approximation of beta expansion in parameter space Jun Wu Huazhong University of Science and Technology Advances on Fractals and Related Fields, France 19-25, September 2015 Outline 1 Background
More informationThe dynamical Sine-Gordon equation
The dynamical Sine-Gordon equation Hao Shen (University of Warwick) Joint work with Martin Hairer October 9, 2014 Dynamical Sine-Gordon Equation Space dimension d = 2. Equation depends on parameter >0.
More informationOn the quantiles of the Brownian motion and their hitting times.
On the quantiles of the Brownian motion and their hitting times. Angelos Dassios London School of Economics May 23 Abstract The distribution of the α-quantile of a Brownian motion on an interval [, t]
More informationON SOME CONSTANTS FOR OSCILLATION AND STABILITY OF DELAY EQUATIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 139, Number 11, November 2011, Pages 4017 4026 S 0002-9939(2011)10820-7 Article electronically published on March 28, 2011 ON SOME CONSTANTS FOR
More informationTesting Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata
Maura Department of Economics and Finance Università Tor Vergata Hypothesis Testing Outline It is a mistake to confound strangeness with mystery Sherlock Holmes A Study in Scarlet Outline 1 The Power Function
More informationIn Memory of Wenbo V Li s Contributions
In Memory of Wenbo V Li s Contributions Qi-Man Shao The Chinese University of Hong Kong qmshao@cuhk.edu.hk The research is partially supported by Hong Kong RGC GRF 403513 Outline Lower tail probabilities
More informationNodal Sets of High-Energy Arithmetic Random Waves
Nodal Sets of High-Energy Arithmetic Random Waves Maurizia Rossi UR Mathématiques, Université du Luxembourg Probabilistic Methods in Spectral Geometry and PDE Montréal August 22-26, 2016 M. Rossi (Université
More informationDetection Theory. Composite tests
Composite tests Chapter 5: Correction Thu I claimed that the above, which is the most general case, was captured by the below Thu Chapter 5: Correction Thu I claimed that the above, which is the most general
More informationA LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION
A LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION O. SAVIN. Introduction In this paper we study the geometry of the sections for solutions to the Monge- Ampere equation det D 2 u = f, u
More informationMagnetic wells in dimension three
Magnetic wells in dimension three Yuri A. Kordyukov joint with Bernard Helffer & Nicolas Raymond & San Vũ Ngọc Magnetic Fields and Semiclassical Analysis Rennes, May 21, 2015 Yuri A. Kordyukov (Ufa) Magnetic
More informationIf there exists a threshold k 0 such that. then we can take k = k 0 γ =0 and achieve a test of size α. c 2004 by Mark R. Bell,
Recall The Neyman-Pearson Lemma Neyman-Pearson Lemma: Let Θ = {θ 0, θ }, and let F θ0 (x) be the cdf of the random vector X under hypothesis and F θ (x) be its cdf under hypothesis. Assume that the cdfs
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationLECTURE 2: LOCAL TIME FOR BROWNIAN MOTION
LECTURE 2: LOCAL TIME FOR BROWNIAN MOTION We will define local time for one-dimensional Brownian motion, and deduce some of its properties. We will then use the generalized Ray-Knight theorem proved in
More informationSTAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots. March 8, 2015
STAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots March 8, 2015 The duality between CI and hypothesis testing The duality between CI and hypothesis
More informationChapter 7. Hypothesis Testing
Chapter 7. Hypothesis Testing Joonpyo Kim June 24, 2017 Joonpyo Kim Ch7 June 24, 2017 1 / 63 Basic Concepts of Testing Suppose that our interest centers on a random variable X which has density function
More informationParacontrolled KPZ equation
Paracontrolled KPZ equation Nicolas Perkowski Humboldt Universität zu Berlin November 6th, 2015 Eighth Workshop on RDS Bielefeld Joint work with Massimiliano Gubinelli Nicolas Perkowski Paracontrolled
More informationLecture 17 Brownian motion as a Markov process
Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is
More informationChapter 7. Basic Probability Theory
Chapter 7. Basic Probability Theory I-Liang Chern October 20, 2016 1 / 49 What s kind of matrices satisfying RIP Random matrices with iid Gaussian entries iid Bernoulli entries (+/ 1) iid subgaussian entries
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationPositive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network. Haifa Statistics Seminar May 5, 2008
Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network Yoni Nazarathy Gideon Weiss Haifa Statistics Seminar May 5, 2008 1 Outline 1 Preview of Results 2 Introduction Queueing
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationBEST TESTS. Abstract. We will discuss the Neymann-Pearson theorem and certain best test where the power function is optimized.
BEST TESTS Abstract. We will discuss the Neymann-Pearson theorem and certain best test where the power function is optimized. 1. Most powerful test Let {f θ } θ Θ be a family of pdfs. We will consider
More informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More informationDiffusion of wave packets in a Markov random potential
Diffusion of wave packets in a Markov random potential Jeffrey Schenker Michigan State University Joint work with Yang Kang (JSP 134 (2009)) WIP with Eman Hamza & Yang Kang ETH Zürich 9 June 2009 Jeffrey
More informationModel Selection and Geometry
Model Selection and Geometry Pascal Massart Université Paris-Sud, Orsay Leipzig, February Purpose of the talk! Concentration of measure plays a fundamental role in the theory of model selection! Model
More informationBranching Brownian motion seen from the tip
Branching Brownian motion seen from the tip J. Berestycki 1 1 Laboratoire de Probabilité et Modèles Aléatoires, UPMC, Paris 09/02/2011 Joint work with Elie Aidekon, Eric Brunet and Zhan Shi J Berestycki
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationLecture 7 Introduction to Statistical Decision Theory
Lecture 7 Introduction to Statistical Decision Theory I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 20, 2016 1 / 55 I-Hsiang Wang IT Lecture 7
More informationWiener Chaos Solution of Stochastic Evolution Equations
Wiener Chaos Solution of Stochastic Evolution Equations Sergey V. Lototsky Department of Mathematics University of Southern California August 2003 Joint work with Boris Rozovskii The Wiener Chaos (Ω, F,
More informationChapter 9: Hypothesis Testing Sections
Chapter 9: Hypothesis Testing Sections 9.1 Problems of Testing Hypotheses 9.2 Testing Simple Hypotheses 9.3 Uniformly Most Powerful Tests Skip: 9.4 Two-Sided Alternatives 9.6 Comparing the Means of Two
More informationSelf-normalized Cramér-Type Large Deviations for Independent Random Variables
Self-normalized Cramér-Type Large Deviations for Independent Random Variables Qi-Man Shao National University of Singapore and University of Oregon qmshao@darkwing.uoregon.edu 1. Introduction Let X, X
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 informationLocal strong convexity and local Lipschitz continuity of the gradient of convex functions
Local strong convexity and local Lipschitz continuity of the gradient of convex functions R. Goebel and R.T. Rockafellar May 23, 2007 Abstract. Given a pair of convex conjugate functions f and f, we investigate
More informationVariational approach to mean field games with density constraints
1 / 18 Variational approach to mean field games with density constraints Alpár Richárd Mészáros LMO, Université Paris-Sud (based on ongoing joint works with F. Santambrogio, P. Cardaliaguet and F. J. Silva)
More informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More informationLecture 21 Representations of Martingales
Lecture 21: Representations of Martingales 1 of 11 Course: Theory of Probability II Term: Spring 215 Instructor: Gordan Zitkovic Lecture 21 Representations of Martingales Right-continuous inverses Let
More informationLecture 21. Hypothesis Testing II
Lecture 21. Hypothesis Testing II December 7, 2011 In the previous lecture, we dened a few key concepts of hypothesis testing and introduced the framework for parametric hypothesis testing. In the parametric
More informationSurveillance of BiometricsAssumptions
Surveillance of BiometricsAssumptions in Insured Populations Journée des Chaires, ILB 2017 N. El Karoui, S. Loisel, Y. Sahli UPMC-Paris 6/LPMA/ISFA-Lyon 1 with the financial support of ANR LoLitA, and
More informationAn Introduction to Malliavin calculus and its applications
An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart
More informationSome Aspects of Universal Portfolio
1 Some Aspects of Universal Portfolio Tomoyuki Ichiba (UC Santa Barbara) joint work with Marcel Brod (ETH Zurich) Conference on Stochastic Asymptotics & Applications Sixth Western Conference on Mathematical
More information