arxiv: v1 [math.pr] 12 Jan 2017

Transcription

1 Central limit theorem for the Horton-Strahler bifurcation ratio of general branch order Ken Yamamoto Department of Physics and arth Sciences, Faculty of Science, University of the Ryukyus, 1 Sembaru, Nishihara, Okinawa 90 01, Japan Abstract The Horton-Strahler ordering method, originating in hydrology, formulates the hierarchical structure of branching patterns using a quantity called the bifurcation ratio The main result of this paper is the central limit theorem for bifurcation ratio of general branch order This is a generalized form of the central limit theorem for the lowest bifurcation ratio, which was previously proved Some useful relations are also derived in the proofs of the main theorems 1 INTRODUCTION arxiv:170101v1 mathpr 1 Jan 017 Branching objects are found very widely 1, ranging from natural patterns like river networks, plants, and dendritic crystals, to conceptual expressions like binary search trees in computer science and phylogenetic trees in taxonomy The topological structure of a branching pattern is modeled by a binary tree if a segment bifurcates does not trifurcate or more) at every branching point Let Ω n denote the set of the different binary trees having n leaves The number of leaves is called the magnitude in research of branching patterns As known well 8, the number of the different binary trees of magnitude n is given by ) 1 n Ω n = n n = n )! n!n)!, which iscalled the n st Catalan number In Fig 1, Ω n for n =,, and are schematically shown Introducing the uniform probability measure P n on Ω n so that each binary tree is assigned equal probability 1/ Ω n ), we obtain the probability space Ω n,p n ) referred to as the random model 6 The formation of real-world branching patterns more or less involves stochastic effects, and the random model is a kind of mathematical simplification of such random factors In hydrology, methods for measuring the hierarchical structure of a river network have been proposed by Horton, Strahler 5, Shreve 6, Tokunaga 7, and other researchers Their methods define how to assign an integer number called the order) to each stream Among all, Strahler s method is currently the most popular because of its simple computation rule Strahler s method is a refinement of Horton s method, so it is sometimes called the Horton-Strahler ordering method The Horton-Strahler method recursively defines the order of each node by the following rules i) The leaf nodes are defined to have order one ii) A node whose children have different order r 1 and r r 1 r ) has order max{r 1,r } iii) A node whose two children have the same order r has order r1 We define a branch of order r as a maximal connected path made by nodes of equal order r A branch here is called a stream in the analysis of river networks) An example of Strahler s ordering is shown in Fig For a binary tree τ Ω n, we let S r,n τ) denote the number of branches of order r in τ By the definition of the order, S 1,n τ) = n and 0 S r,n τ) n/ r r ) Note that S,n τ) 0 if n, because a node of order is produced by the merge of two leaves For the binary tree τ Ω 6 ) in Fig, S 1,6 τ) = 6, S,6 τ) =, S,6 τ) = 1, and S r,6 τ) = 0 for r S r,n is a random variable on Ω n,p n ), and its stochastic property is of main interest in this study For any function f : {0,1,,} R, fs r,n )) is a real-valued random variable on Ω n According to Ref 9, the Ω Ω Ω FIG 1: Ω,Ω, and Ω contain one, two, and five binary trees, respectively

2 branch of order 1 branch of order branch of order 1 1 FIG : A small example of ordering and branches The number on each node represents the order of the node The branches of order and are shown by the dashed rectangles This binary tree consists of six branches of order 1, two branches order, and one branch of order recursive relation between the averages of the rth and r st variables fs r,n ) = n!n)!n )! n )! n/ holds, where denotes the average on the random model The coefficient n!n)!n )! n m n )!n m)!m!m)! n m n m)!m!m)! fs r,m) 1) represents the probability P n S,n = m) In particular, putting r = in q 1), we have fs,n ) = n!n)!n )! n )! n/ n m fm) ) n m)!m!m)! Mathematical properties of S,n have been investigated thoroughly For instance, the average and variance are respectively given by 10 S,n = nn) n ), VarS,n = nn)n )n ) n ) ) n 5) Moreover, from q ), the moment generating function M,n t) of S,n is given by M,n t) := exps,n t) = n!n)!n )! n )! n/ n m n m)!m!m)! emt, and this summation can be expressed using the Gauss hypergeometric function F 11: M,n t) = n n!n)! n e t F, n ),;e t ) n )! The ratio S r1,n τ)/s r,n τ) is called the bifurcation ratio oforderr or simply the rth bifurcation ratio Hydrologists have empirically confirmed that the bifurcation ratios of an actual river network become almost constant for different orders, and this relation is referred to as Horton s law of stream numbers By definition, the bifurcation ratio is always smaller than or equal to 1/ When S r,n τ) = 0, we reasonably define S r1,n τ)/s r,n τ) = 0 The random variable S r1,n /S r,n is also called the rthe bifurcation ratio The lowest bifurcation ratio S,n /S 1,n = S,n /n is relatively easy to deal with, because it is similar to S,n The central limit theorem for S,n /n has been shown by Wang and Waymire 1: Theorem 1 Central limit theorem for the lowest bifurcation ratio) On the random model, S,n n n 1 ) N 0, 1 ), n, 16 where denotes convergence in distribution, and Nµ,σ ) is the normal distribution with mean µ and variance σ

3 It is a simple and natural idea that we extend Theorem 1 to general order r Compared with S,n, however, higherorder branches S r,n for r and the bifurcation ratio of order r is difficult to handle and less studied In this paper, we generalize Theorem 1 in two ways Theorems and in ), and further generalize them Theorem in 6) In 5, we give proofs of lemmas, which are necessary for the main theorems In these proofs, q 1) and its variant are very useful f Sr1,n S r,n ) = n!n)!n )! n )! n/ n m ) n m)!m!m)! Sr,m f 5) S r,m MAIN RSULTS The following two theorems are the main results of the present paper Theorem Central limit theorem for the bifurcation ratio of general order) For any order r = 1,,,, the rth bifurcation ratio S r1,n /S r,n satisfies Sr1,n n 1 ) N 0, r ), n, 6) S r,n Theorem Central limit theorem for the number of branches of general order) For any order r = 1,,,, the number S r1,n of r1st branches satisfies n Sr1,n n 1 r ) N 0, 1 r ) 16 r, n, 7) Remark These two theorems are generalization of Theorem 1 to general order r; they are reduced to Theorem 1 by setting r = 1 Theorem states the property of the bifurcation ratio S r1,n /S r,n, and Theorem states the property of the number of branches S r1,n The limit variance r in Theorem becomes large as r increases, whereas the limit variance in Theorem becomes small as r increases From Theorem, the following property, which can be regarded as Horton s law of stream numbers, is easily derived Corollary 1 Horton s law of stream numbers for the random model) For any order r = 1,,, the rth bifurcation ratio S r1,n /S r,n converges in probability to the common value 1/: where p denotes convergence in probability S r1,n S r,n p 1, Let us introduce the asymptotic equality, since this study mainly focuses on the asymptotic behavior the limit n ) of S r,n Definition 1 The average value fs r,n ) is asymptotically equivalent to g r n) if and this is denoted by For example, from q ), fs r,n ) lim = 1, n g r n) fs r,n ) g r n) S,n = nn) n ) n Theorems and are easily proved by using the following Lemmas 1 and, respectively

4 Lemma 1 For r = 1,, and s = 0,1,,, Sr1,n 1 ) s s)!! n ) s, S r,n s r Sr1,n Lemma For r = 1,, and s = 0,1,,, S r1,n n ) s r s)!! r ) s sr n s, S r1,n n ) s1 r s1)!! r s1)r S r,n 1 ) s 1 r1 5 The odd-power result has a more complicated form than the even-power one ) s1 s1)!! s1 n r ) s ) r s1) n s Proof of Theorem We let ϕ r1 r,n z) denote the characteristic function of the left-hand side of q 6), where the subscript r and superscript r1 respectively correspond to S r,n in the denominator and S r1,n in the numerator in q 6) By definition, ϕ r1 r,n z) is calculated as ϕ r1 r,n exp z) = iz Sr1,n n 1 )) S r,n iz n) k Sr1,n = 1 ) k k! S r,n k=0 iz Sr1,n n) k ) k = = k=0 k! iz n) s s)! 1 S r,n ) s Sr1,n S r,n 1 iz Sr1,n n) s1 1 ) s1, s1)! S r,n where i = At the last equality, we have split the sum into even k k = s) and odd k k = s1) By Lemma 1, the terms of the first sum even k) are On 0 ), whereas the terms of the second sum odd k) are on 0 ) Hence, the second sum can be neglected in the limit n, so that ϕ r1 iz Sr1,n n) s r,n z) 1 ) s s)! S r,n z n) s s)!! n s s)! s r)s = r z ) s 1 s! = exp r z ) Recall that the characteristic function of Nµ,σ ) is expiµz σ z /) Since ϕ r1 r,n converges pointwise to the characteristic function of N0, r ), convergence in distribution in Theorem is proved For the properties of a characteristic function, see Feller 1 for example) Keep in mind that the neglect of the odd-power terms is a crucial point also in the other central limit theorems in this paper Proof of Theorem As with the above proof of Theorem, the characteristic function ϕ r1 1,n z) of the left-hand side

5 5 of q 7) is ϕ r1 1,n exp z) = iz Sr1,n n 1 )) n r iz Sr1,n n) k = 1 ) k k! n r k=0 iz n) s = s)!n s S r1,n n ) s r iz n) s1 s1)!n s1 S r1,n n ) s1 r By Lemma, the second sum is neglected and the dominant terms are calculated to ϕ r1 1,n z) iz n) s s)!n s S r1,n n ) s r z n) s s)!! r ) s s)!n s rs n s = z r ) s 1 16 r s! = exp z r ) 16 r Therefore, the converges in distribution to N0, r )/ 16 r )) is proved We give the proofs of Lemmas 1 and in the following three sections STARTING POINT OF LMMAS 1 AND We show Lemmas 1 and by induction on r In this section, the case of r = 1 in Lemmas 1 and is proved Cor ) Proposition 1 For a two-variable polynomial p, ) of finite degree, S,n ps,n,n) = n ps,n,n) nn ) n ) ps,n,n) 8) Here, S,n ps,n,n) and ps,n,n) are taken over Ω n, whereas ps,n,n) are over Ω n Proof Because of the linearity of, it is sufficient to check the case ps,n,n) = S,nn k l The average is expressed using the moment generating function M,n t) in q ) S,n S,nn k l = n l S,n k1 = n l d k1 dt k1m,nt) t=0 = n ln n!n)! d k1 n n )! dt k1et F, n ),;e t t=0 = n ln n!n)! d k n n )! dt ket F, n ),;e t d n dt F, n ),;e t t=0 Using the derivative of the Gauss hypergeometric function 1 d dz F α,β,γ;z) = α F α1,β,γ;z) F α,β,γ;z) z and the symmetry F α,β,γ;z) = F β,α,γ;z), we have d n dt F, n ),;e t = n n F, n n,;e ) F t, n ),;e t, 9)

6 6 thereby S,n S,n k nl = n ln n!n)! d k n n n )! dt ket F n = n l S,n k nn) n )n ) = n S k,n nl nn ) n ) S k,n nl, n ),;e t n S,n k n n F ), n ),;e t Here we comment on the utility of Prop 1 We can calculate the moments S k,n recursively using q 8): S,n = S,n 1 = n nn ) nn) 1 1 = n ) n ), S,n = S,n S,n = n S,n nn) n ) S,n = nn)n n ), n )n 5) S,n = S,n S,n n = S,n nn) n ) S,n nn)n n 5n n8) =, 8n )n 5)n 7) and so on The first and second moments were individually calculated 10 see q )) Note, however, that Prop 1 provides a systematic calculation method of the kth moment of S,n in a bottom-up way Furthermore, we easily obtain S k,n n/) k for k = 0,1,, using q 8) Corollary Subtracting nps,n,n)/ from q 8), we have S,n n ) ps,n,n) = n ps,n,n) nn ) n ) ps,n,n) Lemma For k = 0,1,,, S,n n ) k k! k1)/ k/! kn k/ That is to say, the asymptotic form of S,n n/) k depends on whether k is even or odd: for s = 0,1,, S,n n ) s s)!! s n s, S,n n ) s1 s1)!! s1 ns, 10) Remark In S,n n/ k = 1), the leading On 1 ) terms are canceled because of S,n n/ Similarly, in general k, On k ),On k ), terms are successively canceled, so that the resultant leading order of S,n n/) k becomes k/ This effect makes the estimation of S,n n/) k difficult Proof The proof is by induction on k The statement is true for k = 0,1, and, because S,n n ) 0 = 1, S,n n ) 1 = nn) n ) n = n n ) 1 8, S,n n ) = S n,n S,n n 16 = nn 7n16) 16n )n 5) n 16

7 Assume that it is true up to k, and we show it is true for k 1 It follows from Cor that S,n n ) k1 = S,n n ) S,n n ) k = n S,n n ) k nn ) n ) S,n n ) k = n S,n n ) k nn ) n ) S,n n 1 ) k = n S,n n ) k nn ) n ) k ) k 1 l l ) S,n n ) k l In the following, we separately investigate odd and even k Case 1: k is odd k = s 1) The dominant terms are l = 0 and 1 in the summation, and the others can be neglected Thus, S,n n ) s11 n S,n n n ) s1 nn ) S,n n ) s1 s1 S,n n ) ) s n ) s1 n) s s1 ) s)!! s n) s s1)!! s1 n s1)!! s1)!! s n s1 Case : k is even k = s) Picking out the terms up to On s ) carefully, we obtain S,n n ) s1 = n S,n n ) s nn ) S,n n ) s s n ) S,n n ) s ss) S,n n ) ) s on s ) 11) 16 We introduce the coefficient a s as S,n n ) s = s)!! s n s a s n s on s ), and expand each term on the right-hand side of q 11) up to On s ) using the induction hypothesis: n S,n n ) s = s)!! s1 n s1 a s ns on s ), nn ) n ) S,n n ) s = s)!! s1 n s1 s1)!! s1 ns a s ns on s ), nn ) n ) nn ) s n ) ss) 16 S,n n ) s S,n n ) s = s)!! s1 sn s on s ), = s)!! s1 sn s on s ) 7 Therefore, S,n n ) s1 = s1)!! s1 ns on s )

8 8 Note that On s1 ) terms and those including a s are all cancelled Therefore, the statement is true for any k Corollary Multiplying q 10) by n k, we have S,n n 1 ) s s)!! s n s, S,n n 1 Remark This result corresponds to the case r = 1 in Lemmas 1 and ) s1 s1)!! s1 n s Using this corollary, we can provide another proof of Theorem 1 as follows The characteristic function ϕ 1,n z) for the lowest bifurcation ratio in Theorem 1 is calculated as ϕ 1,n exp z) = iz S,n n n 1 )) iz n) k S,n = k! n 1 ) k k=0 iz S,n n) k ) k = k=0 k! iz n) s = s)! ) z s 1 s! ) = exp z, n 1 ) s S,n n 1 so the convergence of ns,n /n/) to N0,1/16) is proved iz S,n n) s1 s1)! n 1 ) s1 PROOF OF LMMA 1 We first derive the asymptotic form of S k,n, which is needed in the proof of Lemma 1 Proposition For k = 0,1,,, we have S,n k n ) k Remark Since S,n τ) 0 for any τ Ω n, S k,n surely takes a finite value and S k,n is not divergent Proof Let us introduce the operator d/dt) defined by ) d t ft) := fs)ds, dt where f is integrable on any interval,t) Note that d/dt) is the inverse of d/dt Owing to the property ) k d e mt = m k e mt, dt

9 9 the average of S,n k is expressed by n/ S,n k n!n)!n )! = n )! = n!n)!n )! n )! = n n!n)! n )! d dt n/ n m n m)!m!m)! m k n m n m)!m!m)! ) k n e t F Let us derive a relation between S,n k and S k1),n as in Prop 1 S k1),n = n n!n)! n )! = n n!n)! n )! = n n!n)! n )! = n n!n)! n )! ) k ) d d n e t F dt dt ) k d t n e s F, n dt ) { k d n e s F dt ) { k d n e t F dt = S,n k n n!n)! d n )! dt, n ),;e t t=0,;e )ds s ) k d e mt dt t=0, n ),;e t t=0 t=0, n ) t t,;e s s=, n ) d,;e ) t e t d dt dt F ) k1) e t d n dt F By using the derivative of the hypergeometric function in q 9), whose asymptotic form is or n S k1),n S k1),n Considering S 0,n = 1, we obtain = S k,n n = S,n k nn ) n ) S k1),n S k1),n S k1),n, n ),;e t t=0 nn ) n ) n S k,n S,n k n ) k e s d n ds F n S k1),n S,n k n )ds}, n,;e s t=0, n ) },;e t t=0, S k1),n, Proof of Lemma 1 By induction on r For r = 1, the statement is equivalent to Cor Assume that it is true up to r, and we show it is true for r Using q 5), Sr1,n S r,n 1 ) k = n!n)!n )! n )! n/ n m n m)!m!m)! Sr,m 1 ) k S r,m

10 10 Case 1: k is even k = s) Sr1,n S r,n 1 ) s = n!n)!n )! n )! n!n)!n )! n )! n/ n/ n m n m)!m!m)! Sr,m 1 ) s S r,m n m s)!! m ) s n m)!m!m)! s r = s)!! n/ n!n)!n )! n m s sr ) n )! n m)!m!m)! m s By using q 5) again and Prop, Sr1,n 1 ) s s)!! S r,n s sr ) S,n s s)!! n ) s s)!! n ) s = s sr ) s r Case : k is odd k = s1) Using q 5) and Prop as above, Sr1,n S r,n 1 ) s1 = n!n)!n )! n )! n!n)!n )! n )! n/ n/ s1)!! s1 r) s) S,n s = s1)!! s1 n r ) s n m n m)!m!m)! Sr,m 1 ) s1 S r,m n m s1)!! m ) s n m)!m!m)! s1 r 5 PROOF OF LMMA In the proof of Lemma, we use the following relations Lemma For s,l = 0,1,,, S,n l S,n n ) s n ) l s)!! s n s, S,n l S,n n ) s1 n ) l s1)!! s1 l1)ns Remark The complicated form of the odd-power result in Lemma is actually due to the factor l1) Proof By induction on l For l = 0, the statement is equivalent to Lemma Assume that it is true up to l 0, and we show for l1 Using Prop 1, S,n l1 S,n n ) k = n S,n l S,n n ) k nn ) n ) S,n l = n S,n l S,n n ) k nn ) n ) S,n l = n S l,n S,n n ) k nn ) n ) k p=0 S,n n ) k S,n n ) k 1 p S,n p ) l Case 1: k is even k = s) In the summation, only p = 0 is dominant, so that S,n l1 S,n n ) s n S,n l S,n n ) s nn ) n ) S l,n ) k 1 S,n n ) k p S,n n ) s

11 11 By using the induction hypothesis, S l1,n S,n n ) s n n ) l s)!! s n s n n ) l s)!! n ) l1 s)!! s n) s s n s Case : k is odd k = s1) Note that p = 0 and 1 in the summation are the dominant terms S,n l1 S,n n ) s1 n S,n l S,n n ) s1 nn ) S,n l S,n n ) s1 s1) S,n l S,n n ) ) s n ) n n ) l s1)!! l1)ns s1 n nn ) ) l s1)!! n ) s1 l1)n)s s1) ) l n s)!! s n) s n ) l1 s1)!! s1 l)ns Proof of Lemma By induction on r For r = 1, the statement is equivalent to Lemma Assume that it is true up to r, and we show it is true for r Using q 5) to calculate S r1,n n ) k r = n!n)!n )! n/ n )! = n!n)!n )! n )! = n!n)!n )! n )! n/ n/ n m n m)!m!m)! S r,m n ) k r n m n m)!m!m)! S r,m m r m r n ) k r n m n m)!m!m)! k1 ) k m l r n ) k l r S rm m ) l r We split the summation over l according to the parity of l, and use the induction hypothesis The estimation of the summation is complex compared with the other proofs above

12 1 Case 1: k is even k = s) S r1,n n ) s r = n!n)!n )! n )! s s l n/ ) m r n r ) s l s ) s m l1 n!n)!n )! n )! s ) s 1 l r n m n m)!m!m)! r n r ) s l n/ S rm m r ) l n m n m)!m!m)! S r,m m ) l1 r ) s l m n ) s l l)!! r ) l m l lr) s ) ) s l s 1 l1 r m n ) s l l1)!! r ) l1 m l l1)r) s ) ) s l s 1 l)!! r ) l = l r S lr),n n ) s ls l,n s ) ) s l s 1 l1)!! r ) l1 l1 r S l1)r),n n ) s ls l,n From Lemma, the first summation is On s ), while the second summation is On s ) Thus, we can neglect the second sum, so that S r1,n n ) s r Case : k is odd k = s1) s = s)!ns s! s sr ) s 1 l s s = s)!!ns sr = s)!!ns sr S r1,n n ) s1 r = n!n)!n )! n )! s s1 s r ) s l l)!! r ) r l 1 r r ) s n/ lr) ) l ) s n m n m)!m!m)! ) m r n ) s1 l r S rm m ) l r l ) s1 m l1 r n ) s l r S r,m m r ) l1 ) l 1 l s l)!! s l n s

13 1 n/ n m n!n)!n )! n )! n m)!m!m)! s ) ) s1 l s1 1 l r m n ) s1 l l)!! r ) l m l lr) s ) ) s l s1 1 l1 r m n ) s l l1)!! r ) l l1)r) 1 r r l1) )m l 5 s ) ) s1 l s1 1 l)!! r ) l = l r S lr),n n ) s1 ls l,n s ) ) s l s1 1 l1)!! r ) l l1 r l1)r) 1 r ) r l1) S,n n ) s ls l 5,n Both two sums are On s l ), so we need to consider them By using Lemma, S r1,n n ) s1 r s ) ) s1 l s1 1 l)!! r ) l n ) l s l1)!! l r lr) l1)ns l s l1 s ) ) s l s1 1 l1)!! r ) l l1 r l1)r) 1 r ) r n ) l s l)!! l1) 5 s l n s l s ) = s1)!!ns s r ) l l1) s ) s r ) l r rs1) l 5 l s = s1)!!ns r rs1) 15 r = s1)!!ns rs1) ) s r ) l l5) l r1 ) s 1 5 ) r s1), ) r l1) where we have used s ) s r ) l r ) s =, l s ) s r ) l r ) s r l = s l 6 FURTHR GNRALIZATION Theorems and are further generalized to the central limit theorem as follows Theorem For q,r = 1,,, Sqr,n n 1 ) r ) S q,n r N 0, r q1 Remark This theorem is reduced to Theorem when r = 1 and to Theorem when q = 1; moreover, it is reduced to Theorem 1 when r = q = 1

14 1 Lemma 5 For q,r = 1,, and s = 0,1,,, Sqr,n 1 ) s S q,n r sq)s)!! r ) s sr n s, Sqr,n 1 ) s1 S q,n r sq)s1)!! r s1)r ) s 1 r1 5 ) r s1) n s Remark In comparison with Cor, the effect of q > 1 appears in the form of the factor sq) Proof By induction on q q = 1 is equivalent to Cor Assume that the statement is true for q 1, and we show that it is true for q 1 Using q 5) and Prop, Sq1r,n S q1,n 1 r ) s = n!n)!n )! n )! n/ Sqr,n n m n m)!m!m)! 1 ) s S q,n r n!n)!n )! n/ n m r ) s n )! n m)!m!m)! sq)s)!! sr m s = sq)s)!! r ) s sr S s,n sqs)!! r ) s sr n s Similarly, Sq1r,n S q1,n 1 r ) s1 sqs1)!! s1)r r ) s 1 r1 5 ) r s1) n s Proof of Theorem By using Lemma 5, the characteristic function of the left-hand side, ϕq,n qr z), can be calculated as with Theorems and : ϕq,n qr z) exp 1 r ) r q1z 1 P Ball, Branches Oxford University Press, Oxford, 011) D Knuth, The Art of Computer Programming, vol Addison Wesley, Reading, 197) J D Archibald, Aristotle s Ladder, Darwin s Tree: The volution of Visual Metaphors for Biological Order Columbia University Press, New York, 01) R Horton, Geol Soc Am Bull 56, ) 5 A N Strahler, Trans Am Geophys Un 8, ) 6 R L Shreve, J Geol 75, ) 7 Tokunaga, Geogr Rep Tokyo Metrop Univ 1, ) 8 R P Stanley, numerative Combinatorics, vol Cambridge University Press, Cambridge, 1999) 9 K Yamamoto, J Stat Phys 19, 6 010) 10 C Werner, Canadian Geographer 16, ) 11 K Yamamoto, Phys Rev 78, ) 1 S X Wang and C Waymire, SIAM J Discr Math, ) 1 W Feller, An Introduction to Probability Theory and Its Applications, vol Wiley, New York, 1968) 1 M Abramowitz and I A Stegun, Handbook of Mathematical Functions Dover, New York, 197)