Approximating the negative moments of the Poisson distribution
|
|
- Arron Johns
- 5 years ago
- Views:
Transcription
1 Approximating the negative moments of the Poisson distribution C. Matthew Jones, 2 and Anatoly A. Zhigljavsky School of Mathematics, Cardiff University, CF24 4YH, UK. 2 Cardiff Research Consortium, Cardiff Medicentre, CF4 4UJ, UK. Abstract Through the use of the Stirling numbers of the first kind, a general class of approximations is derived for the negative moments of the Poisson distribution. Keywords: Negative moments; Stirling numbers of the first kind; Poisson distribution Introduction The approximation of negative moments of the Poisson distribution is an important practical problem and has attracted significant attention in literature, see for example Chao and Strawderman (972), Gupta (979), Stancu (968), Stephan (945) and Tiku (964). The main applications of the first negative moment are related to the theory of mixed Poisson distributions: indeed, if η i (i =, 2,... n) are i.i.d.r.v. with variance σ 2 and the sample size n has the positive Poisson distribution, then the variance of the sample mean ni= η n i is σ 2 E (/n), see e.g. Grab and Savage (954). The negative moments of the Poisson distribution are known to be useful in life testing problems, as shown in Bartholomew (957), David and Johnson (956) and Epstein and Sobel (954). A recent application has occurred in the field of multi-centre clinical trials, which involves the simultaneous recruitment of patients to numerous centres. The approximation of negative moments are required in the subsequent calculations of M SE s of the mean treatment effect differences. As demonstrated in Dragalin et al. (2002), three least squared estimators can be defined for the mean treatment difference based on fixed effect models of increasing address: matthew.jones@cardiff-research-consortium.co.uk address: ZhigljavskyAA@cardiff.ac.uk
2 complexity. In particular, the MSE of the Type III estimator is σ 2 N 2 N i= ( + ), () n i2 n i where N is the number of centres, σ 2 is the variance of the response observed at each centre, and n ij is the number of patients of treatment j at centre i. If n ij are assumed to be positive Poisson random variables then () becomes a random variable, with moments that are expressed through the negative moments of the Poisson distribution. Similar phenomena occurs with the Type I and II estimators, however, with a more complicated form. Fedorov et al. (2003) for further details. Let Poisson(λ) denote the Poisson distribution with parameter λ, and let ξ be a random variable, ξ Poisson(λ). Additionally, let ξ + parameter λ; that is, Pr(ξ + = x) = ξ α + See be the so-called positive Poisson r.v. with e λ λ x x! e λ, for x =, 2... The negative moments of the Poisson distribution are defined as the negative moments of ξ + : ( ) µ α = E = e λ λ x for any α > 0. (2) e λ x α x! We develop the following k th order approximations for µ α, with integer α: µ α = µ (k) α = x= k u=α s(u, α) λ u, (3) where s(u, α) are the so-called signless Stirling numbers of the first kind, see Section 3.. The case of non-integer α is considered in Section 3.6. Approximations of the negative moments µ α introduced in Tiku (964) have been cited in a number of reference books; see, for instance, Haight (966), Johnson et al. (992). We shall use these approximations as a benchmark for comparison in subsequent calculations, and we will demonstrate that our approximations are more precise. Tiku s approximations for the first negative moment (see Equation () in Tiku (964)) are defined as follows: ( ) µ () + j r=3 β r () = T j = ; with T (λ )( e λ 0 = ) (λ )( e λ ), (4) 2
3 where β r = a (r) /λ(λ + )(λ + 2)... (λ + r ), r = 3, 4, 5,... with a (3) =, a (4) = 7, a (5) = 43, a (6) = 27, etc. Note that a number of alternative approximations to µ can be constructed through using the fact that µ = e λ (Ei(λ) log λ γ), (5) e λ where γ is Euler s constant and E i (λ) = λ (eu /u) du is the exponential integral, see for example Grab and Savage (954). Note that the relation (5) is often used to approximate values of the exponential integral Ei(λ); consequently, this broadens the scope of our approximation into several other fields. The approximations of Tiku (964) for the higher order integer negative moments (α 2) are defined as µ α = T (α) = (λ )(λ 2)... (λ α). (6) The suggested approximations are far more precise than Tiku s approximations. The criterion by which we determine the accuracy of the approximations is through studying their relative error: Relative Error = Exact V alue Approximate V alue Exact V alue. Note that positive values for the relative error indicate that the respective approximation is underestimating the true value. In this paper we are only interested in moderate and large values of λ. For small λ the moments can be easily computed using the definition (2). In particular, the relative error of the following simple approximation µ α = e λ e λ k= λ k k α k! = e λ e λ 3λ+0 k= λ k k α k! ; is smaller than 0 0 in absolute value for all λ > 0 and α. 3
4 2 Deriving recurrence relations and approximations 2. General Recurrences Denote In particular, A l,m (λ) = µ α = For each l, m > 0, we obviously have x l = x= e λ λ x x l (x + m )!. (7) e λ A α,(λ). (8) x l (x + m) + Through applying Equation (9) to (7) we obtain the basic recurrence: A l, m (λ) = x= and from (7) we have m x l (x + m). (9) e λ λ x x l (x + m)! + me λ λ x x l (x + m)! = A l, m+(λ) + ma l, m+ (λ), (0) x= A 0,m (λ) = ( m λ m x=0 Note that for any fixed m and λ we have e λ λ x ). () x! λ m ( m λ m e λ) A 0,m (λ) λ m (2) so that the difference A 0,m (λ) /λ m is exponentially small as λ. Through performing several iterations of the basic recurrence (0), we obtain A l, m (λ) = A l,m+ (λ) + ma l,m+2 (λ) + m(m + )A l,m+2 (λ) =... N (m + r )! (m + N)! = A l,m+r+ (λ) + (m )! (m )! A l,m+n+(λ). (3) Rather than using the definition (7), we can define A l,m (λ) through the recurrences (3) and the initial conditions (). Let the quantities B l,m (λ), which can be considered as a simplified version of A l,m (λ), be defined through the recurrences (3), that is, B l,m (λ) = N (m + r )! (m + N)! B l,m+r+ (λ) + (m )! (m )! B l,m+n+(λ) for all N > 0, (4) 4
5 and the initial conditions that are similar to () but exclude the exponential terms: B 0,m (λ) =. (5) λm Deriving approximations using B l,m (λ) as opposed to A l,m (λ) implies that we ignore all the terms that are exponentially decreasing when λ. Since B 0,m (λ) = λ m all m, see (2), we have B l,m (λ) A l,m (λ) for all l 0, m. A 0,m (λ) for 2.2 Deriving approximations and determining their precision Let m, l and k be fixed. Consider the expansions of A l,m (λ) and B l,m (λ) of the form: A l,m (λ) = B l,m (λ) = k k a l,m (r) λ r b l,m (r) λ r where a l,m (r) and b l,m (r) are some coefficients. summarized in the following statement. ( ) + o, λ ; (6) λ k ( ) + o, λ, (7) λ k The properties of these coefficients are Theorem. Let l, k 0, m be some integers and let A l,m (λ) and B l,m (λ) be defined through (6) and (7), respectively. Then we have (i) a l,m (r) = b l,m (r) for all l, m, r; (ii) a l,m (r) = 0 for all r < l + m ; (iii) a l,m (r) = for r = l + m ; (iv) a l,m (r) 0 for all r. Proof. The statement (i) follows from the fact that for fixed l and m the difference between A l,m (λ) and B l,m (λ) is exponentially small as λ. In order to prove the other three statements, we use induction in l. For l = 0 the statements obviously follow from either (2) or (5). Assume that (ii) - (iv) hold for l. For all fixed m, l we have for all x : (x + m)... (x + m + l ) x l ( + C ) m, l (x + m)... (x + m + l ) x + m + l 5
6 for some positive constants C m, l ; these constants are C m, l = (m + l + ) ((m + l)!/m! ). Through substituting this into (7) we obtain A 0, m+l (λ) A l, m (λ) A 0, m+l (λ) + C m, l A 0, m+l+ (λ). In view of (2) we obtain as λ this yields (ii) and (iii). ( ) λ + o A m+l λ m+l l, m (λ) λ m+l + C m, l λ m+l ; Assume that (iv) is true for l, through the recurrence relation in (3) it follows that (iv) is true for l; as it implies that each a l, m (r) is a weighted sum of non-negative coefficients a l, u (v) with non-negative weights. In view of (8) and (6), we define the k th order approximations to µ l as µ (k) l = k a l, (r) λ r. (8) The statements (ii) and (iii) of Theorem imply the summation in (8) actually starts at r = l, they also imply that as λ the order of decrease of µ l is λ l ; more precisely: µ l ( ) = o λ l λ l as λ. (9) Statement (iv) yields that for each l, µ l is monotonously increasing as k grows, that is µ (k) l µ (k+) l for all k. Assume that we able to compute the coefficients a l, (r) (which we can, unlike the coefficients a l, m (r) for m > ). Thus, we can compute the k th order approximations µ (k) l with k > l, with Theorem implying that the precision is: ( ) µ (k) l µ l = O λ k+ Formulae (9) and (20) yield that the relative error of µ (k) l µ(k) l µ l = O ( ) λ k+ µ l + o ( ) = λ l λ l O ( λ k l+ ) + o() 6 to µ l as λ. (20) with k > l is: ( ) = O as λ. (2) λ k l+
7 Despite the rate of decrease (as λ ) of the relative error, µ (k) l increases as k grows, we cannot take k arbitrarily large. This is attributable to the constant in O ( λ k l+ ) in (2) rapidly increasing as k grows. This constant, as follows from (2) has the order of at least s(k +, l), which is roughly k!(log k). Numerical results (shown below) demonstrate that one can guarantee extremely precise approximations by taking k up to λ + l. Note also that Theorem implies that to construct a k th order approximation to µ l, we are free to ignore the last terms in (4) and use for any N k (m + l). B l,m (λ) = N (m + r )! B l,m+r+ (λ), (22) (m )! 3 Construction of the approximations In the approximations that follow we will need the so-called signless Stirling numbers of the first kind. 3. Signless Stirling numbers of the first kind Denote the signless Stirling numbers of the first kind as s(n, j), see Pólya et al. (983) and Comtet (974). These numbers satisfy the recurrence s(n, j) = (n )s(n, j) + s(n, j ), n, j, with the initial conditions s(n, j) = 0 if n 0 or j 0, except s(0, 0) =. Values for these numbers can be obtained from numerous combinatorial books, for example Pólya et al. (983). We shall need the following properties of these Stirling numbers, s(n, j). Let H n = be the harmonic numbers, then (see Comtet (974), p. 27) s(n, j) = 0 2 n for n < j and for all n j s(n +, ) = n! ; (23) 7
8 ( s(n +, 2) = n! ) = n!h n ; (24) n s(n +, 3) = n! 2 [ Hn 2 ( )] 2 n 2. (25) Additionally, for n fixed and variable k, s(n+, k+) are the elementary symmetric functions of the first n integers; that is, for p =, 2,... n we have s(n +, n + p) = i <i 2 <...<i p n i i 2... i p. (26) 3.2 Approximating the first negative moment Through using (4) and (5) we have B, (λ) = λ + λ + 2 k 2 λ +... = r! 3 λ + k!b,k+(λ). (27) r+ To construct k th order approximations for µ, we keep the first k terms in Equation (27), ignoring B,k+ (λ), respectively. In this way we obtain µ = B, (λ) = µ (k) k r! = λ = r+ λ + λ + 2 (k )! (28) 2 λ3 λ k Through using (23) it is evident that we can write µ (k) in the form of (3). Table I: Comparison of the relative errors of the approximations for the first negative moment, against different values of λ. λ µ (5) µ ( 0.5λ ) Relative Error µ ( λ ) T () 0 T () 3 T () The accuracy of the approximations used are shown in Table I, which gives the relative error against different values of λ; all of our approximations compare favorably to the more 8
9 complex Tiku approximations T () j defined in Equation (4). Already the simple approximation µ (5) = λ + λ λ λ λ 5, provides very good precision. The approximations µ (3) and µ (4) are also accurate. Table I shows that the approximation µ ( 0.5λ ) is very accurate for λ 0, and always underestimates µ. On the other hand, µ (k) with k = λ always overestimates µ, and it is even more precise than µ (k) with k = λ/2. This yields that any approximation µ (k), with λ k λ provides an extremely accurate approximation Approximating the second negative moment Using Equation (22) with l = 2 and m =, we obtain µ 2 N = r!b,r+2 (λ). (29) We will be constructing approximations of the form µ 2 = µ (k) 2 = k r= a 2, (r) λ r, which we shall obtain from (29) with N large enough (any N > k would suffice), by the repeated application of the approximations in (22), and keeping all the terms that contribute to a 2, (r) s with r k. Table II: Comparison of the relative errors for the approximations of the second negative moment, against different values of λ. λ µ (6) 2 Relative Error µ ( 0.5λ +) 2 µ ( λ +) 2 T (2) Rewrite (22) as B l,m (λ) = r 0 (m + r )! B l,m+r+ (λ), (30) (m )! 9
10 and, respectively, (29) as Hence, from (3) and (30) we obtain µ 2 = r 0 r!b,r+2 (λ) = r 0 r! s 0 µ 2 = r!b,r+2 (λ). (3) r 0 (s + r + )! λ s+r+2 (r + )! Let t = r + s + ; using Equation (24) we obtain µ 2 = r 0 t r+ t! (r + )λ t+ = t t! t λ t+ = r 0 r + = t r + s 0 (s + r + )!. λs+r+2 t!h t λ = s(t, 2). (32) t+ t 2 λ t The k th order approximation to µ 2 is obtained when we keep the first k terms in the right-hand side of (32). In Table II we discover that µ ( 0.5λ +) 2, µ ( λ +) 2 and even the very simple approximation µ (6) 2 compare very favorable against Tiku s approximation T (2). We find that µ λ + 2 overestimates µ 2, whereas µ ( 0.5λ +) 2 understimates, this suggests any approximation µ (k) 2 with λ + k λ + would provide an extremely accurate approximation Approximating the third negative moment To approximate the third negative moment, µ 3, we express B 3, (λ) in a similar way: µ 3 = B3, (λ) = r!b 2,r+2 (λ) = r! (r + s + )! B,r+s+3 (λ) r 0 r 0 s 0 (r + )! = r 0 = r 0 r + r + + s + )! s 0(r t 0 r + s + 2 λ s 0 t 0 Let u = r + s + t + 2 and s = r + s +, so that µ 3 = u 2 u 2 u! λ u+ s =r+ (r + s + t + 2)! λ t+r+s+3 (r + s + 2)! (r + s + t + 2)!. t+r+s+3 u (r + )(s + ). (33) 0
11 Since H 2 u = u u r+ s=0 u 2 = 2 r + we have from (33) and (25) µ 3 = u 2 u 2 s+ = u s=r+ r+ u s=r+ u s+ + [ ] 2 u + r+ r= r r+ s=0 s+ ( s ), (34) 2 u 2 [ u! H 2 2λ u+ u ( )] = s(u+, 3) 2 u 2 u 2 λ u+ = u 3 The k th order approximation µ (k) 3 is obtained by keeping the first k terms in (35). s(u, 3) λ u. (35) Table III: Comparison of the relative errors for approximations of the third negative moment, against different values of λ. λ µ (7) 3 Relative Error µ ( 0.5λ +2) 3 µ ( λ +2) 3 T (3) In Table III we discover that all approximations compare favorably against Tiku s approximation T (3), and if we had and approximation of the form µ (k) 3 with 0.5λ + 2 k λ + 2, the approximation would be extremely accurate. 3.5 General Case For the l th order negative moment we approximate B l, (λ). This can be obtained by applying the same methodology: µ l = µ l = r 0 r + r 2 0 r +r r l 0 r +r r l +l r l 0 (r +r r l +l )! λ r +r r l +l.
12 Let r m = m i= r i + (m ), for 2 r l. Thus, we have µ l = l! r r l l λ r l + r l (l ) r =0 r l (l 2) r 2 =r + r l (l 3) r 3 =r r l (r + )(r 2 + )... (r l + ). r l =r l 2 + Let t = r + and t m = r m +, 2 m l. Hence, using Equation (26) µ l = t l l = t l l t l! λ t l+ t l (l 2) t = s(t l +, l) λ t l+ t l (l ) t l (l 2) t l 2 t l =t l 2... t 2 =t t 3 =t 2 t t 2... t l = s(r, l). (36) r l λ r Similarly to the above, the k th order approximation to µ l is obtained when we keep the first k terms in (36). 3.6 Approximating non-integer negative moments Assume that α > 0 is not an integer. Set l = α and let θ = α l with 0 θ. Accurate approximations for non-integer negative moments can be obtained through using the following approximation: µ (t+α) α = t s(r +, l) + θ [s(r + l +, l + ) s(r +, l)] λ α+r. (37) Table IV clearly demonstrates that all of the approximations for the non-integer negative moments provide very accurate approximations for λ 50, with a relative error less than However, it can also be seen that for larger values of α (e.g. α = 5. and 5.5) the approximations are nor very accurate for small λ. This can be explained by the fact that the approximations for the larger integer negative moments are only accurate for fairly large values of λ. Additionally, we find that the approximations are less accurate as θ departs further away from zero or one. 2
13 Table IV: Comparison of relative errors of the non-integer negative moments for different values of α versus different values of λ, for t = 5. Relative Error λ µ 0. µ 0.5 µ 2. µ 2.5 µ 5. µ Acknowledments We kindly acknowledge the support of this work by GlaxoSmithKline. We are especially grateful to Valerii Fedorov and Byron Jones. We are also very grateful to our Cardiff University colleague Nikolai Leonenko for useful advice. Finally, we thank the referee for useful comments. References Chao M.T. and Strawderman W.E., (972). Negative moments of positive random variables, Journal of the American Statistical Association 67, Comtet, L., (974). Advanced combinatorics: the art of finite and infinite expansions, (D. Reidel Publishing Company, Dordrecht.). David F. N. and Johnson N.L., (956). Reciprocal Bernoulli and Poisson varaibles, Metron 8 (No.-2), Dragalin F., Fedorov V., Jones B. and Rockhold F, (2002). Estimation of the combined response to treatment in multicentre trials, Journal of Biopharmaceutical Statistics, Epstein B. and Sobel M., (954). Some theorems relevant to life testing, Annals of Mathematical Statistics 25, Fedorov V., Jones B., Jones C.M. and Zhigljavsky, (2003). Multicentre clinical trials with random enrolment: theoretical approximations to the moments of the mean squared errors of estimators, GlaxoSmithkline Technical Report, GSK BDS Technical Report Grab E.L. and Savage I.R., (954). Tables of the expected value of /X for positive Bernoulli and Poisson 3
14 variables, Journal of the American Statistical Association 49, Gupta R.C., (979). On negative moments of generalized Poisson distribution, Math. Operationsforsch. Statist. 0, No., Haight F. A., (966). Handbook of the Poisson distribution, (Wiley, N.Y.). Johnson N., Kotz S. and Kemp A., (992). Univariate discrete distributions, (Wiley, N.Y.). Pólya G., Tarjan R. E., Woods. D. R., (983). Notes on introductory combinatorics, (Birkh user, Boston). Stancu D. D.,(968). On the moments of negative order of the positive Bernoulli and Poisson variables, Studia Universitatis Babes Bolyai Series, Math. Phys., Stephan F. F.,(945). The expected value and variance of the reciprocal and other negative powers of a positive Bernoullian variate, Annals of Mathematical Statistics 6, Tiku M. L., (964). A note on the negative moments of a truncated poissson variate, Journal of the American Statistical Association, 59,
arxiv:math/ v1 [math.st] 9 Nov 2005
arxiv:math/511226v1 [math.st] 9 Nov 25 ASYMPTOTIC EXPANSION FOR INVERSE MOMENTS OF BINOMIAL AND POISSON DISTRIBUTIONS MARKO ŽNIDARIČ Abstract. An asymptotic expansion for inverse moments of positive binomial
More informationarxiv: v1 [math.st] 24 Sep 2008
Inverse moments of univariate discrete distributions via the Poisson expansion arxiv:0809.455v [math.st] 24 Sep 2008 Koenraad Audenaert September 24, 2008, : Abstract In this note we present a series expansion
More informationPermutations with Inversions
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 4 2001, Article 01.2.4 Permutations with Inversions Barbara H. Margolius Cleveland State University Cleveland, Ohio 44115 Email address: b.margolius@csuohio.edu
More informationOn the Convergence of the Summation Formulas Constructed by Using a Symbolic Operator Approach
On the Convergence of the Summation Formulas Constructed by Using a Symbolic Operator Approach Tian-Xiao He 1, Leetsch C. Hsu 2, and Peter J.-S. Shiue 3 1 Department of Mathematics and Computer Science
More informationON DIVISIBILITY OF SOME POWER SUMS. Tamás Lengyel Department of Mathematics, Occidental College, 1600 Campus Road, Los Angeles, USA.
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 7 (007, #A4 ON DIVISIBILITY OF SOME POWER SUMS Tamás Lengyel Department of Mathematics, Occidental College, 600 Campus Road, Los Angeles, USA
More informationSimplifying Coefficients in a Family of Ordinary Differential Equations Related to the Generating Function of the Laguerre Polynomials
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Applications and Applied Mathematics: An International Journal (AAM Vol. 13, Issue 2 (December 2018, pp. 750 755 Simplifying Coefficients
More informationA New Identity for Complete Bell Polynomials Based on a Formula of Ramanujan
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 12 (2009, Article 09.3. A New Identity for Complete Bell Polynomials Based on a Formula of Ramanujan Sadek Bouroubi University of Science and Technology
More informationSome Congruences for the Partial Bell Polynomials
3 47 6 3 Journal of Integer Seuences, Vol. 009), Article 09.4. Some Congruences for the Partial Bell Polynomials Miloud Mihoubi University of Science and Technology Houari Boumediene Faculty of Mathematics
More informationON A KIND OF GENERALIZED ARITHMETIC-GEOMETRIC PROGRESSION
ON A KIND OF GENERALIZED ARITHMETIC-GEOMETRIC PROGRESSION Leetsch C. Hsu* Dalian University of Technology, Dalian 116024, China (Submitted August 1995) 1. INTRODUCTION Let (t\h) p denote the generalized
More informationCourse : Algebraic Combinatorics
Course 18.31: Algebraic Combinatorics Lecture Notes # 10 Addendum by Gregg Musiker (Based on Lauren Williams Notes for Math 19 at Harvard) February 5th - 7th, 009 1 Introduction to Partitions A partition
More informationKatz Family of Distributions and Processes
CHAPTER 7 Katz Family of Distributions and Processes 7. Introduction The Poisson distribution and the Negative binomial distribution are the most widely used discrete probability distributions for the
More informationA Note about the Pochhammer Symbol
Mathematica Moravica Vol. 12-1 (2008), 37 42 A Note about the Pochhammer Symbol Aleksandar Petoević Abstract. In this paper we give elementary proofs of the generating functions for the Pochhammer symbol
More informationThe Generating Functions for Pochhammer
The Generating Functions for Pochhammer Symbol { }, n N Aleksandar Petoević University of Novi Sad Teacher Training Faculty, Department of Mathematics Podgorička 4, 25000 Sombor SERBIA and MONTENEGRO Email
More informationEXPONENTIAL POLYNOMIALS, STIRLING NUMBERS, AND EVALUATION OF SOME GAMMA INTEGRALS. Khristo N. Boyadzhiev
January 29, 2010 EXPONENTIAL POLYNOMIALS, STIRLING NUMBERS, AND EVALUATION OF SOME GAMMA INTEGRALS Khristo N. Boyadzhiev Department of Mathematics, Ohio Northern University Ada, Ohio 45810, USA k-boyadzhiev@onu.edu
More informationImpact of Stratified Randomization in Clinical Trials
Impact of Stratified Randomization in Clinical Trials Vladimir V. Anisimov Abstract This paper deals with the analysis of randomization effects in clinical trials. The two randomization schemes most often
More informationModeling Recurrent Events in Panel Data Using Mixed Poisson Models
Modeling Recurrent Events in Panel Data Using Mixed Poisson Models V. Savani and A. Zhigljavsky Abstract This paper reviews the applicability of the mixed Poisson process as a model for recurrent events
More informationOn Certain Sums of Stirling Numbers with Binomial Coefficients
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 18 (2015, Article 15.9.6 On Certain Sums of Stirling Numbers with Binomial Coefficients H. W. Gould Department of Mathematics West Virginia University
More informationExplicit formulas using partitions of integers for numbers defined by recursion
Explicit formulas using partitions of integers for numbers defined by recursion Giuseppe Fera, Vittorino Talamini, arxiv:2.440v2 [math.nt] 27 Feb 203 DCFA, Sezione di Fisica e Matematica, Università di
More informationCongruences for central factorial numbers modulo powers of prime
DOI 10.1186/s40064-016-2028-5 RESEARCH Open Access Congruences for central factorial numbers modulo powers of prime Haiqing Wang * and Guodong Liu *Correspondence: whqlcf@163.com Department of Mathematics,
More informationw = X ^ = ^ ^, (i*i < ix
A SUMMATION FORMULA FOR POWER SERIES USING EULERIAN FRACTIONS* Xinghua Wang Math. Department, Zhejiang University, Hangzhou 310028, China Leetsch C. Hsu Math. Institute, Dalian University of Technology,
More informationMath 320: Real Analysis MWF 1pm, Campion Hall 302 Homework 4 Solutions Please write neatly, and in complete sentences when possible.
Math 320: Real Analysis MWF pm, Campion Hall 302 Homework 4 Solutions Please write neatly, and in complete sentences when possible. Do the following problems from the book: 2.6.3, 2.7.4, 2.7.5, 2.7.2,
More informationON THE MOMENTS OF ITERATED TAIL
ON THE MOMENTS OF ITERATED TAIL RADU PĂLTĂNEA and GHEORGHIŢĂ ZBĂGANU The classical distribution in ruins theory has the property that the sequence of the first moment of the iterated tails is convergent
More informationEVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS
EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS CHARLES HELOU AND JAMES A SELLERS Abstract Motivated by a recent work about finite sequences where the n-th term is bounded by n, we evaluate some classes
More informationGENERALIZATION OF UNIVERSAL PARTITION AND BIPARTITION THEOREMS. Hacène Belbachir 1
#A59 INTEGERS 3 (23) GENERALIZATION OF UNIVERSAL PARTITION AND BIPARTITION THEOREMS Hacène Belbachir USTHB, Faculty of Mathematics, RECITS Laboratory, Algiers, Algeria hbelbachir@usthb.dz, hacenebelbachir@gmail.com
More informationCHAPTER 1 INTRODUCTION
CHAPTER 1 INTRODUCTION 1.0 Discrete distributions in statistical analysis Discrete models play an extremely important role in probability theory and statistics for modeling count data. The use of discrete
More information*!5(/i,*)=t(-i)*-, fty.
ON THE DIVISIBILITY BY 2 OF THE STIRLING NUMBERS OF THE SECOND KIND T* Lengyel Occidental College, 1600 Campus Road., Los Angeles, CA 90041 {Submitted August 1992) 1. INTRODUCTION In this paper we characterize
More informationON SIMILARITIES BETWEEN EXPONENTIAL POLYNOMIALS AND HERMITE POLYNOMIALS
Journal of Applied Mathematics and Computational Mechanics 2013, 12(3), 93-104 ON SIMILARITIES BETWEEN EXPONENTIAL POLYNOMIALS AND HERMITE POLYNOMIALS Edyta Hetmaniok, Mariusz Pleszczyński, Damian Słota,
More informationIDENTITIES INVOLVING BERNOULLI NUMBERS RELATED TO SUMS OF POWERS OF INTEGERS
IDENTITIES INVOLVING ERNOULLI NUMERS RELATED TO SUMS OF POWERS OF INTEGERS PIERLUIGI MAGLI Abstract. Pointing out the relations between integer power s sums and ernoulli and Genocchi polynomials, several
More informationA Symbolic Operator Approach to Several Summation Formulas for Power Series
A Symbolic Operator Approach to Several Summation Formulas for Power Series T. X. He, L. C. Hsu 2, P. J.-S. Shiue 3, and D. C. Torney 4 Department of Mathematics and Computer Science Illinois Wesleyan
More informationa j x j. j=0 The number R (possibly infinite) which Theorem 1 guarantees is called the radius of convergence of the power series.
Lecture 6 Power series A very important class of series to study are the power series. They are interesting in part because they represent functions and in part because they encode their coefficients which
More informationArithmetic properties of lacunary sums of binomial coefficients
Arithmetic properties of lacunary sums of binomial coefficients Tamás Mathematics Department Occidental College 29th Journées Arithmétiques JA2015, July 6-10, 2015 Arithmetic properties of lacunary sums
More informationGeneralized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 16 (2013, Article 1332 Generalized Aiyama-Tanigawa Algorithm for Hypersums of Powers of Integers José Luis Cereceda Distrito Telefónica, Edificio Este
More informationApproximating the Conway-Maxwell-Poisson normalizing constant
Filomat 30:4 016, 953 960 DOI 10.98/FIL1604953S Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Approximating the Conway-Maxwell-Poisson
More informationTakao Komatsu School of Mathematics and Statistics, Wuhan University
Degenerate Bernoulli polynomials and poly-cauchy polynomials Takao Komatsu School of Mathematics and Statistics, Wuhan University 1 Introduction Carlitz [6, 7] defined the degenerate Bernoulli polynomials
More informationA Symbolic Operator Approach to Power Series Transformation-Expansion Formulas
A Symbolic Operator Approach to Power Series Transformation-Expansion Formulas Tian- Xiao He Department of Mathematics and Computer Science Illinois Wesleyan University Bloomington, IL 61702-2900, USA
More informationALL TEXTS BELONG TO OWNERS. Candidate code: glt090 TAKEN FROM
How are Generating Functions used in finding the closed form of sequences involving recurrence relations and in the analysis of probability distributions? Mathematics Extended Essay Word count: 3865 Abstract
More informationOn discrete distributions with gaps having ALM property
ProbStat Forum, Volume 05, April 202, Pages 32 37 ISSN 0974-3235 ProbStat Forum is an e-journal. For details please visit www.probstat.org.in On discrete distributions with gaps having ALM property E.
More informationAlgorithmic Approach to Counting of Certain Types m-ary Partitions
Algorithmic Approach to Counting of Certain Types m-ary Partitions Valentin P. Bakoev Abstract Partitions of integers of the type m n as a sum of powers of m (the so called m-ary partitions) and their
More informationCharacterizations of Levy Distribution via Sub- Independence of the Random Variables and Truncated Moments
Marquette University e-publications@marquette MSCS Faculty Research Publications Mathematics, Statistics Computer Science, Department of 7--5 Characterizations of Levy Distribution via Sub- Independence
More informationNotes 6 : First and second moment methods
Notes 6 : First and second moment methods Math 733-734: Theory of Probability Lecturer: Sebastien Roch References: [Roc, Sections 2.1-2.3]. Recall: THM 6.1 (Markov s inequality) Let X be a non-negative
More informationSolutions to the recurrence relation u n+1 = v n+1 + u n v n in terms of Bell polynomials
Volume 31, N. 2, pp. 245 258, 2012 Copyright 2012 SBMAC ISSN 0101-8205 / ISSN 1807-0302 (Online) www.scielo.br/cam Solutions to the recurrence relation u n+1 = v n+1 + u n v n in terms of Bell polynomials
More informationTHE LOG-BEHAVIOR OF THE SEQUENCE FOR THE PARTIAL SUM OF A LOG-CONVEX SEQUENCE. 1. Introduction
SARAJEVO JOURNAL OF MATHEMATICS Vol.13 (26), No.2, (2017), 163 178 DOI: 10.5644/SJM.13.2.04 THE LOG-BEHAVIOR OF THE SEQUENCE FOR THE PARTIAL SUM OF A LOG-CONVEX SEQUENCE FENG-ZHEN ZHAO Abstract. In this
More informationOn the efficiency of two-stage adaptive designs
On the efficiency of two-stage adaptive designs Björn Bornkamp (Novartis Pharma AG) Based on: Dette, H., Bornkamp, B. and Bretz F. (2010): On the efficiency of adaptive designs www.statistik.tu-dortmund.de/sfb823-dp2010.html
More informationDiscrete Distributions Chapter 6
Discrete Distributions Chapter 6 Negative Binomial Distribution section 6.3 Consider k r, r +,... independent Bernoulli trials with probability of success in one trial being p. Let the random variable
More informationF ( B d > = d I 1 S (d,d-k) F _ k k=0
Canad. Math. Bull. Vol. 26 (3), 1983 ORDERED FIBONACCI PARTITIONS BY HELMUT PRODINGER (1) ABSTRACT. Ordered partitions are enumerated by F n = k fc! S(n, k) where S(n, k) is the Stirling number of the
More informationEXPANSIONS FOR QUANTILES AND MULTIVARIATE MOMENTS OF EXTREMES FOR HEAVY TAILED DISTRIBUTIONS
REVSTAT Statistical Journal Volume 15, Number 1, January 2017, 25 43 EXPANSIONS FOR QUANTILES AND MULTIVARIATE MOMENTS OF EXTREMES FOR HEAVY TAILED DISTRIBUTIONS Authors: Christopher Withers Industrial
More informationA LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS
Applicable Analysis and Discrete Mathematics, 1 (27, 371 385. Available electronically at http://pefmath.etf.bg.ac.yu A LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS H. W. Gould, Jocelyn
More informationPreservation of log-concavity on summation
Preservation of log-concavity on summation Oliver Johnson Christina Goldschmidt October 12, 2005 Abstract We extend Hoggar s theorem that the sum of two independent discrete-valued log-concave random variables
More informationMULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng
Indian J. pure appl. Math., 39(: 37-55, April 008 c Printed in India. MULTISECTION METHOD AND FURTHER FORMULAE FOR π De-Yin Zheng Department of Mathematics, Hangzhou Normal University, Hangzhou 30036,
More informationA196837: Ordinary Generating Functions for Sums of Powers of the First n Positive Integers
Karlsruhe October 14, 2011 November 1, 2011 A196837: Ordinary Generating Functions for Sums of Powers of the First n Positive Integers Wolfdieter L a n g 1 The sum of the k th power of the first n positive
More informationStirling Numbers of the 1st Kind
Daniel Reiss, Colebrook Jackson, Brad Dallas Western Washington University November 28, 2012 Introduction The set of all permutations of a set N is denoted S(N), while the set of all permutations of {1,
More information18.175: Lecture 17 Poisson random variables
18.175: Lecture 17 Poisson random variables Scott Sheffield MIT 1 Outline More on random walks and local CLT Poisson random variable convergence Extend CLT idea to stable random variables 2 Outline More
More informationCalculation of Steady-State Probabilities of M/M Queues: Further Approaches. Queuing Theory and Applications (MA597)
Calculation of Steady-State Probabilities of M/M Queues: Further Approaches A assignment report submitted for the course Queuing Theory and Applications (MA597) by Ambati Narendar Reddy Roll No. 06212301
More informationOn an asymptotic series of Ramanujan
Ramanujan J (2009) 20: 179 188 DOI 10.1007/s11139-009-9169-x On an asymptotic series of Ramanujan Yaming Yu Received: 22 May 2008 / Accepted: 1 April 2009 / Published online: 12 August 2009 The Author(s)
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables THE UNIVERSITY OF MANCHESTER. 21 June :45 11:45
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS 21 June 2010 9:45 11:45 Answer any FOUR of the questions. University-approved
More informationASYMPTOTICS OF THE EULERIAN NUMBERS REVISITED: A LARGE DEVIATION ANALYSIS
ASYMPTOTICS OF THE EULERIAN NUMBERS REVISITED: A LARGE DEVIATION ANALYSIS GUY LOUCHARD ABSTRACT Using the Saddle point method and multiseries expansions we obtain from the generating function of the Eulerian
More informationSeries of Error Terms for Rational Approximations of Irrational Numbers
2 3 47 6 23 Journal of Integer Sequences, Vol. 4 20, Article..4 Series of Error Terms for Rational Approximations of Irrational Numbers Carsten Elsner Fachhochschule für die Wirtschaft Hannover Freundallee
More informationUniversity of Lisbon, Portugal
Development and comparative study of two near-exact approximations to the distribution of the product of an odd number of independent Beta random variables Luís M. Grilo a,, Carlos A. Coelho b, a Dep.
More informationMathematical Statistics 1 Math A 6330
Mathematical Statistics 1 Math A 6330 Chapter 3 Common Families of Distributions Mohamed I. Riffi Department of Mathematics Islamic University of Gaza September 28, 2015 Outline 1 Subjects of Lecture 04
More informationThe Design of a Survival Study
The Design of a Survival Study The design of survival studies are usually based on the logrank test, and sometimes assumes the exponential distribution. As in standard designs, the power depends on The
More informationCONVOLUTION TREES AND PASCAL-T TRIANGLES. JOHN C. TURNER University of Waikato, Hamilton, New Zealand (Submitted December 1986) 1.
JOHN C. TURNER University of Waikato, Hamilton, New Zealand (Submitted December 986). INTRODUCTION Pascal (6-66) made extensive use of the famous arithmetical triangle which now bears his name. He wrote
More informationExplicit formulas for computing Bernoulli numbers of the second kind and Stirling numbers of the first kind
Filomat 28:2 (24), 39 327 DOI.2298/FIL4239O Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Explicit formulas for computing Bernoulli
More informationCombinatorial Proof of the Hot Spot Theorem
Combinatorial Proof of the Hot Spot Theorem Ernie Croot May 30, 2006 1 Introduction A problem which has perplexed mathematicians for a long time, is to decide whether the digits of π are random-looking,
More informationEXPLICIT INVERSE OF THE PASCAL MATRIX PLUS ONE
EXPLICIT INVERSE OF THE PASCAL MATRIX PLUS ONE SHENG-LIANG YANG AND ZHONG-KUI LIU Received 5 June 005; Revised September 005; Accepted 5 December 005 This paper presents a simple approach to invert the
More informationIf the [T, Id] automorphism is Bernoulli then the [T, Id] endomorphism is standard
If the [T, Id] automorphism is Bernoulli then the [T, Id] endomorphism is standard Christopher Hoffman Daniel Rudolph September 27, 2002 Abstract For any - measure preserving map T of a probability space
More informationLengyel's Constant Page 1 of 5
Lengyel's Constant Page 1 of 5 Lengyel's Constant Set Partitions and Bell Numbers Let S be a set with n elements. The set of all subsets of S has elements. By a partition of S we mean a disjoint set of
More informationP i [B k ] = lim. n=1 p(n) ii <. n=1. V i :=
2.7. Recurrence and transience Consider a Markov chain {X n : n N 0 } on state space E with transition matrix P. Definition 2.7.1. A state i E is called recurrent if P i [X n = i for infinitely many n]
More informationSOME RECURRENCE RELATIONS BETWEEN PRODUCT MOMENTS OF ORDER STATISTICS OF A DOUBLY TRUNCATED PARETO DISTRIBUTION
Sankhyā : The Indian Journal of Statistics 1995, Volume 57, Series B, Pt. 1, pp. 1-9 SOME RECURRENCE RELATIONS BETWEEN PRODUCT MOMENTS OF ORDER STATISTICS OF A DOUBLY TRUNCATED PARETO DISTRIBUTION By A.
More informationBias-corrected AIC for selecting variables in Poisson regression models
Bias-corrected AIC for selecting variables in Poisson regression models Ken-ichi Kamo (a), Hirokazu Yanagihara (b) and Kenichi Satoh (c) (a) Corresponding author: Department of Liberal Arts and Sciences,
More informationWeek 9 The Central Limit Theorem and Estimation Concepts
Week 9 and Estimation Concepts Week 9 and Estimation Concepts Week 9 Objectives 1 The Law of Large Numbers and the concept of consistency of averages are introduced. The condition of existence of the population
More informationInvestigating Geometric and Exponential Polynomials with Euler-Seidel Matrices
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 14 (2011), Article 11.4.6 Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices Ayhan Dil and Veli Kurt Department of Mathematics
More informationOn the stirling expansion into negative powers of a triangular number
MATHEMATICAL COMMUNICATIONS 359 Math. Commun., Vol. 5, No. 2, pp. 359-364 200) On the stirling expansion into negative powers of a triangular number Cristinel Mortici, Department of Mathematics, Valahia
More information#A29 INTEGERS 11 (2011) THE (EXPONENTIAL) BIPARTITIONAL POLYNOMIALS AND POLYNOMIAL SEQUENCES OF TRINOMIAL TYPE: PART II
#A29 INTEGERS 11 (2011) THE (EXPONENTIAL) BIPARTITIONAL POLYNOMIALS AND POLYNOMIAL SEQUENCES OF TRINOMIAL TYPE: PART II Hacène Belbachir 1 USTHB, Faculty of Mathematics, El Alia, Bab Ezzouar, Algeria hbelbachir@usthb.dz
More informationThe Number of Inversions in Permutations: A Saddle Point Approach
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 6 (23), Article 3.2.8 The Number of Inversions in Permutations: A Saddle Point Approach Guy Louchard Département d Informatique CP 212, Boulevard du
More informationMonotonicity and Aging Properties of Random Sums
Monotonicity and Aging Properties of Random Sums Jun Cai and Gordon E. Willmot Department of Statistics and Actuarial Science University of Waterloo Waterloo, Ontario Canada N2L 3G1 E-mail: jcai@uwaterloo.ca,
More informationCharacterizations of Weibull Geometric Distribution
Journal of Statistical Theory and Applications Volume 10, Number 4, 2011, pp. 581-590 ISSN 1538-7887 Characterizations of Weibull Geometric Distribution G. G. Hamedani Department of Mathematics, Statistics
More informationSome Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution
Journal of Probability and Statistical Science 14(), 11-4, Aug 016 Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution Teerawat Simmachan
More informationOn a simple construction of bivariate probability functions with fixed marginals 1
On a simple construction of bivariate probability functions with fixed marginals 1 Djilali AIT AOUDIA a, Éric MARCHANDb,2 a Université du Québec à Montréal, Département de mathématiques, 201, Ave Président-Kennedy
More informationInference on reliability in two-parameter exponential stress strength model
Metrika DOI 10.1007/s00184-006-0074-7 Inference on reliability in two-parameter exponential stress strength model K. Krishnamoorthy Shubhabrata Mukherjee Huizhen Guo Received: 19 January 2005 Springer-Verlag
More informationClinical Trial Design Using A Stopped Negative Binomial Distribution. Michelle DeVeaux and Michael J. Kane and Daniel Zelterman
Statistics and Its Interface Volume 0 (0000) 1 10 Clinical Trial Design Using A Stopped Negative Binomial Distribution Michelle DeVeaux and Michael J. Kane and Daniel Zelterman arxiv:1508.01264v8 [math.st]
More informationORDER RESTRICTED STATISTICAL INFERENCE ON LORENZ CURVES OF PARETO DISTRIBUTIONS. Myongsik Oh. 1. Introduction
J. Appl. Math & Computing Vol. 13(2003), No. 1-2, pp. 457-470 ORDER RESTRICTED STATISTICAL INFERENCE ON LORENZ CURVES OF PARETO DISTRIBUTIONS Myongsik Oh Abstract. The comparison of two or more Lorenz
More informationAvoider-Enforcer games played on edge disjoint hypergraphs
Avoider-Enforcer games played on edge disjoint hypergraphs Asaf Ferber Michael Krivelevich Alon Naor July 8, 2013 Abstract We analyze Avoider-Enforcer games played on edge disjoint hypergraphs, providing
More informationHIGHER-ORDER DIFFERENCES AND HIGHER-ORDER PARTIAL SUMS OF EULER S PARTITION FUNCTION
ISSN 2066-6594 Ann Acad Rom Sci Ser Math Appl Vol 10, No 1/2018 HIGHER-ORDER DIFFERENCES AND HIGHER-ORDER PARTIAL SUMS OF EULER S PARTITION FUNCTION Mircea Merca Dedicated to Professor Mihail Megan on
More informationChapter 6 Expectation and Conditional Expectation. Lectures Definition 6.1. Two random variables defined on a probability space are said to be
Chapter 6 Expectation and Conditional Expectation Lectures 24-30 In this chapter, we introduce expected value or the mean of a random variable. First we define expectation for discrete random variables
More information2 The Polygonal Distribution
2 The Polygonal Distribution Dimitris Karlis 1 and Evdokia Xekalaki 2 1 Department of Statistics, Athens University of Economics and Business 2 Department of Statistics, Athens University of Economics
More informationON THE SUM OF POWERS OF TWO. 1. Introduction
t m Mathematical Publications DOI: 0.55/tmmp-06-008 Tatra Mt. Math. Publ. 67 (06, 4 46 ON THE SUM OF POWERS OF TWO k-fibonacci NUMBERS WHICH BELONGS TO THE SEQUENCE OF k-lucas NUMBERS Pavel Trojovský ABSTRACT.
More informationDiscrete Distributions
Discrete Distributions STA 281 Fall 2011 1 Introduction Previously we defined a random variable to be an experiment with numerical outcomes. Often different random variables are related in that they have
More informationPage Max. Possible Points Total 100
Math 3215 Exam 2 Summer 2014 Instructor: Sal Barone Name: GT username: 1. No books or notes are allowed. 2. You may use ONLY NON-GRAPHING and NON-PROGRAMABLE scientific calculators. All other electronic
More informationThe closed-form resolution of integrals is one of the
Construction of Generalized Integral Formulas By Means of Laplace Transformations Adam C. Buss, Department of Physics, Indiana University - Northwest Faculty Mentor: Dr. Axel Schulze-Halberg ABSTRACT We
More informationSPECTRAL GAP FOR ZERO-RANGE DYNAMICS. By C. Landim, S. Sethuraman and S. Varadhan 1 IMPA and CNRS, Courant Institute and Courant Institute
The Annals of Probability 996, Vol. 24, No. 4, 87 902 SPECTRAL GAP FOR ZERO-RANGE DYNAMICS By C. Landim, S. Sethuraman and S. Varadhan IMPA and CNRS, Courant Institute and Courant Institute We give a lower
More informationSection 4.2: Mathematical Induction 1
Section 4.: Mathematical Induction 1 Over the next couple of sections, we shall consider a method of proof called mathematical induction. Induction is fairly complicated, but a very useful proof technique,
More informationA process capability index for discrete processes
Journal of Statistical Computation and Simulation Vol. 75, No. 3, March 2005, 175 187 A process capability index for discrete processes MICHAEL PERAKIS and EVDOKIA XEKALAKI* Department of Statistics, Athens
More informationx arctan x = x x x x2 9 +
Math 1B Project 3 Continued Fractions. Many of the special functions that occur in the applications of mathematics are defined by infinite processes, such as series, integrals, and iterations. The continued
More informationIntroduction to Combinatorial Mathematics
Introduction to Combinatorial Mathematics George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 300 George Voutsadakis (LSSU) Combinatorics April 2016 1 / 57
More informationTHE EXPONENTIAL DISTRIBUTION ANALOG OF THE GRUBBS WEAVER METHOD
THE EXPONENTIAL DISTRIBUTION ANALOG OF THE GRUBBS WEAVER METHOD ANDREW V. SILLS AND CHARLES W. CHAMP Abstract. In Grubbs and Weaver (1947), the authors suggest a minimumvariance unbiased estimator for
More information18.175: Lecture 13 Infinite divisibility and Lévy processes
18.175 Lecture 13 18.175: Lecture 13 Infinite divisibility and Lévy processes Scott Sheffield MIT Outline Poisson random variable convergence Extend CLT idea to stable random variables Infinite divisibility
More informationExistence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin method
Int. J. Nonlinear Anal. Appl. 7 (6) No., 7-8 ISSN: 8-68 (electronic) http://dx.doi.org/.75/ijnaa.5.37 Existence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin
More informationEXPLICIT EXPRESSIONS FOR MOMENTS OF χ 2 ORDER STATISTICS
Bulletin of the Institute of Mathematics Academia Sinica (New Series) Vol. 3 (28), No. 3, pp. 433-444 EXPLICIT EXPRESSIONS FOR MOMENTS OF χ 2 ORDER STATISTICS BY SARALEES NADARAJAH Abstract Explicit closed
More informationResearch Article The Laplace Likelihood Ratio Test for Heteroscedasticity
International Mathematics and Mathematical Sciences Volume 2011, Article ID 249564, 7 pages doi:10.1155/2011/249564 Research Article The Laplace Likelihood Ratio Test for Heteroscedasticity J. Martin van
More informationFrontier estimation based on extreme risk measures
Frontier estimation based on extreme risk measures by Jonathan EL METHNI in collaboration with Ste phane GIRARD & Laurent GARDES CMStatistics 2016 University of Seville December 2016 1 Risk measures 2
More information