SUPPLEMENT TO POSTERIOR CONTRACTION AND CREDIBLE SETS FOR FILAMENTS OF REGRESSION FUNCTIONS. North Carolina State University

Similar documents
Piecewise Smooth Solutions to the Burgers-Hilbert Equation

DISCUSSION: COVERAGE OF BAYESIAN CREDIBLE SETS. By Subhashis Ghosal North Carolina State University

Introduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued

Metric Spaces and Topology

Measure and Integration: Solutions of CW2

The Existence of the Riemann Integral

Harmonic Analysis Homework 5

Math 328 Course Notes

Approximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle

REAL ANALYSIS II: PROBLEM SET 2

Refining the Central Limit Theorem Approximation via Extreme Value Theory

Solution of the 7 th Homework

NOTES ON EXISTENCE AND UNIQUENESS THEOREMS FOR ODES

Solutions to Problem Set 5 for , Fall 2007

FIRST YEAR CALCULUS W W L CHEN

SUPPLEMENT TO EMPIRICAL BAYES ORACLE UNCERTAINTY QUANTIFICATION FOR REGRESSION. VU Amsterdam and North Carolina State University

Economics 204 Fall 2011 Problem Set 2 Suggested Solutions

Concentration behavior of the penalized least squares estimator

Chapter 1. Optimality Conditions: Unconstrained Optimization. 1.1 Differentiable Problems

Properties of the Scattering Transform on the Real Line

Analysis Finite and Infinite Sets The Real Numbers The Cantor Set

CONTINUITY OF A QUEUEING INTEGRAL REPRESENTATION IN THE M 1 TOPOLOGY. Department of Industrial Engineering and Operations Research Columbia University

Analysis III. Exam 1

1 Preliminaries on locally compact spaces and Radon measure

Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011

The deterministic Lasso

FUNDAMENTALS OF REAL ANALYSIS by. IV.1. Differentiation of Monotonic Functions

for 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

Mathematical Analysis Outline. William G. Faris

d(x n, x) d(x n, x nk ) + d(x nk, x) where we chose any fixed k > N

Bayesian Regularization

A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem

Nonparametric Function Estimation with Infinite-Order Kernels

Proof. We indicate by α, β (finite or not) the end-points of I and call

(Part 1) High-dimensional statistics May / 41

NETWORK-REGULARIZED HIGH-DIMENSIONAL COX REGRESSION FOR ANALYSIS OF GENOMIC DATA

HW #1 SOLUTIONS. g(x) sin 1 x

HOMEWORK ASSIGNMENT 6

Existence and Comparisons for BSDEs in general spaces

Introduction to Empirical Processes and Semiparametric Inference Lecture 13: Entropy Calculations

Fixed Point Theorem and Sequences in One or Two Dimensions

Robust Estimation and Inference for Extremal Dependence in Time Series. Appendix C: Omitted Proofs and Supporting Lemmata

M.S.N. Murty* and G. Suresh Kumar**

Appendix. Online supplement to Coordinated logistics with a truck and a drone

Examination paper for TMA4145 Linear Methods

Analysis Comprehensive Exam Questions Fall 2008

Learning Theory. Ingo Steinwart University of Stuttgart. September 4, 2013

Supplement to Quantile-Based Nonparametric Inference for First-Price Auctions

Homework If the inverse T 1 of a closed linear operator exists, show that T 1 is a closed linear operator.

Functional Analysis (2006) Homework assignment 2

Index Models for Sparsely Sampled Functional Data. Supplementary Material.

4 Sobolev spaces, trace theorem and normal derivative

Take-Home Final Examination of Math 55a (January 17 to January 23, 2004)

Oberwolfach Preprints

Existence Results for Multivalued Semilinear Functional Differential Equations

Putzer s Algorithm. Norman Lebovitz. September 8, 2016

ORTHOGONAL SERIES REGRESSION ESTIMATORS FOR AN IRREGULARLY SPACED DESIGN

CS229T/STATS231: Statistical Learning Theory. Lecturer: Tengyu Ma Lecture 11 Scribe: Jongho Kim, Jamie Kang October 29th, 2018

First Variation of a Functional

A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES

ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM

Lecture 3: September 10

Functional Analysis F3/F4/NVP (2005) Homework assignment 3

MA2223 Tutorial solutions Part 1. Metric spaces

4th Preparation Sheet - Solutions

Herz (cf. [H], and also [BS]) proved that the reverse inequality is also true, that is,

Introduction to the Numerical Solution of IVP for ODE

l(y j ) = 0 for all y j (1)

Bootcamp. Christoph Thiele. Summer 2012

Richard S. Palais Department of Mathematics Brandeis University Waltham, MA The Magic of Iteration

MATH 51H Section 4. October 16, Recall what it means for a function between metric spaces to be continuous:

Theoretical Statistics. Lecture 14.

OPTIMAL POINTWISE ADAPTIVE METHODS IN NONPARAMETRIC ESTIMATION 1

Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.

Multivariate Least Weighted Squares (MLWS)

Consistency and Convergence

Normal approximation of Poisson functionals in Kolmogorov distance

Euler Equations: local existence

A Lower Bound Theorem. Lin Hu.

Nonlinear Control Systems

Math 201 Handout 1. Rich Schwartz. September 6, 2006

An introduction to Birkhoff normal form

Numerical Solutions to Partial Differential Equations

Bootcamp. Christoph Thiele. Summer As in the case of separability we have the following two observations: Lemma 1 Finite sets are compact.

EXISTENCE OF SOLUTIONS FOR A RESONANT PROBLEM UNDER LANDESMAN-LAZER CONDITIONS

Houston Journal of Mathematics. c 2015 University of Houston Volume 41, No. 4, 2015

A MATROID EXTENSION RESULT

The Contraction Mapping Theorem and the Implicit Function Theorem

UPPER AND LOWER HENSTOCK INTEGRALS

Institut für Mathematik

at time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))

arxiv: v1 [math.oc] 22 Sep 2016

S chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.

Generalization theory

Atomic decompositions of square-integrable functions

AN EFFECTIVE METRIC ON C(H, K) WITH NORMAL STRUCTURE. Mona Nabiei (Received 23 June, 2015)

2.2 Some Consequences of the Completeness Axiom

Existence and Uniqueness

On variable bandwidth kernel density estimation

On the Intrinsic Differentiability Theorem of Gromov-Schoen

Transcription:

Submitted to the Annals of Statistics SUPPLEMENT TO POSTERIOR CONTRACTION AND CREDIBLE SETS FOR FILAMENTS OF REGRESSION FUNCTIONS By Wei Li and Subhashis Ghosal North Carolina State University The supplementary file contains proofs for Lemma 5.7, 8., and 8.4. Proof of Lemma 5.7. Let be an arbitrary point on ˆL. Thus Υ = ˆL and Υ t L for some t >. Note that inf y L y Υ Υ t t as V =. Let D Υ,V ft := d dt fυ t = fυ t T V Υ t and D 2 Υ,V ft := fυ t, V Υ t, V Υ t = fυ t T V Υ tv Υ t + V Υ t T HfΥ tv Υ t, where the second line is due to the chain rule. A Taylor epansion yields that D Υ,V ft D Υ,V ft = t t D 2 Υ,V f t for some t between and t. In particular, since D Υ,V ft =, letting t =, we obtain D Υ,V f = t D 2 Υ,V f t. Furthermore, D Υ,V f = D Υ,V f D Υ, ˆV ˆf = fυ T V Υ ˆfΥ T ˆV Υ sup f T V ˆf T ˆV sup Csup f T V ˆV + sup f ˆf T ˆV V ˆV + sup f ˆf C sup V ˆV. By the uniform continuity of f, V, V and Hf and the continuity of Υ t in t, without loss of generality, we can make t small enough imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

2 LI, W. AND GHOSAL, S. see comments after the proof. Thus we have DΥ 2,V f t DΥ 2,V ft < η/2. By Assumption A3, DΥ 2,V ft > η, and hence DΥ 2,V f t DΥ 2,V ft DΥ 2,V f t DΥ 2,V ft > η/2. Therefore, inf y L y t C η sup V ˆV. Thus d ˆL L C η sup V ˆV. Similarly, dl ˆL C η sup V ˆV. Therefore, HausL, ˆL C η sup [,] 2 V ˆV. Recall that V = Gdf 2. Now since G is a fied Lipschitz continuous function, it is easy to get the upper bound for sup [,] 2 V ˆV in terms of the supremum distance of the derivatives of f ˆf. In the proof, t can be made arbitrarily small in the limit. To see this, if Assumption A5 holds for the f or more precisely, Υ, then Υ t ˆΥ ˆt = Υ t = Υ t Υ > C G t. Since ˆΥ ˆt Υ t can be made arbitrarily small due to the closedness in supremum norm see previous theorems, t can be made arbitrarily small. Proof of Lemma 8.. These results can be directly adapted from [] by noting the highest degree of derivatives in each epression. We shall show the rates for sup Hf Hf F, sup d 2 f d 2 f F and sup V V F. Notice that Hf Hf F = f 2, f 2, 2 + 2f, f, 2 + f,2 f,2 2 /2, 4 /2 f 2, f 2, + 2 f, f, + f,2 f,2. The rate for sup Hf Hf F then follows easily. Also, the contraction rates for d 2 f d 2 f F follows from the inequality d 2 f d 2 f F 6 /2 r: r =3 Dr f D r f. imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

Since V V = Gd 2 f d 2 f Gd 2 f d 2 f, which is 2 by 2 matri, straightforward calculation gives its, element G,, d 2 ff 3, G,, d 2 f f 3,, + G,, d 2 ff 2, G,, d 2 f f 2,, + G,, d 2 ff,2 G,, d 2 f f,2. The absolute value of the first summand is bounded by the sum of G,, d 2 ff 3, G,, d 2 ff 3,, 3 and G,, d 2 ff 3, G,, d 2 f f 3,. Noting that d 2 f contracts to d 2 f uniformly in, hence {d 2 f : [, ] 2 } Q δ with posterior probability tending to, the first term is bounded by a constant multiple of f 3, f 3,, in view of the result 4 of Remark 8.. By the same remark and assumption f α, <, the second term is bounded by d 2 f d 2 f f 3,. Using similar arguments for the second and third summand, one can see that the absolute value of the, element of V V is bounded by f 3, f 3, + f 2, f 2, + f,2 f,2. Dealing the rest elements of V V similarly, we can see that V V F r: r =3 Dr f D r f. Proof of Lemma 8.4. We shall consider only two separate cases i r = 2, and ii r =,. i. For r = 2, similarly for r =, 2, a, t = t t J t j G sb j Υ, sb j2 Υ 2, sds j Υ, s B j2 Υ 2, sds j Υ, s ds, imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

4 LI, W. AND GHOSAL, S. the last line follows by J 2 j 2 = B j 2 2 = and noticing that B j,q = q q 2B j,q 2 t j t j q +2t j t j q + + q q 2B j,q 2 t j t j q +t j t j q + q q 2B j,q 2 t j t j q +2t j t j q + q q 2B j 2,q 2 t j 2 t j q +t j t j q, which implies J j = B j Υ, s 4J 2. Therefore, a, t J 2. Net, let S j = [t,j q 2, t,j ], S j2 = [t 2,j2 q 2, t 2,j2 ] and j,j 2 s := {s : Υ s S j S j2 }. Turn to a, t 2, which is = = t G s G s 2 G sb j Υ, sb j2 Υ 2, sds B B j Υ, s B j2 Υ 2, s dsds j Υ, sb j2 Υ 2, s j Υ, s B j2 Υ 2, s B j Υ, s B j2 Υ 2, s dsds j,j 2 s j,j 2 s B j Υ, s B j2 Υ 2, s j Υ, s B j2 Υ 2, s dsds. The last equality is obtained as follows. Since for any fi j, j 2, B j is supported on S j and B j2 is supported on S j2. So j,j 2 s j,j 2 s = {s, s : Υ s S j S j2 and Υ s S j S j2 }. Notice that for n large, S j S j2 J and Υ s S j S j2 and Υ s S j S j2 implies that Υ s Υ s J. By Assumption A5, C G s s s s Υ s Υ s /s s J, and hence s s J. Therefore, for n large enough, above quantity can be further bounded by a constant multiple of { s s < CJ } B j Υ, s B j2 Υ 2, s j Υ, s B j2 Υ 2, s dsds. imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

Noting B j2 2 and J 2 j 2 = B j 2 =, it can be further bounded by { s s < CJ } B j Υ, s B j2 Υ 2, s B j Υ, s dsds J J 2 J t { s s < CJ } j j = j Υ, s ds. j Υ, sb j Υ, s dsds From argument used in bounding a, t, we have J j = B j Υ, s J 2. This completes the proof for a, t 2 J 3. For the third result, we write where a, t, t := and a 2 t,, := a, t a, t 2 a, t, t 2 + a 2 t,, 2, t First, note that a, t, t 2 = GΥ sb 2, J,J 2 Υ sds GΥ sb 2, J,J 2 Υ sds t J 2 J 4J 2 ds t = J 3 t t, t 5 GΥ sb 2, J,J 2 Υ sds, GΥ sb 2, J,J 2 Υ sds. 2 GΥ sb j Υ, sb j2 Υ 2, sds where the second line follows by a similar argument used to bound a, t 2. Net, a 2 t,, 2 is given by [ GΥ sb j Υ, sb j2 Υ 2, s GΥ sb j Υ, sb j2 Υ 2, sds] 2, imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

6 LI, W. AND GHOSAL, S. which is bounded up to a multiple constant by + + 2 GΥ s B j Υ, s B j Υ, s B j2 Υ 2, sds 2 GΥ s B j2 Υ 2, s B j2 Υ 2, s B j Υ, sds GΥ s GΥ s B j Υ, sb j2 Υ 2, sds 2. Bound the first term in the right hand side of above epression as B j2 Υ 2, sds 2 2 GΥ s B j Υ, s B j Υ, s B j2 Υ 2, sds GΥ s j Υ, s Υ, s Υ, s 2 2 J 2 2 J 5. GΥ 2 s B j Υ, s B j2 Υ 2, sds j Υ, s ds The second term is bounded as 2 2 GΥ s B j2 Υ 2, s B j2 Υ 2, s B j Υ, sds GΥ 2 s Υ 2, s Υ 2, s B j 2 Υ 2, s B j Υ, s GΥ 2 s B j 2 Υ 2, s B j Υ, s imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

2 { s s CJ } 2 B j 2 Υ 2, s j Υ, s B j 2 Υ 2, s j Υ, s dsds J 5 2, { s s CJ }J 2 j = 7 j Υ, s j Υ, s dsds where the second line follows from the mean value theorem, the third line from the Lipschitz continuity of Υ in Remark 8. and fourth line by a similar argument used to bound a, t 2. The third term GΥ s GΥ s 2 2 J 3, B j Υ, sb j2 Υ 2, sds Υ s Υ s B j Υ, sb j2 Υ 2, sds 2 B j Υ, s B j2 Υ 2, sds where the second line holds by mean value theorem and the third line holds due to Lipschitz continuity of Υ in Remark 8. and last line holds by similar argument for a, t 2. In summary, we have that a 2 t,, 2 J 5 2 and a, t, t 2 J 3 t t. ii. Now turning to the case r =,. By similar argument we have a, t = t t J t j = G sb j Υ, sb j 2 Υ 2, sds B j Υ, s B j 2 Υ 2, s ds B j Υ, s ds J 2 2 2 imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

8 LI, W. AND GHOSAL, S. Likewise, a, t 2 = J J2 t G sb j Υ, sb j 2 Υ 2, sds can be bounded by J { s s < CJ } B j Υ, sb j 2 Υ 2, s B j Υ, s B j 2 Υ 2, s dsds B j Υ, s dsds J 2 t J 3 t t t J 3 J t j = { s s < CJ } B j Υ, s B j 2 Υ 2, s { s s < CJ } B j Υ, sb j Υ, s dsds { s s < CJ } B j Υ, s dsds j = j B j Υ, s ds which is bounded by a constant multiple of J 3. The third result a, t, t 2 J 3 t t and a 2 t,, 2 J 5 2 can be derived in a similar manner, since a, t, t 2 = t J 2 J 4J 2 ds t t 2 GΥ sb j Υ, sb j 2 Υ 2, sds which equals to J 3 t t ; where the second line follows by a similar argument used in bounding a, t 2. Net, to bound a 2 t,, 2, we need to estimate the following three terms 2 2, GΥ s B j Υ, s B j Υ, s B j 2 Υ 2, sds imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28

2, GΥ s B j 2 Υ 2, s B j 2 Υ 2, s B j Υ, sds GΥ 2. s GΥ s B j Υ, sb j 2 Υ 2, sds This can be done by similar argument used for the case r = 2, and we omit the details. References. [] Yoo, W. W. and Ghosal, S. 26. Supremum norm posterior contraction and credible sets for nonparametric multivariate regression. The Annals of Statistics 44 69 2. 9 Department of Statistics North Carolina State University Raleigh, North Carolina 27695-823 USA E-mail: wli35@ncsu.edu Department of Statistics North Carolina State University Raleigh, North Carolina 27695-823 USA E-mail: sghosal@stat.ncsu.edu imsart-aos ver. 24//6 file: Filaments_AOS_supp_9Mar28.te date: March 9, 28