arxiv:math/ v1 [math.pr] 5 Aug 2006

Transcription

1 Preprint SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION OF SUMS OF LOCALLY DEPENDENT RANDOM VARIABLES arxiv:math/68138v1 [math.pr] 5 Aug 26 By Adrian Röllin Stein s method is used to approximate sums of discrete and locally dependent random variables by a centered and symmetric Binomial distribution. Under appropriate smoothness properties of the summands, the same order of accuracy as in the Berry-Essen Theorem is achieved. The approximation of the total number of points of a point processes is also considered. The results are applied to the exceedances of the r-scans process and to the Matérn hardcore point process type I. 1. Introduction. The approximation of sums of dependent random variables by the standard normal distribution has been investigated in a large variety of settings. The accuracy of approximation is most often measured by the Kolmogorov and Wasserstein metrics. The use of stronger metrics typically requires that some smoothness -condition must be satisfied. In this paper, under the assumption of a general local dependence structure, we study the approximation of sums of discrete random variables by a symmetric and centered Binomial distribution. This distribution serves as replacement for the normal distribution in a discrete setting. Under some general smoothness property of the summands, the same order of accuracy as in the Berry-Essen Theorem can be achieved, but now for the total variation metric. We also examine another metric, from which local limit approximations can be obtained. In the setting of independent summands, approximation by a centered Poisson distribution has been successfully adopted by Čekanavičius and Vaǐtkus [11] and Barbour and Čekanavičius [4]. However, for dependent summands, applications were limited to simple examples; first attempts were made by Barbour and Xia [5] and Čekanavičius and Vaǐtkus [11]. In contrast, the results in this paper are of general nature and allow a wide range of applications. The proofs are based on Stein s method for distributional approximation. A main idea, introduced in Röllin [2], is to use interpolation functions to represent the Stein operator of a discrete distribution as the Stein operator of a continuous distribution. In the case of the Binomial, this then allows the application of standard techniques in Stein s method for normal approximation. A careful analysis of the remainder terms then shows how a suitable smoothness condition can be exploited, to obtain total variation error bounds. The paper is organized as follows. In the next section, we introduce the main technique in the simple case of independent summands. In section 3 these results are extended to locally dependent summands and section 4 shows their application in some examples. Section 5 contains some technical lemmas. Partially supported by Schweizerischer Nationalfondsprojekt /1 AMS 2 subject classifications: Primary 6F5 Keywords and phrases: Stein s method; total variation metric; Binomial distribution; local dependence 1

2 2 A. RÖLLIN 1.1. Notation. Denote by Bin, p) the Binomial distribution with n trials of probability p each. Denote by Bin, p) the centered Binomial distribution, i.e. a Binomial distribution shifted by np. Note that this distribution does not necessarily lie on the integers, but on a lattice of Ê with span 1. Throughout the paper, we shall be concerned with two metrics for probability distributions, the total variation metric d TV and the local limit metric d loc, where, for two probability distributions P and Q, ) d TV P, Q := sup PA) QA), A Ê ) ) ) d loc P, Q := sup P [x, x 1) Q [x, x 1). x Ê For simplicity, we will often use the notation d l, where l = 1 will stand for d TV and l = 2 for d loc. We denote by the supremum norm if applied to functions, and the variation norm if applied to measures. Let δ x denote the unit mass at x Ê, and the convolution of measures. Define for any measure µ and any l Æ := 1, 2,... } Note that for measures µ and λ, D l µ) = µ δ 1 δ ) l. D 1 µ ) ) 1.1) = 2 d TV µ, µ δ1, 1.2) D 2 µ λ) D 1 µ)d 1 λ). Furthermore, define x := x x to be the fractional part of x Ê, and x) = x Basic setup. Consider a sum of the form W = i J ξ i, where W takes its values in a lattice of Ê with span 1. The expectation of W has no influence on the quality of the approximation, and we therefore assume without loss of generality that W = ; this can always be accomplished by subtracting the expectation from each individual summand. Each of the summands may now take its values on a different lattice; this, however, will result in no further complications. To approximate W by a centered binomial distribution, we have to choose n in such a way that the variance of Bin, 1/2) is as close to the variance of W as possible. As n has to be integer, this is only possible up to a rounding error. However, the symmetric and centered Binomial distribution thus chosen will in general take its values on a different lattice from W and the total variation distance will become 1. To circumvent this problem, we introduce an additional parameter t and approximate W by a centered Binomial distribution with success probability 1/2 t instead t being small), to be able to match not only the variance but also the lattice. Hence, to put the above in a rigorous form, we will make the following assumptions if not otherwise stated: Assumptions G: Let J be a finite set and let ξ i, i J} be a collection of random variables with ξ i = for all i J and assume that there are numbers a i Ê; i

3 SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION 3 J} such that almost surely ξ i a i. Let W = i J ξ i; then W = and almost surely W a for a := i J a i. Assume that σ 2 := VarW > 1. Define now δ := 4σ 2 and t := a 2σ 2 δ/2 /4σ 2 δ). Clearly, 4σ 2 δ = 4σ 2, and by definition the distribution Bi 4σ 2, 1/2 t ) has expectation ; it is also easy to check that it takes values in a. From the above definition, we see that t is only of order Oσ 2 ), which is rather small in the setting that we are concerned with; Corollary 2.3 shows how to obtain results without t, using Lemma Sum of Independent Random Variables. First, we examine the case of independent discrete summands. Previous work on total variation approximation has been concerned with the compound Poisson distribution see Le Cam [17] and Roos [21] and references therein), the signed compound Poisson distribution see Čekanavičius [9] and references therein), the Poisson distribution see Barbour et al. [7]), the centered Poisson distribution see Čekanavičius [1], Čekanavičius and Vaǐtkus [11], Barbour and Xia [5] and Barbour and Čekanavičius [4]) and some more general distributions see Brown and Xia [8]). We present the theorem below to demonstrate the main technique in a simple setting, noting that it also follows as a consequence of Theorem 3.1. Theorem 2.1. Let ξ i ; i J} be independent and satisfy Assumptions G. Then, if the ξ i have finite third moments, d l L W), Bi 4σ 2, 1/2 t )) σ 2 ) c l,i ρ i 1.75, l = 1, 2, where ρ i = σ 3 i 1 2 ξ i 3, σ 2 i = Varξ i and c l,i = D l L W ξ i ) ). It is clear that the above bound is useful only if the c l,i are small. In the case of n identically distributed random variables, we need c 1,i = o1) as n for asymptotic approximation in total variation, and in order to deduce a local limit theorem we must have c 2,i = on 1/2 ). This is however always the case if D 1 X 1 ) < 2 this corresponds to the usual condition in the LLT that X 1 must not be concentrated on a lattice with span greater than 1), as can be seen from 5.9) 5.1), and we then even have c l,i = On l/2 ) for l = 1, 2. Before proving the theorem, we start with a short summary of Stein s method for Binomial approximation; for details see also Stein [22] and Ehm [16]. Denote by FM) the set of all real valued measurable functions on some given measure space M. A Stein operator B : F ) F ) for the Binomial distribution Bin, p) is characterized by the fact that, for any integer valued random variable W, 2.1) Bg)W) = for all bounded g F ) W Bin, p), and a possible choice is 2.2) Bg)z) = qzgz 1) pn z)gz), for all z, i J

4 4 A. RÖLLIN where, as usual, we put q = 1 p. Let h F ) be a bounded function. Then, the solution g = g h to the Stein equation 2.3) Bg)z) = I[ z n] hz) hy ) }, for all z, where Y Bin, p), is also bounded. If the functions h are of the form hz) = h A z) = I[z A], A, we have the uniform bound 2.4) g A 1 pn1 q n1, n 1)pq where gz) := gz 1) gz), and the same bound holds for g b}, b ; see Ehm [16]. Now, for all z, we can write I[z A] È[Y A] = Bg A )z) I[z /...n}] I[z A] È[Y A] ), and thus, for any integer valued random variable V, ) d TV L V ), Bin, p) = sup È[V A] È[Y A] A 2.5) sup BgA )V ) [ È V n/2 > n/2 ]. A We now construct a Stein operator for the centered Binomial distribution Bin, p) on the lattice np. For any function g F ) define the function ĝ F np) by ĝw) := gw np) for w np. Then the Stein operator is defined as 2.6) ˆBĝ)w) := Bg)w np) = pw np)gw np) qw np)gw 1 np) npgw np) = w pĝw) qĝw 1) ) npq ĝw 1). for all w np. Thus, for W = V np, an inequality corresponding to 2.5) holds, namely ) d TV L W), Bin, p) 2.7) sup B np ˆBĝB )W) È [ W np 1/2) > n/2 ]. An equivalent inequality holds for the d loc metric, but the supremum is taken only over the sets b}, b np. Under the assumptions of the theorem, n = 4σ 2 = 4σ 2 δ and p = 1/2 t, and 2.6) becomes 2.8) ˆBĝ)w) = wθĝw 1) σ 2 ĝw 1) t 2 4σ 2 δ) wt δ/4 ) ĝw 1), where Θĝw) := 1 2ĝw 1) ĝw) ). Since σ 2 > 1, the bound 2.4) simplifies to 2.9) ĝ B 1 σ 2.

5 SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION 5 To see this, note that t < 1/ 4σ 2 = 1/n and n = 4σ 2 5. Then from 2.4) we have ĝ B 1 n 1)pq = 1 n 1)1/4 t 2 ) 4n 2 n 1)n 2 4) 4 n 1 σ 2. Lemma 2.2. Assume the conditions of Theorem 2.1. Define A : F a) F a) by Aĝ)w) := wθĝw 1) σ 2 ĝw 1), w a, ĝ F a). Then, 2.1) Aĝ)W) ĝ ) c 1,i ρ i ĝ c 2,i ρ i ). i J i J Proof. For every w a and x [, 1) define 2.11) fw x) := Θĝw 1) x ĝw 1) 1 2 x2 2 ĝw 1). One easily checks that f C 1 and fw) = Θĝw 1) and f w) = ĝw 1), hence 2.12) Aĝ)w) = wfw) σ 2 f w), for all w a. Furthermore, f is absolutely continuous, hence f exists almost everywhere. Choose f to be the function 2.13) f w x) = 2 ĝw 1) for all w a, x < 1. We can now apply the usual Taylor expansion cf. Reinert [19], Theorem 2.1), but with a refined estimate of the remainder terms. Write W i = W ξ i, i J; then ξ i fw) = ξ i fw i ) ξ 2 i f W i ) ξ 3 i σ 2 i f W) = σ 2 i f W i ) ξ i σ 2 i s)f W i sξ i )ds, f W i sξ i )ds, and hence, using the independence of ξ i and W i and that ξ i =, 2.14) ξi fw) σ 2 i f W) } ξ 3 i 1 ξ i σ 2 i 1 s) [ f W i sξ i ) ξi ] ds 1 [ f W i sξ i ) ξi ] ds. Note now that for any real valued random variable U taking values on a lattice with span 1, we obtain together with 2.13) 2.15) f U z) ) ĝ D 1 L U) )) ĝ D 2 L U) )),

6 6 A. RÖLLIN for all z Ê. Thus, from 2.14) and 2.15), ξi fw) σi 2 f W) } 2.16) ĝ D 1 L W i ) ) σ 3 i 1 2 ξ i 3)) Now, using 2.12) we have AĝW) } i J ĝ D 2 L W i ) ) σ 3 i 1 2 ξ i 3)). ξi fw) σ 2 i f W) } and with 2.16) the lemma is proved. Proof of Theorem 2.1. Recall that, by Assumptions G, the distributions L W) and Bi 4σ 2, 1/2 t) ) are concentrated on the same lattice. Thus, using 2.7) and the form 2.8) of the Stein operator, and applying the left side of the minimum in 2.1) to the first part of 2.8) with the bound 2.9) gives d TV L W), Bi4σ 2 δ, 1/2 t) ) 2.17) i J c 1,iρ i σ 2 t2 4σ 2 δ) σt δ/4 σ 2 [ È W 2σ 2 1 ]. To bound the middle part of 2.17) note that t < 4σ 2 δ) 1 and δ < 1. Thus, recalling that σ 2 > 1, we obtain the simple bounds t 2 4σ 2 δ) < 4σ 2 δ) 1 1/4, σt σ/4σ 2 δ) 1/4, δ/4 1/4. Applying Chebyshev s inequality on the last term of 2.17) we obtain È [ W 2σ 2 1 ] σ 2 2σ 2 1) 2 1 σ 2. The d loc case is analogous, using the right side of the minimum in 2.1) instead and the remark after 2.4). Note that in the next corollary we do not assume that the ξ i have expectation zero. Corollary 2.3. Let W be the sum of independent and integer valued random variables ξ i, i J} with σ 2 i = Varξ i and Then, if σ 2 > 1, v i = min 1/2, 1 d TV L ξi ), L ξ i 1) )}. d TV L W), Bi 4σ 2, 1/2) δ s ) 2 σ 3 i 1 2 ξ i 3) σ 2 V v ) 1/2 d loc L W), Bi 4σ 2, 1/2) δ s ) 8 σ 3 i 1 2 ξ i 3) σ 2 V 4v ) σ 1.25σ 2, σ σ 1 σ 2, where s := µ 4σ 2 /2, µ = W, V = i J v i and v = max i J v i.

7 SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION 7 Proof. Define W = W µ, and let t be defined with respect to W, taking a = µ. Then, as the metrics d l are shift invariant, d l L W), Bi 4σ 2, 1/2) δ s ) = dl L W ), Bi 4σ 2, 1/2) δ s µ ) d l L W ), Bi 4σ 2, 1/2 t )) d l Bi 4σ 2, 1/2 t), Bi 4σ 2, 1/2) ) =: R l 1 R l 2, since Bi 4σ 2, 1/2 t) δ s δ µ = Bi 4σ 2, 1/2 t ). Applying Lemma 5.2 to R l 2 with the fact that t 4σ 2 δ) 1 gives 2.18) R 1 2 σ 1 1 2σ) 1 4σ 2 ) 1), R 2 2 σ σ) 1). Define now c l = max i J D l L W ξ i ) )}. Application of 5.9)-5.1) yields 2.19) c 1 2 V v ) 1/2, c 2 8 V 4v ). Thus, application of Theorem 2.1 to R1 l proves the corollary. 3. Locally dependent random variables. In this section we present the main results of the paper. We exploit a finite local dependence structure as presented in Chen and Shao [12]. In the context of Stein s method for normal approximation, it has been successfully applied to a variety of problems; see for example Barbour et al. [6], Dembo and Rinott [15] and Barbour and Xia [3]. Note that Barbour et al. [6] use a slightly more general dependence structure, often yielding crucial improvements when approximating sums of dissociated random variables by the normal distribution. The generalization of Theorem 3.1 is straightforward, yet somewhat tedious, and we therefore use the simpler dependence structure of Chen and Shao [12]; see the Appendix for the more general version, but without proof. Let ξ i ; i J} be a collection of random variables satisfying Assumptions G. For convenience, let ξ A denote ξ i ; i A} for every subset A J. Assume further the following dependence structure: For every i J there are subsets A i B i J such that ξ i is independent of ξ A c i, and ξ Ai is independent of ξ B c i. Define η i = j A i ξ j and τ i = j B i ξ j. Theorem 3.1. With W as above, 3.1) d l L W), Bi 4σ 2, 1/2 t )) ) σ 2 ϑ l,i 1.75, l = 1, 2, i J where 3.2) ϑ l,i = 1 2 ξ i η 2 i Dl L W ξ i, η i ) )} ξ i η i τ i η i ) D l L W ξ i, η i, τ i ) )} ξ i η i τ i D l L W τ i ) )} If further there are constants c l,i such that almost surely 3.3) D l L W ξ Bi ) ) c l,i,

8 8 A. RÖLLIN then 3.4) ϑ l,i c l,i 1 2 ξ iη 2 i ξ i η i τ i η i ) ξ i η i τ i ). Proof. Estimate 3.4) is immediate. Following the proof of Theorem 2.1 and using Lemma 3.2 below, 3.1) is proved. Note that Theorem 2.1 follows from Theorem 3.1 with the choices A i = B i = i}. Lemma 3.2. Assume the conditions of Theorem 3.1. Define A : F a) F a) as in Lemma 2.2. Then, 3.5) Aĝ)W) ĝ ) ϑ 1,i ĝ ϑ 2,i ). i J i J Proof. We follow the proof of Lemma 2.2 right up to the end of the paragraph of 2.13). Note now that 3.6) σ 2 = i J ξ i η i } and that, by Taylor expansion, almost surely 3.7) ξ i fw) = ξ i fw η i ) ξ i η i f W η i ) ξ i η 2 i ξ i η i f W η i ) = ξ i η i f W τ i ) ξ i η i τ i η i ) 1 1 f W η i sη i )ds, f W η i sτ i η i ))ds, 1 ξ i η i }f W) = ξ i η i }f W τ i ) ξ i η i }τ i f W sτ i )ds. Now, using the facts that ξ i =, that ξ i is independent of W η i and that η i is independent of W τ i, we obtain from 3.6) and 3.7) that WfW) σ 2 f W) } = ξ i fw) ξ i η i }f W) } i J = 1 ξ i ηi 2 1 s) f W η i sη i ) } ξi, η i ds i J ξ i η i τ i η i ) 1 f W τ i sτ i η i )) ξi, η i, τ i } ds } 1 ξ i η i τi f W τ i sτ i )) } τi ds }. With 2.12) and 2.15) the lemma follows.

9 SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION 9 We now give a point process version of Theorem 3.1, exploiting mainly the same dependency structure as before. Theorem 3.3. Let Φ be a simple point process on a Polish space J with mean measure µ. For all points α J, assume that there are measurable subsets A α B α J, such that for every α J 3.8) 3.9) 3.1) L Φ α A c α )) = L ΦA c α )), Φ α A α ) and Φ α B c α) are independent, ΦA α ) and ΦBα c ) are independent, where Φ α denotes the Palm process at point α. Then, for W = ΦJ) µj) and if σ 2 > 1, 3.11) d l L W), Bi 4σ 2, 1/2 t )) σ 2 ϑ l α)µdα) 1.75σ 2, l = 1, 2. α J where, with Φ := Φ µ and Φ α := Φ α µ, 3.12) ϑ l α) = Φ α A α ) Φ B α ) D l ΦB c α ) ΦBα ) )} 1 2 Φ αa α ) 2 D l Φ α A c α) Φα A α ) )} 1 2 Φ A α ) 2 D l ΦA c α) ΦAα ) )} Φ α A α)φ α B α \ A α ) D l Φ α B c α ) Φ α A α ), Φ α B α ) )} Φ A α )Φ B α \ A α ) D l ΦB c α ) ΦAα ), ΦB α ) )}. Proof. Following the proof of Theorem 2.1 and Lemma 2.2, it is clear that we only have to bound WfW) σ 2 f W) } for f defined as in 2.11). In what follows, all integrals are taken over α J} if not otherwise stated. Note first that, because of 3.8), σ 2 = ΦJ)Φ J)} = µdα) Φ α A α) Φ α Ac α )} = µdα) Φ α A α) and hence with Taylor expansion σ 2 f W) = µdα) Φ α A α) f Φ Bα c )) µdα) Φ αa α ) Φ B α ) =: R 1 R 2. 1 } f Φ Bα) c tφ B α ))dt

10 1 A. RÖLLIN Now, again by Taylor, WfW)} = µdα) [ fφ α J)) fφ J)) ] = µdα) [ fφ αa c α)) fφ A c α)) ] µdα) [ Φ αa α )f Φ αa c α)) } Φ A α )f Φ A c α)) }] [ 1 µdα) Φ α A α) 2 =: R 3 R 4 R 5 1 t)f Φ α Ac α ) tφ α A α) ) } dt 1 Φ A α ) 2 1 t)f Φ A c α ) tφ A α ) ) } ] dt and R 4 = µdα) [ Φ α A α)f Φ α Bc α ))} Φ A α )f Φ B c α ))}] [ µdα) Φ αa α )Φ αb α \ A α ) =: R 6 R 7. 1 Φ A α )Φ B α \ A α ) 1 f Φ αbα) c tφ αb α \ A α ) ) } dt f Φ Bα) c tφ B α \ A α ) ) } ] dt Using 3.8) 3.1), we see that R 3 = and R 1 = R 6, hence WfW) σ 2 f W) } R2 R 5 R 7. With 2.15) we finally obtain R 2 ĝ µdα) Φ α A α ) Φ B α ) D 1[ L ΦBα) c ΦBα ) )]}, R ĝ µdα) [ Φ α A α) 2 D 1[ L Φ α A c α ) Φ α A α ) )]} R 7 ĝ Φ A α ) 2 D 1[ L ΦA c α) ΦAα ) )]}], µdα) [ Φ αa α )Φ αb α \ A α ) D 1[ L Φ α Bα) c Φα A α ), Φ α B α ) )]} Φ A α )Φ B α \ A α ) D 1[ L ΦB c α) ΦAα ), ΦB α ) )]}]. To obtain ϑ 2, just replace g by g and D 1 by D 2 in the above bounds. Corollary 3.4. Let Φ be a simple point process satisfying 3.8) 3.1). If there is further a function c l α), such that for µ-almost every α J almost surely 3.13) D l[ L ΦJ) Φ Bα )], D l [ L Φ α J) Φ α Bα )] cl α), l = 1, 2,

11 then 3.12) satisfies SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION ) Φ ϑ l α) c l α)[ α A α ) Φ B α ) 1 2 Φ αa α ) Φ A α ) 2 ] Φ α A α)φ α B α \ A α ) Φ A α )Φ B α \ A α ) c l α) [1.5 Φ α A α )Φ α B α ) } 1.5 ΦA α )ΦB α ) } ] 6µA α )µb α ) 4µB α ) Φ α B α ). 4. Applications. In what follows, we calculate only rough bounds, leaving much scope for improvement. In particular, we replace the moments in the estimates by almost sure bounds Exceedances of the r-scans process. We follow the notation of Dembo and Karlin [14]. Let X 1, X 2,..., X nr 1 be independent and identically distributed random variables with distribution function F. Define the r-scan process R i = r 1 k= X ik, i = 1, 2,...,n and further W i = I[R i a] for a Ê. We are interested in the number N = n i=1 W i, that is the number of R i not exceeding a. With p = W i = È[R 1 a], we have N = np and r 1 ) 4.1) σ 2 = Var W = np 1 p 2 1 d/n)ψd), where ψd) = È[R d1 a R 1 < a] p. Poisson approximations for the r-scan process have been extensively studied by Dembo and Karlin [14]. Normal approximation has been considered by Dembo and Rinott [15]; in particular they show, that, for fixed r and a, N converges in the Kolmogorov metric to the normal distribution with rate On 1/2 ). In the next theorem we achieve the same rate in total variation, and also a rate for the corresponding local limit approximation. Theorem 4.1. Assume that F is continuous, F) =, and Fx) < Fy) for all x < y, and let a > be fixed. Then, for all n such that σ 2 > 1, d=1 d l L N np), Bi 4σ 2, 1/2 t )) C l n l/2, l = 1, 2, where the constants C 1 and C 2 are independent of n and can be extracted from the proof. Proof. We apply Theorem 3.1 for W = n i=1 ξ i = n i=1 W i p). We can set A i = i r 1,...,i r 1} 1,...,n}, B i = i 2r 2,...,i 2r 2} 1,...,n}.

12 12 A. RÖLLIN Then, as A i 2r 1, B i 4r 3 and B i \ A i 2r 2, the following rough bounds are obvious: thus ξ i η 2 i 2r 1) 2, ξ i η i τ i η i ) 2r 1)2r 2), ξ i η i τ i 2r 1)4r 3), 4.2) ϑ l,i c l,i 16r 2 2r 6 ) Consider now the block B 1 = 3r 2 i=1 W i, and assume that the values B 1 = X 1,..., X r 1 ) and B 2 = X 3r 1,..., X 4r 2 ) are given. Define the events A := a/r < X r,...,x 2r 2, X 2r1,..., X 3r 2 ar 1)/r 2, < X 2r a/2r 2 ) } A := a/r < X 2r 1 ar 1)/r 2}, A 1 := < X 2r 1 a/2r 2 ) }. Due to the conditions on F and independence it is clear that p j := È[A A j ] > for j =, 1. Note now that R r = 2r 1 i=r X i > a on A A, R r < a on A A 1. Note further that R s < a for all s = r 1,...,2r 1 on A A A 1 ). Hence 2r 1 i=r W i = r 1 on A A, 2r 1 i=r W i = r on A A 1. It easy to see now by a coupling argument that 1 2 D1 L B 1 ) ) 1 p p 1 ) < 1. Noting that by sequentially stringing together blocks like B 1, we can have m := n/3r 2) such blocks, which are independent given all the borders B i. Furthermore, for every i, the R j in B i depend on the X k of at most two such blocks. Therefore, defining Z = B 1,..., B m ) and using 5.11) and 5.12), D 1 L W B i, i = 1,...,m) ) D 2 L W B i, i = 1,...,m) ) 2 min1/2, p, p 1 }m 2) ) 1/2 =: c 1,i, 8 min1/2, p, p 1 }m 4) =: c 2,i. Clearly, c l,i = On l/2 ). Hence, putting this, 4.1) and 4.2) into 3.1), the theorem follows.

13 SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION Matérn hard-core process type I. We approximate the total number of points of the Matérn hard-core process type I introduced by Matérn [18]. We use rectangular instead of the usual circular neighborhoods. Let Φ be the process on the d-dimensional cube J = [, 1) d Ê d defined as ΦB) = τ I[X i B]I[X j / K r X i ) for all j = 1,...,τ, j i], i=1 where τ Poλ) and X i ; i Æ} is a sequence of independent and uniformly distributed random variables on J and where, for x = x 1,..., x d ) J and r >, K r x) denotes the d-dimensional closed cube with center x and side length r. To avoid edge effects, we treat J as a d-dimensional torus, thus identifying any point outside J by the point in J which results in coordinate-wise shifting by 1. The process Φ is thus a thinned Poisson point process with rate λ having all points deleted which contain another point in their K r neighborhood. For the mean measure µ of Φ we obtain 4.3) dµx) dx = λe c. We are now interested in the distribution of ΦB) when r is small an λ large. Theorem 4.2. Put W := ΦJ) µj) and let a > be a fixed real number. Then, for every λ and r such that λr d = a and σ 2 := VarW > 1, d l L ΦJ) µj)), Bi 4σ 2, 1/2 t )) C l λ l/2, l = 1, 2, for constants C 1 and C 2 which are independent of λ and can be extracted from the proof. Proof. We apply Corollary 3.4. We can take A x = K 2r x) and B x = K 4r x) and check that the conditions 3.8) 3.1) are fulfilled. Some calculations show that the reduced second factorial moment measure M satisfies dmx) dx if x K r ), = λ 2 e λ Kr) Krx) if x K 2r ) \ K r ), λ 2 e 2a if x / K 2r ), compare with Daley and Vere-Jones [13, pp. 367, 373]. Thus, MJ) λ 2 e 2a 1 r d ) and 4.4) σ 2 = λe a MJ) µj) 2 λe a 1 ae a ). Since we can have at most 7 d points of Φ in B x, we obtain from 3.14) the rough estimate 4.5) ϑ l x) 26 7 d c l x),

14 14 A. RÖLLIN 5r }} K 3r x) x K 6r x) } } 6r Fig 1. Matérn hard-core process type I: Given that the process Φ is known on the borders l M K 6r x l ) \ K 5r x l ) grey area), the boxes Φ K6r x l ), l M, are independent. where c l ) is as in 3.13). To estimate c l x) write K r = K r ). We have È [ ΦK r ) = ΦK7r\K 5r ] Po λ Kr ) } = e a =: p, È [ ΦK r ) = 1 Φ K7r\K 5r ] Po λ K3r \ K r ) } Po λ K r ) 1} = ae 3da =: p 1. Hence, by a coupling argument, 4.6) 1 2 D1 L ΦK 6r ) Φ K7r\K 5r ) ) 1 p p 1 ) < 1. Let now x be arbitrary. Divide the space J into boxes of side length 6r, centered around x see Figure 1). With m := 1/6r), we can have m d such boxes plus a remainder. Denote this remainder by J R and denote by x l, l 1,..., m} d =: M the centers of the boxes where x 1,...,1 = x. Note now that, given Φ on all the borders K 6r x l ) \K 5r x l ), l M grey area in Figure 1), the random variables ΦK 6r x l )), l M, are independent and satisfy inequality 4.6). Furthermore, Φ J\K6rx) is independent of Φ Bx, and therefore, defining Z = ) Φ K6rx l )\K 5rx l )) l M, Φ J R and using 5.11) and 5.12), we obtain 4.7) 4.8) D 1 L ΦJ) Φ Bx )) 2 min1/2, p, p 1 }m d 2) ) 1/2 =: c 1 x), D 2 L ΦJ) Φ Bx )) 8 min1/2, p, p 1 }m d 3) =: c 2 x). Noting that almost surely L [ Φ x J \K6r x) ) Φ x Bx ] = L [ Φ J \K6r x) )], we see that 4.7) and 4.8) hold also for Φ x, thus c l x) satisfies 3.13). Now, recalling that a = λr d is constant, we have c l x) = Oλ l/2 ). Hence, putting this and 4.3) 4.5) into 3.11), the theorem follows.

15 5. Appendix. SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION Properties of the solution to the Stein equation. Lemma 5.1. For any indicator function hz) = I[z A], z, A, the solution g = g h to the Stein equation 2.3) satisfies 5.1) g 1 npq) 1/2. Proof. We apply the generator method introduced by Barbour [1]. For any function f :,...,n } Ê, define 5.2) Af)z) = B f) ) z) = qzfz 1) qzpn z) ) fz)pn z)fz1), which is the infinitesimal generator of a pure jump Markov process. A solution g to 2.3) is now given through ψz) = h Y z t) ) hy ) } dt, for z,...,n}, and gz) = ψz) for z,...,n 1} and gz) = else, where Y z is a Markov process with generator A starting at point z, and Y is a random variable having the stationary distribution Bin, p). Now, we have for z,...,n 1}, 5.3) ψz) = h Y z t) ) h Y z1 t) )} dt. We now fix z and construct a coupling of Y z and Y z1 to bound 5.3). Let thereto X i) k t), k 1,..., n}, i, 1}, be independent Markov processes with state space, 1}, starting in point i and having jump rate p if the process is in and q otherwise. It is easy to see by the Kolmogorov differential equations that 5.4) X 1) k t) Bep qe t ), X ) k t) Bep pe t ) where Bep) denotes the Bernoulli distribution with success probability p. Let τ be the minimum of the first jump times of the two processes X ) z1 and X1) z1, and define a new process X 1) z1 if τ > t, Xt) = X ) z1 if τ t, describing the well-known Doeblin coupling. Then, let z n 5.5) Y z = X 1) k X ) k, Y z1 = Y z X ) z1 Xt), k=1 k=z1 and one proves that Y z and Y z1 are Markov processes with generator 5.2). Hence, we can write 5.3) as 5.6) ψz) = e t hy z ) } dt,

16 16 A. RÖLLIN since τ is exponentially distributed with rate 1. The bound g 1 is now immediate from 5.6), thus we may assume that npq > 1. Note that, from 5.4) and 5.5), L Y z ) = Biz, p qe t ) Bin z, p pe t ), and hence, from Barbour and Jensen [2, Lemma 1], 5.7) D 1 L Y z ) ) VarY k ) 1/2 zp qe t )q qe t ) n z)p pe t )q pe t ) ) 1/2 npq1 e t ) ) 1/2. Note also that for h := h 1/2 5.8) hyz ) } = hyz ) } D 1 L Y z ) ) /2. Thus, applying 5.8) on 5.6) and using 5.7), s ψ e t dt 1 2 npq s e t 1 e t dt. Choosing s = ln 1 npq) 1) and computing the integrals proves the lemma Change of the success probabilities. Lemma 5.2. For every n Æ, < p < 1 and 1 p) < t < p ) n d TV Bin, p t), Bin, p) t p t ) p t)q t) pq pq pq n ) ) 1 p t p t)q t) d loc Bin, p t), Bin, p) t pq pq n Proof. We use Stein s method. If W Bin, p t), we obtain from 2.1) and 2.2) 1 p)wgw 1) pn W)gW) } = tw gw 1) tngw) } for every bounded function g F ). The left side is just the Stein operator for Bin, p) hence, taking g = g A obtained by solving 2.3) for Bin, p), with the bounds 2.4) and 5.1) the d TV -bound follows, noting also that W W VarW. With the remark after 2.4), the d loc -bound is proved Smoothing properties of independent random variables. In several parts of this paper, we have the situation that we need to estimate D m U), m = 1, 2, for some integer valued random variable U, being a sum of some other random variables. If the U is a sum of independent random variables, we can proceed as follows. Assume that U = n i=1 X i, where the X i are independent. Defining v i = min 1 2, 1

17 SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION D1 X i )} and V = i v i we obtain from Barbour and Xia [5, Proposition 4.6] the bound 5.9) D 1 U) 2 V 1/2. Define further v = max i v i. Now it is always possible to write U = U 1) U 2) in such a way that the analogously defined numbers V 1) and V 2) satisfy V k) V/2 v, k = 1, 2. Using 1.2) and 5.9), we obtain 5.1) D 2 U) D 1 U 1)) D 1 U 2)) 4 V 1) V 2)) 1/2 8 V 2v ) Smoothing properties of conditional independent random variables. In most applications, U is a sum of dependent summands and we can not apply 5.9) and 5.1) directly. However, assuming that there is a random variable Z on the same probability space as U such that L U Z = z) can be represented as a sum of independend summands, say X z) i, i = 1,...,n z, for each z that Z can attain, we can still apply 5.9) and 5.1), and we obtain 5.11) 5.12) D 1 U) [D 1 U) Z] } 2 V 1/2 Z }, D 2 U) [D 1 U) Z] } 8 V Z 2vZ ) where, for each z, V z and v z are the corresponding values as defined in subsection 5.3 with respect to the X z) i. 6. Appendix. We now give a generalization of Theorem 3.1. The proof is omitted, because it runs analogously to the proof of Theorem 3.1; see also Barbour et al. [6]. Suppose that a random variable W satisfies Assumptions G and assume that there are sets K i J, i I, and square integrable random variables Z i, Z ik and V ik, k K i and i I, as follows: }, 6.1) 6.2) 6.3) W = W i Z i, i I, where W i is independent of ξ i, Z i = Z ik, k K i W i = W ik V ik, i I, k K i, where W ik is independent of the pair X i, Z ik ). Theorem 6.1. With W as above, 6.4) d l L W), Bi 4σ 2, 1/2 t )) σ 2 i I ) ϑ l,i 1.75, l = 1, 2,

18 18 A. RÖLLIN where 6.5) ϑ l,i = 1 2 ξ i Z 2 i D l L W i ξ i, Z i ) )} k K i ξ i Z ik V ik D l L W ik ξ i, Z ik, V ik ) )} ξi Z ik } Zi V ik D l L W ik Z i, V ik ) )}. k K i Acknowledgments. I thank A. D. Barbour, D. Schuhmacher and B. Nietlispach for many helpful discussions. REFERENCES [1] A. D. Barbour 1988). Stein s Method and Poisson Process Convergence. J. Appl. Probab. 25A, [2] A. D. Barbour and J. L. Jensen 1989). Local and Tail Approximation Near the Poisson Limit. Scand. J. Statist. 16, [3] A. D. Barbour and A. Xia 21). The number of two-dimensional maxima. Adv. in Appl. Probab. 33, [4] A.D. Barbour and V. Čekanavičius 22). Total variation asymptotics for sums of independent integer random variables. Ann. Probab. 3, [5] A.D. Barbour and A. Xia 1999). Poisson perturbations. ESAIM Probab. Statist. 3, [6] A.D. Barbour, M. Karoński, and A. Ruciński 1989). A central limit theorem for decomposable random variables with applications to random graphs. J. Combin. Theory Ser. B 47, [7] A.D. Barbour, L. Holst, and S. Janson 1992). Poisson Approximation. Oxford University Press. [8] T. C. Brown and A. Xia 21). Stein s method and birth-death processes. Ann. Probab. 29, [9] V. Čekanavičius 1997). Asymptotic expansions in the exponent: a compound Poisson approach. Adv. in Appl. Probab. 29, [1] V. Čekanavičius 1998). Poisson approximations for sequences of random variables. Statist. Probab. Lett. 39, [11] V. Čekanavičius and P. Vaǐtkus 21). Centered Poisson approximation by the Stein method. Lithuanian Math. J. 41, [12] L. H. Y. Chen and Q.-M. Shao 24). Normal approximation under local dependence. Ann. Probab. 32, [13] D. J. Daley and D. Vere-Jones 1988). An introduction to the theory of point processes. Springer Series in Statistics. Springer-Verlag, New York. [14] A. Dembo and S. Karlin 1992). Poisson approximations for r-scan processes. Ann. Appl. Probab. 2, [15] A. Dembo and Y. Rinott 1996). Some examples of normal approximations by Stein s method. In Random Discrete Structures, IMA Vol. Math. Appl. 76, Springer, New York.

19 SYMMETRIC AND CENTERED BINOMIAL APPROXIMATION 19 [16] W. Ehm 1991). Binomial approximation to the Possion Binomial distribution. Statist. Probab. Lett. 11, 7 16, [17] L. Le Cam 1965). On the distribution of sums of independent random variables. In Proc. Internat. Res. Sem., Statist. Lab., Univ. California, Berkeley, Calif., Springer- Verlag, New York. [18] B. Matérn 196). Spatial variation: Stochastic models and their application to some problems in forest surveys and other sampling investigations. Meddelanden Fran Statens Skogsforskningsinstitut, Band 49, Nr. 5. Stockholm. [19] G. Reinert 1998). Coupling for normal approximations with Stein s method. DIMACS Ser. Discrete Math. Theoret. Comput. Sci. 41, [2] A. Röllin 25). Approximation of sums of conditionally independent variables by the translated Poisson distribution. Bernoulli 11, [21] B. Roos 23). Kerstan s method for compound Poisson approximation. Ann. Probab. 31, [22] C. Stein 1986). Approximate computation of expectation. IMS Lecture Notes Monogr. Ser. 7. Universität Zürich, Winterthurerstrasse 19, 852 Zürich, Switzerland adrian.roellin@math.unizh.ch