arxiv: v1 [math.pr] 9 Feb 2011

Transcription

1 Multivariate piecewise linear interpolation of a random field Konrad Abramowicz, Oleg Seleznev, Department of Mathematics and Mathematical Statistics Umeå University, SE Umeå, Sweden arxiv: v1 [math.pr] 9 Feb 011 February 10, 011 Abstract We consider a multivariate piecewise linear interpolation of a continuous random field on a d-dimensional cube. The approximation performance is measured by the integrated mean square error. Multivariate piecewise linear interpolator is defined by N field observations on a locations gridor design. We investigate the class of locally stationary random fields whose local behavior is like a fractional Brownian field in mean square sense and find the asymptotic approximation accuracy for a sequence of designs for large N. Moreover, for certain classes of continuous and continuously differentiable fields we provide the upper bound for the approximation accuracy in the uniform mean square norm. Keywords: approximation, random field, sampling design, multivariate piecewise linear interpolator 1 Introduction LetarandomfieldXt, t [0,1] d, withfinitesecond momentbeobservedat finitenumberof points. Suppose further that the points are vertices of hyperrectangles generated by a grid in a unit hypercube. At any unsampled point we approximate the value of the field by a piecewise linear multivariate interpolator, which is a natural extension of a conventional one-dimensional piecewise linear interpolator. The approximation accuracy is measured by the integrated mean squared error. This paper aims modelling random fields with given accuracy based on a finite number of observations. Following Berman 1974, we extend the concept of local stationarity for random fields and focus on fields satisfying this condition. For quadratic mean q.m. continuous locally stationary random fields, we derive the exact asymptotic behavior of the approximation error. A method is proposed for determining the asymptotically optimal knot sample points distribution between the mesh dimensions. We also study optimality of knot allocation along coordinates of the sampling grid. Additionally, for q.m. continuous and continuously differentiable fields satisfying Hölder type conditions, we determine asymptotical upper bounds for the approximation accuracy. The problem of random field approximation arises in many research and applied areas, like Gaussian random fields modellingadler and Taylor, 007; Brouste et al., 007, environmental and geoscienceschristakos, 199; Stein, 1999, sensor networks Zhang and Wicker, 005, and image processing Pratt, 007. The upper bound for the approximation error for isotropic random fields satisfying Hölder type conditions is given in Ritter et al Müller-Gronbach 1998 consider affine linear approximation methods and hyperbolic cross designs for fields with covariance function of tensor type. An optimal allocation of the observations for Gaussian random fields with product type kernel is investigated in Müller-Gronbach and Schwabe Su 1997 studies limit behavior of the piecewise constant estimator for random fields with a particular form of covariance function. Benhenni 001 investigates exact asymptotics of stationary spatial process approximation based on an equidistant sampling. The approximation complexity and the curse of dimensionality 1

2 for additive random fields are broadly discussed in Lifshits and Zani 008. In one-dimensional case, the piecewise linear interpolation of continuous stochastic processes is considered in, e.g., Seleznev Results for approximation of locally stationary processes can be found in, e.g., Seleznev 000; Hüsler et al. 003; Abramowicz and Seleznev 011. Ritter 000 contains a very detailed survey of various random process and field approximation problems. For an extensive studies of approximation problems in deterministic setting, we refer to, e.g., Nikolskii 1975; de Boor et al. 008; Kuo et al The paper is organized as follows. First we introduce a basic notation. In Section, we consider a piecewise multivariate linear approximation of continuous fields which local behavior is like a fractional Brownian field in mean square sense. We derive exact asymptotics and a formula for the optimal interdimensional knot distribution. In the second part of this section, we provide an asymptotical upper bound for the approximation accuracy for q.m. continuous and differentiable fields satisfying Hölder type conditions. In Section 3, we present the results of numerical experiments, while Section 4 contains the proofs of the statements from Section. 1.1 Basic notation Let X = Xt,t D := [0,1] d, be a random field defined on a probability space Ω,F,P. Assume that for every t, the random variable Xt lies in the normed linear space L Ω = L Ω,F,P of random variables with finite second moment and identified equivalent elements with respect to P. We set ξ := Eξ 1/ for all ξ L Ω and consider the approximation based on the normed linear spaces of q.m. continuous and continuously differentiable random fields denoted by CD and C 1 D, respectively. We define the norm for any X CD by setting X p := D Xt p dt 1/p, 1 p <, and X := max t D Xt. For p =, we call the norm integrated mean squared norm and the corresponding measure of approximation accuracy the integrated mean squared error IMSE. Now we introduce the classes of random fields used throughout this paper. For k d, let l = l 1,...,l k be a vector of positive integers such that k =1 l = d, and let L i := i =1 l,i = 0,...k, L 0 = 0, be the sequence of its cumulative sums. Then the vector l defines the l-decomposition of D into D 1 D...D k, with the l -cube D = [0,1] l, = 1,...,k. For any s D, we denote the coordinates vector corresponding to the -th component of the decomposition by s, i.e., s = s l := s L 1 +1,...,s L D, = 1,...,k. For a vector α = α 1,...,α k, 0 < α <, = 1,...,k, and the decomposition vector l = l 1,...,l k, we define s α := with the Euclidean norms s, = 1,...,k. s α =1 for all s D For a random field X C[0,1] d, we say that i X Cl α[0,1]d,c if for some α, l, and a positive constant C, the random field X satisfies the Hölder condition, i.e., Xt+s Xt C s α for all t,t+s [0,1] d, 1

3 ii X B α l [0,1]d,c if for some α, l, and a vector function ct = c 1 t,...,c k t, t [0,1] d, the random field X is locally stationary, i.e., Xt+s Xt k =1 c kt s α 1 as s 0 uniformly in t [0,1] d, with positive and continuous functions c 1,...,c k. We assume additionally that for = 1,...,k, the function c is invariant with respect to coordinates permutation within the -th component. For the classes C α l and B α l, the withincomponent smoothness is defined by the vector α = α 1,...,α k. We denote the vector describing the smoothness for each coordinate by α = α 1,...,α d, where α i = α, i = L 1 +1,...,L, = 1,...,k. Example 1. Let m = m 1,...,m k be a decomposition vector of [0,1] m, and m = k =1 m. Denote by B β,m t, t [0,1] m, β = β 1,...,β k, 0 < β <, = 1,...,k, an m-dimensional fractional Brownian field with covariance function rt,s = 1 t β + s β t s β. Then B β,m has stationary increments, B β,m t+s B β,m t = s β, t,t+s [0,1] m, and therefore, B β,m B β md,c with local stationarity functions c 1 t =... = c k t = 1, t [0,1] m. In particular, if k = 1, then B β,m t, t [0,1] m, 0 < β <, m Æ, is an m-dimensional fractal Brownian field with covariance function rt,s = 1 t β + s β t s β, t,t+s [0,1] m. 3 For X C 1 [0,1] d, we write X t, t [0,1]d, to denote a q.m. partial derivative of X with respect to the -th coordinate, and say that X C 1,α [0,1] d,c if there exist a vector α = α 1,...,α d and a positive constant C such that each partial derivative X is Hölder continuous with respect to the -th coordinate, i.e., if for all t,t+s [0,1] d, X t 1,...,t +s,...,t d X t 1,...,t,...,t d C s α = 1,...,d. 4 Moreover, we say that X C 1,α l [0,1] d,c with α = α 1,...,α k if X C 1,α [0,1] d,c and for a given partition vector l, α i := α L i 1 +1 =... = α L i, i = 1,...,k. Let X be sampled at N distinct design points T N. We consider cross regular sequences of sampling designs T N := {t i = t 1,i1,...,t d,id : i = i 1,...,i d, 0 i k n k, k = 1,...,d} defined by the one-dimensional grids t,i 0 h vdv = i n, i = 0,1,...,n, = 1,...,d, where h s, s [0,1], = 1,...,d, are positive and continuous density functions, say, withindimensional densities, and let h t := h 1 t 1,...,h d t d. The interdimensional knot distribution is determined by a vector function π : Æ Æ d : π N := n 1 N,...,n d N, 3

4 where lim N n N =, = 1,...,d, and the condition d n N+1 = N =1 is satisfied. We suppress the argument N for the sampling grid sizes n = n N, = 1,...,d, when doing so causes no confusion. Cross regular sequences are one of the possible extensions of the well known regular sequences introduced by Sacks and Ylvisaker The introduced classes of random fields have the same smoothness and local behavior for each coordinate of components generated by a decomposition vector l. Therefore in the following, we use only approximation designs with the same within- and interdimensional knot distributions within the components. Formally, for the partition generated by a vector l = l 1,...,l k, we consider cross regular designs T N, defined by the functions h := h 1,...,h k and πn := n 1 N,...,n k N, as follows: h i h, n i = n, i = L 1 +1,...,L, = 1,...,k. We call the functions h 1,...,h k and πn withincomponent densities and intercomponent knot distribution, respectively. The corresponding property of a design T N is denoted by: T N is crsh,π,l. For a given cross regular sampling design, the hypercube D is partitioned into M = d =1 n disoint hyperrectangles D i, i = i 1,...,i d, 0 i k n k 1, k = 1,...,d. Let 1 d = 1,...,1 denote a d-dimensional vector of ones. The hyperrectangle D i is determined by the vertex t i = t 1,i1,...,t,id and the main diagonal r i = t i+1d t i, i.e., D i := { t : t = t i +r i s,s [0,1] d}, where denotes the coordinatewise multiplication, i.e., for x = x 1,...,x d and y = y 1,...,y d, x y := x 1 y 1,...,x d y d. For a random field X CD, define a multivariate piecewise linear interpolator MPLI with knots T N X N t := X N X,T N t = E η Xt i +r i η, t D i, t = t i +r i s, where η = η 1,...,η d and η 1,...,η d are auxiliary independent Bernoulli random variables with means s 1,...,s d, respectively, i.e., η Bes, = 1,...,d. Such defined interpolator is continuous and piecewise linear along all coordinates. Example. Let d =, N = 4, D = [0,1]. Then t = s, r = 1,1, X N t = E η Xη = X0,01 t 1 1 t +X1,0t 1 1 t +X0,11 t 1 t +X1,1t 1 t, and X N is a conventional bilinear interpolator see, e.g., Lancaster and Šalkauskas, We introduce some additional notation used throughout the paper. For sequences of real numbers u n and v n, we write u n v n if lim n u n /v n 1 and u n v n if there exist positive constants c 1,c such that c 1 u n v n c u n for n large enough. Results Let B β,m t,t R m +, 0 < β <, m N, denote an m-dimensional fractional Brownian field with covariance function 3. For any u R m +, we denote b β,m u := B β,m u s E η B β,m u η ds, [0,1] m 4

5 where η = η 1,...,η m, and η 1,...,η m are independent Bernoulli random variables η Bes, = 1,...,m. Then b β,m u is the squared IMSE of approximation for B β,m u t,t [0,1] m, by the MPLI with m observations in the vertices of unit hypercube. In the following theorem, we provide an exact asymptotics for the IMSE of a local stationary field approximation by MPLI when a cross regular sequence of sampling designs is used. Theorem 1 Let X B α l D,c be a random field approximated by the MPLI X NX,T N, where T N is crsh, π, l. Then where X X N v = D =1 v n α and H t := 1/h t L 1 +1,...,1/h t L, = 1,...,k. > 0 as N, c tb α,l H t dt > 0, Remark 1 If for the -th component, the uniform withincomponent knot distribution is used, i.e., h s = 1, s [0,1], then the asymptotic constant is reduced to where b α,l := b α,l 1 l. v = b α,l D c tdt, In Theorem 1, the approximation accuracy is determined by the sampling grid sizes n. The next theorem provides the asymptotically optimal intercomponent knot distribution for a given total number of observation points N. Denote by ρ := l i α i=1 i 1 = d α i=1 i 1 1, κ := k =1 v l /α, where d ρ is the harmonic mean of the smoothness parameters α, = 1,...,d. Theorem Let X B α l D,c be a local stationary random field approximated by the MPLI X NX,T N, where T N is crsh,π,l. Then X X N κρ k as N. 5 Nρ Moreover, for the asymptotically optimal intercomponent knot allocation, n,opt Nρ/α v 1/α κ ρ/α as N, = 1,...,k, 6 the equality in 5 is attained asymptotically. The above result agrees with the intuition that more points should be distributed in directions with lower smoothness parameters. Note that the optimal intercomponent knot distribution leads to an increased approximation rate. 5

6 Remark Let X B α l D,c with k = d and α i α for some i, = 1,...,d, and α := min i=1,...,d α i, i.e., ρ > α. Consider the approximation with uniform intercomponent knot distribution, n 1 = = n d N 1/d. Then by Theorem 1, we have X X N On the other hand, the sampling distribution 6 gives 1 N α/d. X X N 1 N ρ/ < 1 N α/d. Example 3. Let d = k =, α 1 = /3, α = 5/3. Then for n 1 = n, the approximation rate is N α/d = N 1/6 while using the asymptotically optimal intercomponent distribution we obtain the rate N ρ/ = N 1/4. < N 1/6. In general setting, numerical procedures can be used for finding optimal densities. However, in practice such methods are very computationally demanding. We present a simplification of the asymptotic constant expression for one-dimensional components. Further, in this case, we provide the exact formula for the density minimizing the asymptotic constant. For a random field X Bl α D,c, define the integrated local stationarity functions C t L := c tdt 1...dt L 1dt L +1...dt d, t L [0,1], = 1,...,k. [0,1] d 1 Moreover, for 0 < β <, let a β := β +1β Proposition 1 Let X Bl αd,c be a random field approximated by the MPLI X NX,T N, where T N is crsh,π,l. If for some, 1 k, l = 1, then for any regular density h, we have The density minimizing v is given by 1 v = a α C t L h t L α dt L. 0 h,opt t L = C t L γ 1 0 C, t τ L γ L [0,1], dτ L where γ := 1/1+α. Furthermore, for such chosen density, we get v,opt = a α C γ. In the subsequent proposition, we give an upper bound for the approximation error together with expressions for generating densities minimizing this upper bound, called suboptimal densities. Proposition Let X B α l D,c be a random field approximated by the MPLI X NX,T N, where T N is crsh, π, l. Then X X N =1 w n α as N, 6

7 where w = l 1+α / a α C t L h t L α dt, = 1,...,k. 6 0 The density minimizing w is given by C t L γ h,subopt t L = 1 0 C, t τ L γ L [0,1], dτ L where γ := 1/1+α, = 1,...,k. Furthermore, for such chosen densities, we get w,subopt = l 1+α / a α + 1 C 6 γ, = 1,...,k. Now we focus on random fields satisfying the introduced Hölder type conditions. In this case, we provide results for the uniform mean square norm of approximation error X X N. The following proposition provides an upper bound for the accuracy of MPLI for Hölder classes of continuous and continuously differentiable fields. Proposition 3 Let X CD be a random field approximated by the MPLI X N X,T N, where T N is crsh, π, l. i If X C α l D,C, then for positive constants c 1,...,c k. ii If X C 1,α l D, C, then for positive constants d 1,...,d k. X X N C X X N C =1 =1 c n α / d n 1+α / 7 8 Remark 3 It follows from the proof of Proposition 3 that 7 holds if c = α l 1+α / D α, = 1,...,k, where D := 1/min s [0,1] h s, = 1,...,k. Therefore the constants depend only on the parameters of the Hölder class and the corresponding sampling design. Similar formulas can be obtained for d 1,...,d k in 8. In addition, we provide the intercomponent knot distribution leading to an increased rate of the upper bounds obtained in Proposition 3. Remark 4 Let X CD be a random field approximated by the MPLI X N X,T N, where T N is crsh,π,l. i If X C α l D,C and n N ρ 0/α, = 1,...,k, where ρ 0 = k i=1 l i/α i 1, then X X N = ON ρ 0/ as N. ii If X C 1,α l D,C and n N ρ 1/+α, = 1,...,k, where ρ 1 = k i=1 l i/+α i 1, then X X N = ON ρ 1/ as N. 7

8 The approximation rates obtained in the above remark are optimal in a certain sense, i.e., the rate of convergence can not be improved in general for random fields satisfying Hölder type condition see, e.g., Ritter, 000. Moreover, these rates correspond to the optimal approximation rates for anisotropic Nikolskii-Hölder classes see, e.g., Yanie and Yongping, 000, which are deterministic analogues of the introduced Hölder classes. 3 Numerical Experiments In this section, we present some examples illustrating the obtained results. For given knot densities and covariance functions, first the pointwise approximation errors are found analytically. Then numerical integration is used to evaluate the approximation errors on the entire unit hypercube. Let δ N h,πt = δ N X,X N,T N h,π,lt := Xt X N X,T N h,π,lt, t [0,1] d, be the deviation field for the approximation of X by the MPLI with N knots, where T N is crsh,π,l, and write e N h,π := δ N h,π for the corresponding IMSE. We write h uni, to denote the vector of withincomponent uniform densities. Analogously, by π uni we denote the uniform interdimensional knot distribution, i.e., n 1 =... = n k. Example 4. Let D = [0,1] 3 and Xt = B α,l t, where α = 1/,3/ and l = 1,. Then X Bl α[0,1]3,c, with ct = 1,1, t [0,1] 3, k =, α = 1/,3/,3/. We compare behavior of e N h uni,π uni and e N h uni,π opt, where π opt given by Theorem. Observe that by using the asymptotically optimal intercomponent distribution, we obtain a gain in the rate of approximation. Figure 1 shows the fitted values of the squared IMSEs e N h uni,π uni and e N h uni,π opt in a log-log scale. In such scale, the slopes of fitted lines correspond to the rates of approximation. These 1.5 π uni π opt logn Figure 1: The fitted plots of e N h uni,π uni solid line, e N h uni,π opt dash line versus N in a log-log scale. plots represent the following asymptotic behavior: e Nh uni,π uni N 1/ N 1/ N 1/6, e Nh uni,π opt 0.445N 3/10 as N. 8

9 Example 5. Let D = [0,1] and define Xt = Xt 1,t to be a zero mean Gaussian field with covariance function 1 1 CovXt,Xs = t +0.1 s +0.1 exp t s. Then X Bl α[0,1],c with ct = c 1 t = / t + 0.1, t [0,1], α = 1, α = 1,1, l =, and k = 1. The field has one component, hence the uniform interdimensional knot distribution is used. Theorem provides the formula for the suboptimal withincomponent density. Figure a shows the fitted values of the squared IMSEs e N h uni,π uni and e N h subopt,π uni. Figure b demonstrates the convergence of the scaled squared approximation error N 0.5 e N h subopt,π uni to the asymptotic constant obtained in Theorem 1. Note that utilizing the suboptimal withincomponent density leads to a significant reduction of the asymptotic h subopt h uni logn a N b Figure : a The fitted plots of e N h uni,π uni dashed line and e N h subopt,π uni solid line versus N in a log-log scale. b The convergence of N 0.5 e N h subopt,π uni solid line to the asymptotic constant dashed line. constant, as compared to the uniform withincomponent knot distribution. 4 Proofs Proof of Theorem 1. First we investigate the asymptotic behavior of the approximation error e N t := Xt X N t for any t D i, i I, where I := {i = i 1,...,i d, 0 i k n k 1, k = 1,...,d}, when the number of knots N tends to infinity. Further, we find the asymptotic form of the IMSE 1/ e N := e N t dt for any positive continuous densities h 1,...,h k. We start by observing that e N t = EXt X N t = EE η Xt i +r i η Xt = E η,ξ EXt i +r i η XtXt i +r i ξ Xt D = 1 E η,ξe Xt i +r i η Xt +Xt i +r i ξ Xt Xt i +r i η Xt i +r i ξ, 9 where ξ is an independent copy of η. Further, the property together with the uniform continuity and positiveness of local stationarity functions c 1,...,c k imply that e N t = 1 r c t i E η,ξ i η s α r + i ξ s α r i η ξ α 1+q N,i t, 10 =1 9

10 where ε N := max{ q N,i t,t D i,i I} = o1 as N cf. Seleznev, 000. It follows from the definition and the mean integral value theorem that r i = h 1 w 1,i 1 n, 1 h w,i n,..., h d w d,i d n, w,i [t,i,t,i +1], = 1,...,d. d Denote by w i := w 1,i1,...,w d,id. Now the definition of crsh,π,l implies r i = 1 1 n h w L 1 +1,i,..., = 1 H w L 1 n +1 h w L,i L n i, = 1,...,k, where H t := 1/h t L 1 +1,...,1/h t L, = 1,...,k. Consequently, e N t = 1 =1 n α c t i E η,ξ H w i η s α + H w i ξ s α H w i η ξ α 1+o1 as N. Applying the uniform continuity of h yields e N t = 1 n α c t i E η,ξ H t i η s α + H t i ξ s α = =1 =1 H t i η ξ α 1+o1 n α c t i C α,l s ;H t i 1+o1 as N, where C α,l s ;H t i :=1 E η,ξ H t i η s α + H t i ξ s α H t i η ξ α = B α,l H t i s E η B α,l H t i η. Let D i = Di 1 Dk i and denote by D i the volume of hyperrectangle D i. Then e N = e N t dt = n α c t i C α,l s ;H t i dt 1+o1 i I D i i I D i =1 = n α c t i C α,l s ;H t i ds D i 1+o1 i I =1 D = c t i b α,l H t i D i 1+o1 as N. =1 n α i I Now the Riemann integrability of the functions c tb α,l H t, = 1,...,k, gives e N = =1 n α D c tb α,l H t dt 1+o1 = 10 =1 v n α 1+o1 as N.

11 Note that for any u R m +, b β,m u > 0, otherwise the fractional Brownian field is degenerated cf. Seleznev, 000. Consequently, v > 0, = 1,...,k. This completes the proof. Proof of Theorem. Note that by the inequality for the arithmetic and geometric means, with equality if only if 1 k =1 v n α ν 1 = v n α k =1 v n α 1/k, = 1,...,k. Hence, the equality is attained for ñ = νv 1/α, = 1,...,k. Let The total number of observations satisfies n = ñ νv 1/α as N. 11 d N = n 1 +1 n d +1 n i = i=1 k =1 n l = M as N. This implies that for the asymptotically optimal intercomponent knot distribution N M ν 1/ρ k =1 v l /α, and therefore, ν N ρ κ ρ as N. By equation 11, the asymptotically optimal intercomponent knot distribution is n Nρ/α v 1/α κ ρ/α as N, = 1,...,k. Moreover, with such chosen knot distribution, the equality in 5 is attained asymptotically. This completes the proof. Proof of Proposition 1. The proof is a straightforward implication of the assumptions and equation 10. The exact constant and the expression for the optimal density are due to Seleznev 000. Proof of Proposition. The first steps of the proof repeat those of Theorem 1. By 10, we have e N t = 1 r c t i E η,ξ i η s α r + i ξ s α r i η ξ α 1+o1 =1 1 r c t i E η,ξ i η s α r + i ξ s α 1+o1 =1 r = c t i E η i η s α 1+o1 as N. =1 11

12 For any nonnegative numbers a 1,...,a k and any α R +, the inequality holds, and consequently, e N t = =1 =1 c t i l α / c t i l α / α a i k α a α i 1 L i=1 m=l 1 +1 L m=l 1 +1 i=1 E η r i,m η m s m α 1+o1 r α i,m 1 sm α s m +1 s m s α m 1+o1. By the mean value theorem and the uniform continuity of withincomponent densities, we obtain e N t =1 c t i l α / L n α m=l 1 +1 Proceeding now to the calculation of the IMSE, we get e N = i I where D i e N t dt i I =1 h t i,m α 1 s m α s m +1 s m s α m 1+o1 as N. c t i l α / n α α +1α + = 1 0 L m=l 1 +1 h t i,m α α +1α + D i 1+o1, 1 s α s+1 ss α ds. Now the Riemann integrability of c th t m α, = 1,...,k, together with the definition of integrated local stationarity functions imply that e N =1 = = =1 =1 1 n α 1 n α 1 n α l α / a α + 1 L 6 m=l 1 +1 l α / a α l 1+α / a α L m=l 1 +1 c t i h t i,m α D i 1+o1 i I D c th t m α dt 1+o1 1 C t L h t L α dt L 1+o1 as N. The expression for the suboptimal density is due to Seleznev 000. This completes the proof. 0 Proof of Proposition 3. We start by proving i. Let X C α l [0,1]d,C and consider t D i, i I. Applying the Hölder condition 1 to equation 9 yields e N t = 1 E η,ξe Xt i +r i η Xt +Xt i +r i ξ Xt Xt i +r i η Xt i +r i ξ CE η r i η α = CE η ri η α C =1 L l α / =1 m=l 1 +1 E η r i,m η m s m α, 1

13 where the last inequality follows from 1. Furthermore, since max s [0,1] 1 s α s+1 ss α = α, we obtain e N t C L l α / =1 m=l 1 +1 r α i,m 1 sm α s m +1 s m s α m α l α / =1 L m=l 1 +1 By the regularity of the generating densities, we have that r i,m 1/n mmin s [0,1] h ms,i I,m = 1,...,d. Moreover, the definition of crsh, π, l implies the following uniform bound for the squared approximation accuracy X X N = max t D e Nt C =1 α l 1+α / D n α, with D = 1/min s [0,1] h s, = 1,...,k. Finally, we obtain the required assertion X X N C c =1 n α / where c := α l 1+α / D α > 0, = 1,...,k. For the smooth case, we use the multivariate Taylor formula to obtain the following representation of the deviation field δ n t := Xt X N t = E η 1 0, r α i,m. d X t i +ur i η sr i, η s du, t D i, t = t i +s r i, =1 where η = η 1,...,η d and η 1,...,η d are independent Bernoulli random variables, η Bes, = 1,...,d. Introducing an auxiliary uniform random variable U U 0,1 we get δ n t = = d E η,u X t i +Uη s r i r i, η s =1 d E η,u X t i,1 +Uη 1 s 1 r i,1,...,t i, +Uη s r i,,...,t i,d +Uη d s d r i,d =1 since for any = 1,...,d, X t i,1 +Uη 1 s 1 r i,1,...,t i,,...,t i,d +Uη d s d r i,d η s, E η X t i,1 +Uη 1 s 1 r i,1,...,t i,,...,t i,d +Uη d s d r i,d η s = E η1,...,η 1,η +1,...,η d X t i,1 +Uη 1 s 1 r i,1,...,t i,,...,t i,d +Uη d s d r i,d E η η s = 0. The triangle inequality and the condition 4 imply that e N t d =1 CV r 1+α / i,, for some positive constants V, = 1,...,d. Analogously to i, the required assertion follows from the regularity of the generating densities and the definition of crsh, π, l. This completes the proof. Acknowledgments The second author is partly supported by the Swedish Research Council grant and the proect Digital Zoo funded by the European Regional Development Fund. 13

14 References Abramowicz, K., Seleznev, O., 011. Spline approximation of a random process with singularity. J. Statist. Plann. Inference 141, Adler, R., Taylor, J., 007. Random fields and geometry. Springer, New York. Benhenni, K., 001. Reconstruction of a stationary spatial process from a systematic sampling. Lecture Notes-Monograph Series 37, Berman, S.M., Soourns and extremes of Gaussian process. Ann. Probab., de Boor, C., Gout, C., Kunoth, A., Rabut, C., 008. Multivariate approximation: theory and applications. An overview. Numer. Algorithms 48, 1 9. Brouste, A., Istas, J., Lambert-Lacroix, S., 007. On fractional Gaussian random fields simulations. J. Stat. Soft. 3, 1 3. Christakos, G., 199. Random field models in earth sciences. Academic Press, London. Hüsler, J., Piterbarg, V., Seleznev, O., 003. On convergence of the uniform norms for Gaussian processes and linear approximation problems. Ann. Appl. Probab. 13, Kuo, F., Wasilkowski, G., Woźniakowski, H., 009. On the power of standard information for multivariate approximation in the worst case setting. J. Approx. Theory 158, Lancaster, P., Šalkauskas, K., Curve and surface fitting: an introduction. Academic Press, London. Lifshits, M., Zani, M., 008. Approximation complexity of additive random fields. J. Complexity 4, Müller-Gronbach, T., Hyperbolic cross designs for approximation of random fields. J. Statist. Plann. Inference 66, Müller-Gronbach, T., Schwabe, R., On optimal allocations for estimating the surface of a random field. Metrika 44, Nikolskii, S., Approximation of functions of several variables and imbedding theorems. Springer-Verlag, Berlin. Pratt, W., 007. Digital image processing: PIKS scientific inside. Wiley-Interscience, New York. Ritter, K., 000. Average-case analysis of numerical problems. Springer-Verlag, Berlin. Ritter, K., Wasilkowski, G.W., Woźniakowski, H., Multivariate integration and approximation for random fields satisfying Sacks-Ylvisaker conditions. Ann. Appl. Probab. 5, Sacks, J., Ylvisaker, D., Designs for regression problems with correlated errors. Ann. Math. Statist. 37, Seleznev, O., Large deviations in the piecewise linear approximation of Gaussian processes with stationary increments. Adv. Appl. Prob. 8,

15 Seleznev, O., 000. Spline approximation of stochastic processes and design problems. J. Statist. Plann. Inference 84, Stein, M., Interpolation of spatial data. Springer-Verlag, New York. Su, Y., Estimation of random fields by piecewise constant estimators. Stochastic Process. Appl. 71, Yanie, J., Yongping, L., 000. Average widths and optimal recovery of multivariate Besov classes in L p R d. J. Approx. Theory 10, Zhang, X., Wicker, S.B., 005. On the optimal distribution of sensors in a random field. ACM Trans. Sen. Netw. 1,