Next Article in Journal
Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam-Type Stability
Previous Article in Journal
Has the COVID-19 Pandemic Affected Mathematics Achievement? A Case Study of University Students in Social Sciences
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiple q-Integral and Records from Geometrically Distributed Sequences

by
Fernando López-Blázquez
*,†,‡ and
Begoña Salamanca-Miño
Departamento de Estadística e I.O., Universidad de Sevilla, 41012 Seville, Spain
*
Author to whom correspondence should be addressed.
Current address: Facultad de Matemáticas, Avda. Reina Mercedes s/n, 41012 Sevilla, Spain.
These authors contributed equally to this work.
Mathematics 2022, 10(13), 2313; https://doi.org/10.3390/math10132313
Submission received: 1 April 2022 / Revised: 29 June 2022 / Accepted: 30 June 2022 / Published: 1 July 2022
(This article belongs to the Section Probability and Statistics)

Abstract

:
We study the distribution theory of some statistics related to records from a sequence of geometrically distributed random variables. The main novelty in this study is the introduction of q-calculus techniques. We obtain representations for the probability functions of some variables of interest, such as, for instance, record times, number of records, inter-record times, record indicators, etc., in terms of multiple q-integrals. It is remarkable to note that the expressions thus obtained are q-analogues of the corresponding well-known results in the absolutely continuous case. We also explore a duality between the results in the case of weak and ordinary records.

1. Introduction

The origins of q-calculus date back to the works of [1] in connection to infinite series in his famous book Introductio in Analysin Infinitorum. In the last few decades, there has been a considerable interest in the topic due to its numerous applications in different areas of mathematics and physics, such as number theory, special functions, orthogonal polynomials, operator theory, combinatorics, probability, quantum mechanics, relativity, string theory, etc. An extensive treatment on q-calculus, its history and applications can be found in [2].
In this paper, we explore a connection between q-calculus and the distribution theory of records from sequences of random variables geometrically distributed.
Ref. [3] started the probabilistic study of record values and gave many of their basic distributional properties. Let us recall some definitions and results. Let { X n } n 1 be a sequence of independent and identically distributed (iid) random variables with the common distribution function F. The sequence of (upper) ordinary record times, { T m } m 0 , is defined as T m = min { j > T m 1 : X j > X T m 1 } , m 1 with T 0 = 1 . The m-th (ordinary) record of the sequence is the random variable R m = X T m , m 0 . For n 2 , the record indicators are the binary random variables
I n = 1 , if   X n > max { X 1 , , X n 1 } 0 , otherwise
and I 1 = 1 with probability one. In the case in which F is continuous, the basic distribution theory of the statistics above defined is very well-known; see, for instance, chapter 2 of [4]. In the continuous case, the distributions of the record indicators, the record times and the inter-record times do not depend on the parent distribution, and we have results like the following (for further references, see [5]):
  • The record indicators are independent Bernoulli random variables with P ( I n = 1 ) = 1 / n , n 1 . Then
    P I 1 = 1 , I 2 = 1 , , I n = 1 = 1 n ! .
  • The joint probability mass function (pmf) of the first m (non-trivial) record times is
    P ( T 1 = t 1 , T 2 = t 2 , , T m = t m ) = 1 ( t 1 1 ) ( t 2 1 ) ( t m 1 ) t m ,
    for integers 1 < t 1 < t 2 < < t m ; also
    P ( T m = n ) = | s n 1 , m | n ! ,
    with s n , k being the Stirling numbers of the first kind defined by the equality
    x ( x 1 ) ( x n + 1 ) = k 0 s n , k x k .
  • The pmf of the inter-record times, S m = T m T m 1 , m 1 satisfies (see [6])
    P ( S m > s ) = 1 ( m 1 ) ! 0 x m 1 e x ( 1 e x ) s d x , s 0 ,
    so then, for s 1 ,
    P ( S m = s ) = P ( S m > s 1 ) P ( S m > s ) = 1 ( m 1 ) ! 0 x m 1 e 2 x ( 1 e x ) s 1 d x = 1 ( m 1 ) ! = 0 s 1 s 1 ( 1 ) 0 x m 1 e ( + 2 ) x d x = = 0 s 1 s 1 ( 1 ) ( + 2 ) m .
  • The number of records in a sequence of length n, denoted N n , has pmf
    P ( N n = m ) = | s n , m | n ! , m = 1 , , n ,
    and its expected value is E N n = j = 1 n 1 j .
If the parent distribution, F, is discrete, the previous results are no longer true, and in general, the distribution of record times, inter-record times, record indicators and number of records depend on F. Thus, the distribution theory for records from discrete parents is more complicated than in the continuous case. In fact, there exist in the literature two definitions of records, ordinary and weak records, which are equivalent in the case of continuous parents, but not so in the discrete case. The weak record time sequence is defined as T ˜ m = min { j > T ˜ m 1 : X j X T ˜ m 1 } , m 1 , and T ˜ 0 = 1 . Then, the m t h weak record of the sequence { X n } n 1 (or from the parent distribution F) is the random variable R ˜ m = X T ˜ m , m 0 . Note that the difference with ordinary records is that ties with the previous records are also considered as records. Obviously, for continuous distributions, weak records coincide (almost surely) with ordinary records. In the case of discrete distributions, ordinary and weak records lead to different definitions.
We say that a random variable, X, follows a geometric distribution with parameter p ( 0 , 1 ) if its pmf is P ( X = x ) = q x p for x = 0 , 1 , and q = 1 p . We denote X Ge ( p ) .
The number of records in geometrically distributed sequences has been studied in connection with a data structure called skiplist; see [7,8] or [9]. Related studies can be found in [10], which established the asymptotic normality of the number of records from geometric parents and also studied large deviations, local limit theorems and approximations. Some of these results have been extended to other discrete parents, so [11] established a CLT for the number of records from discrete distributions; rates of growth for the number of records and laws of large numbers can be found in [12,13].
Most of the results above deal with limit properties, and not very much attention has been paid to the exact distribution of statistics related to the number of records. Our main contribution in this paper is the use of techniques of q-calculus for the study of the distribution theory of records from geometric parents.
In Section 2, we present the Jackson integral or q-integral and the extension of this concept to multiple q-integration in a convex polytope; see [14]. We also present in this section some useful lemmas for the calculation of multiple q-integrals (proofs are deferred to the Appendix A) and the non-central q-Stirling numbers; see [15]. No previous knowledge of q-calculus is necessary to understand this section. In Section 3, we study some statistics related to the sequence of weak records from a sequence of geometrically distributed random variables. Theorem 1 is a central result in this section. It gives a multiple q-integral representation of the joint probability function of the (weak) record times. As a consequence of this result, we show that the distribution of record times and that of the number of records in a sample of size n can be written in terms of the central and non-central q-Stirling numbers. We also study the distribution of inter-(weak) record times and we prove its asymptotic log-normal behavior. In Section 4, it becomes clear that the results obtained for a geometric parent, Ge ( p ) , are q-analogues of the corresponding results in the case of an absolutely continuous distribution, so we can say that the absolutely continuous case is the limit case of geometric parents when q = 1 p 1 . Finally, in Section 5, we consider the case of ordinary records. We show that there is the following duality: the formulas in the case of ordinary records can be obtained from the corresponding ones in the case of weak records by replacing q = 1 p with Q : = 1 / q .
Throughout this paper, we will use the usual conventions with sums and products: for a > b , a b · = 0 and a b · = 1 , and more generally, the summation (product) over an empty set of indexes is equal to zero (one). Additionally, 0 0 = 1 .

2. The Jackson Integral and q -Analogues

In the following, we present some basic concepts of q-calculus. The interested reader can find more details in [16], or, for a deeper treatment of the subject, see [2]. Regardless, this section contains all the results of q-calculus that we need in this paper.
The most basic concept in q-calculus is perhaps that of the q-number. Let q R { 0 } . The q-numbers are defined as
x q : = ( 1 q x ) / ( 1 q ) , q 1 x , q = 1
for x R . For a positive integer, m, its q-factorial is m q ! = j = 1 m j q and 0 q ! = 1 . Throughout this paper, Q : = 1 / q . Thus, for instance, for x R , we have x Q = x q / q x 1 .
Given a function f defined on the non-negative semi-axis of the real line, a , b > 0 and q ( 0 , 1 ) , the Jackson integral or q-integral of f in the interval ( a , b ) is defined as
a b f ( y ) d q y : = ( 1 q ) = 0 f ( q b ) b q f ( q a ) a q
provided that the series in the right hand side (rhs) of (5) is absolutely convergent.
The Jackson integral has some resemblances to the Riemann integral; for instance, for α > 1 ,
0 x y α d q y = x α + 1 α + 1 q .
Among the differences with respect to the ordinary Riemann integrals, we mention that the formula of the change of variables is no longer valid, so, for natural α ,
1 α + 1 q 0 1 ( 1 y ) α d q y = j = 0 α α j ( 1 ) j 1 j + 1 q .
The definition of a q-integral can be extended to a multiple q-integral in a convex polytope; see [14]. Let us denote as S m the set of permutations of { 1 , , m } . Let π = ( π 1 π m ) S m . For a function of m variables, f ( y 1 , , y m ) and a convex polytope C R m , the q-integral of f over C with respect to the order π of integration is defined as
C f ( y 1 , , y m ) d q y π 1 d q y π m = min ( C m ) max ( C m ) min ( C 1 ) max ( C 1 ) f ( y 1 , , y m ) d q y π 1 d q y π m
where C i , i = 1 , , m 1 is the set of real numbers depending on the values of y π i + 1 , , y π m with
C i : = C i ( y π i + 1 , , y π m ) = { x π i : ( x 1 , , x m ) C , x π j = y π j , for j > i } ,
and C m : = { x π m : ( x 1 , , x m ) C } (implicitly, hereafter we assume the absolute convergence of the involved series).
It follows from the definition that when the polytope C is a hyper-rectangle, i.e., C = [ a 1 , b 1 ] × × [ a m , b m ] , the ordering π is irrelevant; that is to say, we can change the order of integration, and so, for any π S m ,
[ 0 , 1 ] m f ( y 1 , , y m ) d q y π 1 d q y π m : = ( 1 q ) m 1 , , m 0 f ( q 1 , , q m ) q i = 1 m i ,
for f, a function in m variables, such that the series in the rhs of (8) is absolutely convergent. Thus, for example,
[ 0 , 1 ] m j = 1 m y j α j d q y π 1 d q y π m = [ 0 , 1 ] m j = 1 m y j α j d q y 1 d q y m = j = 1 m 1 α j + 1 q
for any π S m and α j > 1 , j = 1 , , m .
The next Lemmas are useful in connection with iterated and multiple Jackson integrals.
Lemma 1.
Let f be a real function of m variables such that the series appearing below are absolutely convergent. Then, for any π S m ,
(a) 
( 1 q ) m 0 x 1 x m < f ( q x 1 , q x 2 , , q x m ) q i = 1 m x i
= [ 0 , 1 ] m f ( y 1 , y 1 y 2 , , y 1 y m ) j = 1 m y j m j d q y π 1 d q y π m .
(b) 
( 1 q ) m 0 x 1 < < x m < f ( q x 1 , q x 2 , , q x m ) q i = 1 m x i
= q m 2 [ 0 , 1 ] m f ( y 1 , q y 1 y 2 , , q m 1 y 1 y m ) j = 1 m y j m j d q y π 1 d q y π m .
In general, when the polytope of integration is not a hyper-rectangle, the order of integration is important. We illustrate this fact with the following example:
0 y 2 y 1 1 d q y 2 d q y 1 = 0 1 0 y 1 d q y 2 d q y 1 = 0 1 y 1 d q y 1 = 1 2 q ,
but,
0 y 2 y 1 1 d q y 1 d q y 2 = 0 1 y 2 1 d q y 1 d q y 2 = 0 1 ( 1 y 2 ) d q y 2 = 1 1 2 q .
Throughout this paper we will use the reverse ordering, d q y ( m ) : = d q y m d q y m 1 d q y 1 . For q ( 0 , 1 ] , we define the polytopes
D m , q : = { ( y 1 , , y m ) R n : 0 y 1 1 ; 0 y j q y j 1 , j = 2 , , m } , D m : = D m , 1 = { ( y 1 , , y m ) R n : 0 y m y m 1 y 1 1 } ,
and from the definition of a multiple q-integral given in (7),
D m f ( y 1 , , y m ) d q y ( m ) = 0 1 0 y 1 0 y m 1 f ( y 1 , , y m ) d q y ( m ) ,
D m , q f ( y 1 , , y m ) d q y ( m ) = 0 1 0 q y 1 0 q y m 1 f ( y 1 , , y m ) d q y ( m ) .
Lemma 2.
Let f be a real function of m variables such that the series appearing below are absolutely convergent. Then,
(a) 
( 1 q ) m 0 x 1 x m < f ( q x 1 , , q x m ) q j = 1 m x j = D m f ( y 1 , , y m ) d q y ( m ) .
(b) 
( 1 q ) m 0 x 1 < < x m < f ( q x 1 , , q x m ) q j = 1 m x j = D m , q f ( y 1 , , y m ) d q y ( m ) .
Comparing the results in Lemmas 1 and 2, we obtain
D m f ( y 1 , , y m ) d q y ( m ) = [ 0 , 1 ] m f ( y 1 , y 1 y 2 , , y 1 y m ) j = 1 m y j m j d q y ( m )
D m , q f ( y 1 , , y m ) d q y ( m ) = q m 2 [ 0 , 1 ] m f ( y 1 , q y 1 y 2 , , q m 1 y 1 y m ) j = 1 m y j m j d q y ( m ) ,
and using these results in combination with (6), it can be confirmed that for m 1 and β 1 , , β m , real numbers such that β j + + β m + m j > 1 , for j = 1 , , m ,
D m j = 1 m y j β j d q y ( m ) = j = 1 m 1 β j + + β m + m j + 1 q ,
D m , q j = 1 m y j β j d q y ( m ) = q m 2 j = 1 m q ( j 1 ) β j β j + + β m + m j + 1 q = Q β 1 + + β m j = 1 m 1 β j + + β m + m j + 1 Q ,
and in particular, for β 1 = = β m = 0 ,
D m d q y ( m ) = 1 m q ! ,
D m , q d q y ( m ) = 1 m Q ! ,
and for β 1 = = β m 1 = 1 and β m = b > 1 ,
D m y m b j = 1 m 1 y j d q y ( m ) = 1 b + 1 q m ,
D m , q y m b j = 1 m 1 y j d q y ( m ) = Q b m + 1 b + 1 Q m .
Given a (formal) power series A ( z ) = k 0 a k z k , we denote [ z k ] A ( z ) = a k , k 0 . For non-negative integers, n and r, the non-central q-Stirling numbers of the first kind (see [15]) are defined as
s q ( n , k ; r ) = [ z k ] = r n + r 1 ( z q ) , k = 0 , , n ,
and by convention s q ( n , k ; r ) = 0 if k < 0 or k > n or n < 0 . The case r = 0 corresponds to the q-Stirling numbers of the first kind; see [2], Ch. 5. We denote s q ( n , k ) : = s q ( n , k ; 0 ) . The non-central q-Stirling numbers satisfy the recurrence
s q ( n + 1 , k ; r ) = s q ( n , k 1 ; r ) r + n q s q ( n , k ; r ) , n 0 ; k = 0 , , n + 1 , s q ( 0 , 0 ; r ) = 1 .
Some useful relationships for our purposes are the following:
| s q ( n , k ; r ) | = ( 1 ) n k s q ( n , k ; r ) ,
s q ( k , k ; r ) = 1 , k 0 ,
s q ( n , k ; 1 ) = s q ( n 1 , k ; 2 ) s q ( n 1 , k ; 2 ) ,
s q ( n , k ) = s q ( n 1 , k 1 ; 1 ) .
By equating coefficients of z in both sides of
1 + z r q 1 + z n + r 1 q = r 1 q ! n + r 1 q ! ( z + r q ) ( z + n + r 1 q )
we obtain
r c 1 < < c n + r 1 1 c j q = ( 1 ) n r 1 q ! s q ( n , ; r ) n + r 1 q ! .
Lemma 3.
For the integers m 1 , k 0 and q ( 0 , 1 ) ,
(a) 
α 1 , , α m 0 α 1 + α m = k D m j = 1 m y j α j d q y ( m ) = ( 1 ) k s q ( m + k 1 , m 1 ; 2 ) m + k + 1 q ! .
(b) 
α 1 , , α m 0 α 1 + α m = k D m , q j = 1 m y j α j d q y ( m ) = ( 1 ) k s Q ( m + k 1 , m 1 ; 2 ) m + k + 1 Q ! , where Q : = 1 / q .
Lemma 4.
For m > 1 and f, a non-negative function defined in the non-negative real numbers,
(a) 
s 1 , , s m 1 1 D m f ( y m ) j = 1 m 1 ( 1 y j ) s j 1 d q y ( m ) = D m f ( y m ) y 1 y m 1 d q y ( m ) .
(b) 
s 1 , , s m 1 1 D m , q f ( y m ) j = 1 m 1 ( 1 q y j ) s j 1 d q y ( m ) = 1 q m 1 D m , q f ( y m ) y 1 y m 1 d q y ( m ) .
Many variants of Lemma 4 (all of them with similar proofs) are possible; for instance, for a non-negative function f i , i = 1 , , m ,
s 1 D m f 1 ( y 1 ) f m ( y m ) ( 1 y m ) s 1 d q y ( m ) = D m f 1 ( y 1 ) f m ( y m ) y m d q y ( m ) .
Note that results like (25) or those in Lemma 4 justify the formal interchange of the summation and integration symbols.
Suppose that A ( q ) is an expression involving the parameter q. Then, it is said that A ( q ) is a q-analogue of A if A ( q ) A when q 1 . q-analogues appear in different areas of mathematics and physics, such as, for instance, in quantum groups, string theory, fractals, dynamical systems, elliptic curves, etc. It is clear that q-numbers are q-analogues of the corresponding real numbers. Expression (6) is a q-analogue of the well-known formula for Riemann integrals 0 x y α d y = x α + 1 / ( α + 1 ) , α > 1 . Similarly, (9), (14)–(17) are q-analogues of the corresponding expressions obtained by substituting q-integrals for Riemann integrals and q-numbers for real numbers. Then, the following lemma is immediate.
Lemma 5.
If R ( y 1 , , y m ) is a polynomial with real coefficients in the variables y 1 , , y m , then the multiple q-integral
[ 0 , 1 ] m R ( y 1 , , y m ) d q y ( m ) or D m R ( y 1 , , y m ) d q y ( m )
is a q-analogue of the multiple Riemann integral
[ 0 , 1 ] m R ( y 1 , , y m ) d y m d y m 1 d y 1 or D m R ( y 1 , , y m ) d y m d y m 1 d y 1 .
Finally, we show a connection between multiple q-integrals and expected values involving geometric samples. Let X i Ge ( p ) , i 1 be independent random variables with geometric distributions of the parameter p ( 0 , 1 ) (i.e., P ( X i = k ) = ( 1 p ) k p , k = 0 , 1 , ). Let q = 1 p . Let f : R m R be a function such that E f ( q X 1 , , q X m ) < + , and thus, from (8) we get
E f ( q X 1 , , q X m ) = [ 0 , 1 ] m f ( y 1 , , y m ) d q y ( m ) .

3. Weak Records from Geometric Parent

Let X = { X n } n 1 be a sequence of iid random variables with a common distribution Ge ( p ) with p ( 0 , 1 ) and q = 1 p . The index n in the process X = { X n } n 1 can be regarded as discrete time. A trajectory of the process X is a sequence of non-negative integers. A trajectory up to time n is a finite sequence of length n of non-negative integers. We will denote a finite trajectory of length n 1 by a sequence ( k 1 , , k n ) where k i , i 1 are non-negative integers. Consecutive equal integers in a trajectory are denoted with exponents; for instance, ( k ) t is the sequence ( k , , k ) of length t. A sequence of length 1 formed with a non-negative integer less than k is denoted by ( < k ) . Similar notations are used in the cases ≤, >, ≥, etc. As an example of these notations, consider the trajectory ( 3 , 1 , 1 , 5 , 4 , 2 , 3 , 1 , 8 ) . Whenever it is convenient, this trajectory could be denoted as: ( 3 ) ( 1 ) 2 ( 5 ) 5 ( 8 ) or ( 3 ) 3 ( 3 ) 2 ( 8 ) 4 , or ( 3 ) ( 1 ) 8 ,…etc.
Theorem 1.
For m 1 ,
(a) 
For 1 = t 0 < t 1 < < t m , the joint pmf of ( T ˜ 1 , , T ˜ m ) from X admits the following q-integral representation:
P T ˜ 1 = t 1 , , T ˜ m = t m = D m y m j = 1 m ( 1 y j ) t j t j 1 1 d q y ( m ) .
(b) 
The joint probability generating function of ( T ˜ 1 , , T ˜ m ) is
G ( z 1 , , z m ) = D m y m j = 1 m z j j + 1 1 ( 1 y j ) z j z m d q y ( m ) ,
for | z j | 1 , j = 1 , , m .
Proof. 
(a) For the integers 1 = t 0 < t 1 < < t m , m 1 , the event { T ˜ 1 = t 1 , , T ˜ m = t m } corresponds to finite trajectories of length t m m + 1 of the process X , which are of the form:
( x 1 ) ( < x 1 ) t 1 2 ( x 2 ) ( < x 2 ) t 2 t 1 1 ( x m ) ( < x m ) t m t m 1 1 ( x m + 1 ) ,
with 0 x 1 x 2 x m + 1 . The probability associated with the trajectory given in (28) is j = 1 m q x j p ( 1 q x j ) t j t j 1 1 q x m + 1 p . Then, using Lemma 2(a),
P T ˜ 1 = t 1 , , T ˜ m = t m = 0 x 1 x 2 x m + 1 < j = 1 m q x j p ( 1 q x j ) t j t j 1 1 q x m + 1 p = D m + 1 j = 1 m ( 1 y j ) t j t j 1 1 d q y ( m + 1 ) = D m j = 1 m ( 1 y j ) t j t j 1 1 0 y m d q y m + 1 d q y m = D m y m j = 1 m ( 1 y j ) t j t j 1 1 d q y ( m ) .
(b) Firstly, observe that for the integers s 1 , , s m , and m 1 ,
j = 1 m z j s 1 + + s j + 1 = z 1 2 z 2 3 z m m + 1 j = 1 m l = j m z s j 1
For | z j | 1 and s j = t j t j 1 , j = 1 , , m , ( t 0 = 1 ) . Using Theorem 1(a) and (29), we have
G ( z 1 , , z m ) = 1 = t 0 < t 1 < < t m < + P T ˜ 1 = t 1 , , T ˜ m = t m j = 1 m z j t j = s 1 , , s m 0 P T ˜ 1 = s 1 + 1 , T ˜ 2 = s 1 + s 2 + 1 , , T ˜ m = s 1 + + s m + 1 j = 1 m z j s 1 + + s j + 1 = z 1 2 z 2 3 z m m + 1 s 1 , , s m 0 D m y m j = 1 m ( 1 y j ) = j m z s j 1 d q y ( m ) = j = 1 m z j j + 1 D m y m s 1 , , s m 0 j = 1 m ( 1 y j ) z j z m s j 1 d q y ( m ) = D m y m j = 1 m z j j + 1 1 ( 1 y j ) z j z m d q y ( m ) ,
where the interchange of the summation and integral symbols can be justified with a similar argument to the one in the proof of Lemma 4.  □
Theorem 2.
For m 1 and n m + 1 ,
(a) 
P ( T ˜ m = n ) = = m 1 n 2 n 2 s q ( , m 1 ; 2 ) + 2 q !
(b) 
P ( T ˜ m n ) = = m 1 n 2 n 1 + 1 s q ( , m 1 ; 2 ) + 2 q !
Proof. 
(a)
From Theorem 1(b), the probability-generating function of T ˜ m is
G T ˜ m ( z ) = G ( 1 , , 1 , z ) = z m + 1 D m y m j = 1 m 1 1 ( 1 y j ) z d q y ( m ) , | z | 1 .
If | z | < 1 / 2 , then | z / ( 1 z ) | < 1 . Therefore, if 0 < y j < 1 , i = 1 , , m then we have that y j z 1 z < 1 , and the following expansion holds:
j = 1 m 1 1 ( 1 y j ) z = ( 1 z ) m j = 1 m 1 + y j z 1 z 1 = ( 1 z ) m j = 1 m α j 0 ( 1 ) α j y j z 1 z α j = k = 0 ( 1 ) k α 1 , , α m 0 α 1 + α m = k j = 1 m y j α j z k ( 1 z ) m k .
On the other hand, it is not difficult to check that
[ z n m 1 ] z k ( 1 z ) m k = n 2 m + k 1 ,
for n 2 , m 1 and 0 k n m 1 .
Thus, for n m + 1 , using (30), (31), Lemma 3 and (32),
P ( T ˜ m = n ) = [ z n ] G T ˜ m ( z ) = [ z n m 1 ] D m y m j = 1 m 1 1 ( 1 y j ) z d q y ( m ) = k = 0 ( 1 ) k α 1 , , α m 0 α 1 + α m = k D m y m j = 1 m y j α j d q y ( m ) [ z n m 1 ] z k ( 1 z ) m k = k = 0 n m 1 n 2 m + k 1 s q ( m + k 1 , m 1 ; 2 ) m + k + 1 q ! = = m 1 n 2 n 2 s q ( , m 1 ; 2 ) + 2 q ! .
(b)
For n m + 1 , using Theorem 2(a),
P ( T ˜ m n ) = k = m + 1 n P ( T ˜ m = k ) = k = m + 1 n = m 1 k 2 k 2 s q ( , m 1 ; 2 ) + 2 q ! = = m 1 n 2 k = + 2 n k 2 s q ( , m 1 ; 2 ) + 2 q ! = = m 1 n 2 n 1 + 1 s q ( , m 1 ; 2 ) + 2 q ! .
 □
Some other statistics related to weak record indicators can be expressed in terms of q-integrals and q-numbers as shown below.
Theorem 3.
For n 2 ,
(a) 
P ( I ˜ n = 1 ) = 1 q = 0 n 1 n 1 ( 1 ) + 1 q .
(b) 
P ( I ˜ 2 = 1 , , I ˜ n = 1 ) = 1 n q ! .
(c) 
P ( I ˜ 2 = 0 , , I ˜ n = 0 ) = = 0 n 1 n 1 ( 1 ) + 1 q .
Proof. 
(a) For n 2 , firstly observe that
0 = ( 1 1 ) n 1 = = 0 n 1 n 1 ( 1 ) = 1 1 q = 0 n 1 n 1 ( 1 ) 1 q + 1 + 1 q
Then,
= 0 n 1 n 1 ( 1 ) 1 + 1 q = q = 0 n 1 n 1 ( 1 ) q + 1 q .
The event { I ˜ n = 1 } consists of all the trajectories of the process X with a weak record at position n, i.e., trajectories of the form ( x ) n 1 x ; then, using the definition of the Jackson integral (see (5) and (33))
P ( I ˜ n = 1 ) = x 0 ( 1 q x + 1 ) n 1 q x p = 0 1 ( 1 q y ) n 1 d q y = = 0 n 1 ( 1 ) q 0 1 y d q y = = 0 n 1 n 1 ( 1 ) q + 1 q = 1 q = 0 n 1 n 1 ( 1 ) + 1 q .
(b) The trajectories corresponding to the event { I ˜ 1 = 1 , I ˜ 2 = 1 , , I ˜ n = 1 } are of the form ( x 1 ) ( x 2 ) ( x n ) , with 0 x 1 x 2 x n < + . Then, using Lemma 2 and (16),
P ( I ˜ 1 = 1 , I ˜ 2 = 1 , , I ˜ n = 1 ) = 0 x 1 x 2 x n < + j = 1 n q x j p = D n d q y ( n ) = 1 n q ! .
(c) The trajectories corresponding to the event { I ˜ 2 = 0 , , I ˜ n = 0 } are of the form ( x ) ( < x ) n 1 , with 1 x < + . Then
P ( I ˜ 2 = 0 , , I ˜ n = 0 ) = x 1 q x p ( 1 q x ) n 1 = q 0 q p ( 1 q + 1 ) n 1 = q 0 1 ( 1 q y ) n 1 d q y = = 0 n 1 n 1 ( 1 ) + 1 q ,
(see the proof of (a) for the last step).  □
The weak record indicators from geometric parents, as opposed to the continuous case, are not independentl for instance,
P ( I ˜ 2 = 1 , I ˜ 3 = 1 ) = 1 ( 1 + q ) ( 1 + q + q 2 ) 1 1 + q × 1 + q 2 ( 1 + q ) ( 1 + q + q 2 ) = P ( I ˜ 2 = 1 ) × P ( I ˜ 3 = 1 ) .
The number of weak records in a sample of size n is N ˜ n = m = 1 n I ˜ m (note that the 0-th record, R 0 = X 1 , is counted).
Theorem 4.
For n 1 ,
(a) 
P ( N ˜ n m ) = = m 1 n 1 n 1 s q ( 1 , m 2 ; 2 ) + 1 q ! , m = 2 , , n .
(b) 
P ( N ˜ n = m ) = = m n n 1 1 s q ( , m ) q ! , m = 1 , , n .
(c) 
E N ˜ n = = 1 n n ( 1 ) 1 q 1 q .
Proof. 
(a)
Observe that N ˜ n m if T ˜ m 1 n , m = 2 , , n . Then, the result is immediate using Theorem 2(b).
(b)
For m = 2 , , n , using (a) and the properties (22) and (23)
P ( N ˜ n = m ) = P ( N ˜ n m ) P ( N ˜ n m + 1 ) = j = m 1 n 1 n 1 j s q ( j 1 , m 2 ; 2 ) j + 1 q ! j = m n 1 n 1 j s q ( j 1 , m 1 ; 2 ) j + 1 q ! = n 1 m 1 s q ( m 2 , m 2 ; 2 ) m q ! + j = m n 1 n 1 j s q ( j 1 , m 2 ; 2 ) s q ( j 1 , m 1 ; 2 ) j + 1 q ! = n 1 m 1 s q ( m 1 , m 1 ; 1 ) m q ! + j = m n 1 n 1 j s q ( j , m 1 ; 1 ) j + 1 q ! = = j = m 1 n 1 n 1 j s q ( j , m 1 ; 1 ) j + 1 q ! = = m n n 1 1 s q ( , m ) q ! .
For m = 1 , note that if in a sample of size n, there is only one weak record. Then, it is in the first position, so, using Theorem 3(c),
P ( N ˜ n = 1 ) = P ( I ˜ 2 = 0 , , I ˜ n = 0 ) = j = 0 n 1 n 1 j ( 1 ) j j + 1 q = = 1 n n 1 s q ( , 1 ) q ! ,
where in the last equality we have used the property s q ( , 1 ) = ( 1 ) 1 1 q ! .
(c)
In the proof of Theorem 3(a), we obtained the integral representation
P ( I ˜ m = 1 ) = 0 1 ( 1 q y ) m 1 d q y , m 2 .
Clearly, this formula is also valid for m = 1 . Therefore,
E N ˜ n = m = 1 n E I ˜ m = m = 1 n P ( I ˜ m = 1 ) = 0 1 m = 1 n ( 1 q y ) m 1 d q y = 0 1 1 ( 1 q y ) n q y d q y = 0 1 = 1 n n ( 1 ) 1 ( q y ) 1 d q y = = 1 n n ( 1 ) 1 q 1 q .
 □
The inter-(weak)-record times are defined as S ˜ m = T ˜ m T ˜ m 1 , m 1 .
Theorem 5.
For m 1 ,
(a) 
The joint pmf of the inter-record times from a Ge(p) parent is
P S ˜ 1 = s 1 , S ˜ 2 = s 2 , , S ˜ m = s m = D m y m j = 1 m ( 1 y j ) s j 1 d q y ( m )
for s j 1 , j = 1 , , m .
(b) 
The pmf of the m-th inter-weak record is
P ( S ˜ m = s ) = = 0 s 1 s 1 ( 1 ) + 2 q m , s 1 .
(c) 
P ( S ˜ m > s ) = [ 0 , 1 ] m ( 1 y 1 y m ) s d q y ( m ) , s 1
Proof. 
(a) It is straightforward from Theorem 1(a).
(b) For s 1 (integer), from (34), using Lemma 4(a),
P ( S ˜ m = s ) = s 1 , , s m 1 1 P S ˜ 1 = s 1 , S ˜ 2 = s 2 , , S ˜ m = s = s 1 , , s m 1 1 D m y m ( 1 y m ) s 1 j = 1 m 1 ( 1 y j ) s j 1 d q y ( m ) = D m y m ( 1 y m ) s 1 y 1 y m 1 d q y ( m ) = = 0 s 1 s 1 ( 1 ) D m y 1 1 y m 1 1 y m + 1 d q y ( m ) = = 0 s 1 s 1 ( 1 ) + 2 q m ,
where the last equality follows from (18). (c) Using (36) and (12), for s 1 ,
P ( S ˜ m > s ) = t > s P ( S ˜ m = t ) = t > s D m y m ( 1 y m ) t 1 y 1 y m 1 d q y ( m ) = D m ( 1 y m ) s y 1 y m 1 d q y ( m ) = [ 0 , 1 ] m ( 1 y 1 y m ) s d q y ( m ) .
 □
The next result proves the asymptotic log-normality of the inter-(weak) record times from geometric parents.
Theorem 6.
The inter-weak record times from a Ge ( p ) parent, p ( 0 , 1 ) , satisfy
log S ˜ m + m q p log q m q p log q d Z N ( 0 , 1 ) , as m .
Proof. 
Let X i Ge ( p ) , i 1 , be independent random variables. The CLT yields,
Z m : = i = 1 m X i m q p m q p d Z N ( 0 , 1 ) , as m .
For x R , let us define v m ( x ) : = m q p + x m q p . From Theorem 5(c) and (26),
P ( S ˜ m > s ) = E 1 q i = 1 m X i s = E 1 q v ( Z m ) s
for s = 0 , 1 , .
Let us define the functions g m ( z , t ) : = ( 1 q v m ( z ) ) 1 / q v m ( t ) , z , t R , m 1 . These functions satisfy the following properties:
(a)
lim m + g m ( z , t ) = 0 , if z < t e 1 , if z = t 1 , if z > t
(b)
For a fixed value of t > m q , the function z ( m q , + ) g m ( z , t ) is non-decreasing.
Then, for t R ,
P log S ˜ m + m q p log q m q p log q > t = P S ˜ m > q v m ( t ) = E 1 q v m ( Z m ) q v m ( t ) = E g m ( Z m , t ) = k 0 g m ( z m , k , t ) P ( Z m = z m , k ) ,
where z m , k = ( k m q / p ) / m q / p .
Let Φ be the cumulative distribution function of a standard normal random variable. If we show that
lim m + E I ( Z m > t ) g m ( Z m , t ) = 0 , t R
then we can conclude
lim m + E g m ( Z m , t ) = lim m + P ( Z m > t ) = 1 Φ ( t ) , t R ,
and this proves (37).
In order to prove (38), for η > 0 , we have
E I ( Z m > t ) g m ( Z m , t ) = k 0 I ( z m , k > t ) g m ( z m , k , t ) P ( Z m = z m , k ) = k 0 z m , k > t + η + k 0 z m , k < t η + k 0 | z m , k t | η I ( z m , k > t ) g m ( z m , k , t ) P ( Z m = z m , k ) = ( I ) + ( I I ) + ( I I I ) ,
where
( I ) = k 0 z m , k > t + η 1 g m ( z m , k , t ) P ( Z m = z m , k ) 1 g m ( t + η , t ) P ( Z m > t + η ) 1 g m ( t + η , t ) 0 , as m + ; ( I I ) = k 0 z m , k < t η g m ( z m , k , t ) P ( Z m = z m , k ) g m ( t η , t ) P ( Z m < t η ) g m ( t η , t ) 0 , as m + ; ( I I I ) = k 0 | z m , k t | η I ( z m , k > t ) g m ( z m , k , t ) P ( Z m = z m , k ) 2 P ( t η Z m t + η ) m 2 Φ ( t + η ) Φ ( t η ) ,
and this quantity can be made arbitrarily small choosing an η > 0 yet small enough.  □
A by-product of Theorem 6 is the following weak law of large numbers:
S ˜ m 1 / m P q q p , m .
When q 1 , q q / p e , and this is the limit in probability for Δ m 1 / m as m + , where Δ m denotes the waiting-time between the ( m 1 ) -th and the m-th upper record in a sequence of iid random variables with a continuous distribution; see [6]. Then, (39) can be seen as a q-analogue of Theorem 1 in [6]. Similarly, our Theorem 6 is a q-analogue of Theorem 2 in [6].

4. Weak Records from Geometric Parents Are q -Analogues of Records from Absolutely Continuous Parents

In this section, we show that the distributions and statistics related to weak records from a geometric parent, Ge ( p ) are q-analogues ( q = 1 p ) of the corresponding distributions and statistics related to ordinary records from absolutely continuous parents.
Theorem 7.
The joint pmf of T ˜ 1 , , T ˜ m from a Ge ( p ) parent given in (27) is a q-analogue of the joint pmf of ( T 1 * , , T m * ) from an absolutely continuous parent given in (2).
Proof. 
Without loss of generality, suppose that ( T 1 * , , T m * ) are the record times of a sequence of iid random variables uniformly distributed in the interval ( 0 , 1 ) . For integers t 0 = 1 < t 1 < < t m , the event { T 1 * = t 1 , , T m * = t m } corresponds to trajectories of length t m m + 1 of the form
( x 1 ) ( < x 1 ) t 1 2 ( x 2 ) ( < x 2 ) t 2 t 1 1 ( x m ) ( < x m ) t m t m 1 1 ( x m + 1 ) ,
with 0 < x 1 < x 2 < < x m + 1 < 1 . Then,
P ( T 1 * = t 1 , , T m * = t m ) = 0 1 x 1 1 x m 1 j = 1 m x j t j t j 1 1 d x m + 1 d x 2 d x 1 = 0 1 x 1 1 x m 1 1 ( 1 x m ) j = 1 m x j t j t j 1 1 d x m d x 2 d x 1 = 0 1 0 y 1 0 y m 1 y m j = 1 m ( 1 y j ) t j t j 1 1 d y m d y 2 d y 1 ,
where the last step is the result of the change in variables y j = 1 x j , j = 1 , , m .
Finally, the result follows by comparing (40) with the q-integral expression given in (27) and by using Lemma 5.  □
The addition of q-analogues is a q-analogue of the corresponding sum. As an application of this property, as a corollary of Theorem 7, we have that the pmf of T ˜ m (given in Theorem 2) is a q-analogue of (3). Clearly, the pmf of inter-weak-record times from geometric parents (see (35)) is a q-analogue of (4). Similarly, as q 1 , it can be checked that
P ( I ˜ n = 1 ) = 1 q = 0 n 1 n 1 ( 1 ) + 1 q 1 n P ( I ˜ 2 = 1 , , I ˜ n = 1 ) = 1 n q ! 1 n ! . P ( N ˜ n = m ) = = m n n 1 1 s q ( , m ) q ! | s n , m | n ! . E N ˜ n = = 1 n n ( 1 ) 1 q 1 q j = 1 n 1 j ,
so all of these quantities in the case of a geometric parent are q-analogues of the corresponding quantities in the absolutely continuous case.

5. Ordinary Records from Geometric Parents

Following similar arguments as in the previous sections, it is possible to derive results for ordinary records from a Ge ( p ) parent, p ( 0 , 1 ) . As before, let q = 1 p and Q = 1 / q . We quote without proof the following results.
Joint distribution of (ordinary) record times:
P T 1 = t 1 , T 2 = t 2 , , T m = t m = D m , q q y m j = 0 m ( 1 q y j ) t j t j 1 1 d q y ( m ) ,
for m 1 and 1 = t 0 < t 1 < < t m . The random vector ( T 1 , , T m ) has a probability generating function
G ( z 1 , , z m ) = D m , q q y m j = 1 m z j j + 1 1 ( 1 q y j ) z j z m d q y ( m ) .
Probability function of T m :
P ( T m = n ) = = m 1 n 2 n 2 s Q ( , m 1 ; 2 ) + 2 Q !
for m 1 and n m + 1 .
Joint distribution of inter-record times:
P S 1 = s 1 , S 2 = s 2 , , S m = s m = D m , q j = 0 m ( 1 q y j ) s j d q y ( m ) ,
for s j 0 , j = 1 , , m .
Probability function of inter-record times:
P ( S m = s ) = = 0 s 1 s 1 ( 1 ) + 2 Q m , s 1 .
for m 1 , s 0 .
Record indicators:
P ( I n = 1 ) = 1 Q = 0 n 1 n 1 ( 1 ) + 1 Q .
P ( I 2 = 1 , , I n = 1 ) = 1 n Q ! .
P ( I 2 = 0 , , I n = 0 ) = = 0 n 1 n 1 ( 1 ) + 1 Q .
Number of ordinary records:
P ( N n = m ) = = m n n 1 1 s Q ( , m ) Q ! , m = 1 , , n , P ( N n m ) = = m 1 n 1 n 1 s Q ( 1 , m 2 ; 2 ) + 1 Q ! m = 1 , , n . E N n = = 1 n n ( 1 ) 1 Q 1 Q .
Comparing these results with those in Section 3, we observe that the formulas for ordinary records are obtained by changing q by Q : = 1 / q in the formulas for weak records. Additionally, as in the case of weak records, the above quantities are Q-analogues of the corresponding ones in the absolutely continuous case.

6. Conclusions

In the case of a geometric parent, we have obtained representations of the joint distribution of weak and ordinary record times and inter-record times in terms of Jackson integrals. These representations are useful for obtaining the marginal distribution of record times, inter-record times, record indicators, number of records, etc. In particular, the marginal distribution of record times can be expressed in terms of non-central q-Stirling numbers. Additionally, we show that statistics related to weak records from geometric parents, Ge ( p ) , behave as q-analogues ( q = 1 p ) of the corresponding statistics in the case of absolutely continuous parents. There is also a duality between weak and ordinary records from geometric parents: it is possible to obtain results from ordinary records by changing q by Q = 1 / q in the corresponding formulas or results for weak records.

Author Contributions

These authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

Grant PID2020-114594GB-C21 funded by MCIN/AEI 10.13039/50110001033.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Lemma 1.
(a) By changing the summation variables, 1 = x 1 , 2 = x 2 x 1 , …, m = x m x m 1 , then
( 1 q ) m 0 x 1 x m < f ( q x 1 , q x 2 , , q x m ) q i = 1 m x i = ( 1 q ) m 1 , , m 0 f ( q 1 , q 1 + 2 , , q 1 + , m ) j = 1 m q j m j i = 1 m q i = [ 0 , 1 ] m f ( y 1 , y 1 y 2 , , y 1 y m ) j = 1 m y j m j d q y π 1 d q y π m
where the last equality is justified by (8).
(b) By using the change of variables 1 = x 1 , 2 = x 2 x 1 1 , …, m = x m x m 1 1 , we obtain
( 1 q ) m 0 x 1 < < x m < f ( q x 1 , q x 2 , , q x m ) q i = 1 m x i = q m 2 ( 1 q ) m 1 , , m 0 f ( q 1 , q 1 + 2 + 1 , , q 1 + , m + ( m 1 ) ) j = 1 m q j m j i = 1 m q i = q m 2 [ 0 , 1 ] m f ( y 1 , q y 1 y 2 , , q m 1 y 1 y m ) j = 1 m y j m j d q y π 1 d q y π m .
 □
Proof of Lemma 2.
(a) If m = 1 , the result is immediate by using the definition of Jackson integral given in (5). Thus, let us suppose in the following that m > 1 . Let us define
A m ( x 1 , , x m 1 ) : = ( 1 q ) x m = x m 1 f ( q x 1 , , q x m ) q x m , A j ( x 1 , x j 1 ) : = ( 1 q ) x j = x j 1 q x j A j + 1 ( x 1 , , x j ) , j = m 1 , , 2 A 1 : = ( 1 q ) x 1 = 0 q x 1 A 2 ( x 1 ) .
First of all, observe that by construction, A 1 coincides with the left-hand side of the statement in Lemma 2 ( a ) . On the other hand, the definition of the Jackson integral yields
A m ( x 1 , , x m 1 ) = ( 1 q ) q x m 1 = 0 f ( q x 1 , , q x m 1 , q x m 1 q ) q = 0 q x m 1 f ( q x 1 , , q x m 1 , y m ) d q y m ,
and by (back) recurrence, it is not difficult to show that
A j ( x 1 , x j 1 ) = 0 q x j 1 0 y j 0 y m 1 f ( q x 1 , , q x j 1 , y j , , y m ) d q y m d q y j ,
for j = m , , 2 . Finally, from (A1) and (A2),
A 1 = ( 1 q ) x 1 = 0 q x 1 0 q x 1 0 y j 0 y m 1 f ( q x 1 , y 2 , , y m ) d q y m d q y 2 = 0 1 0 y 1 0 y m 1 f ( y 1 , , y m ) d q y ( m )
and the result follows from (10).
(b) The proof is similar to the one in ( a ) , so we omit it.  □
Proof of Lemma 3.
We only prove (a). The proof of (b) is similar. For m 1 and k = 0 , using (14), the change of variables c j = α j + + α m + m j + 2 , j = 1 , , m 1 and (24), we obtain
α 1 , , α m 0 α 1 + α m = k D m y m j = 1 m y j α j d q y ( m ) = α 1 , , α m 0 α 1 + α m = k j = 1 m 1 α j + + α m + m j + 2 q = 1 m + k + 1 q 2 c 1 < < c m 1 m + k i = 1 m 1 1 c j q = ( 1 ) k s q ( m + k 1 , m 1 ; 2 ) m + k + 1 q !
where the last equality follows from (24).  □
Proof of Lemma 4.
For simplicity, we only give the proof in the case that m = 2 ; for larger values the proof is similar but tedious. Thus, we have
s 1 1 0 1 0 y 1 f ( y 2 ) ( 1 y 1 ) s 1 1 d q y 2 d q y 1 = s 1 1 0 1 ( 1 y 1 ) s 1 1 0 y 1 f ( y 2 ) d q y 2 d q y 1 = s 1 1 ( 1 q ) = 0 ( 1 q ) s 1 1 0 q f ( y 2 ) d q y 2 q = ( 1 q ) = 0 0 q f ( y 2 ) d q y 2 q s 1 1 ( 1 q ) s 1 1 = ( 1 q ) 0 1 f ( y 2 ) d q y 2 + = 1 0 q f ( y 2 ) d q y 2 q s 1 1 ( 1 q ) s 1 1 = ( 1 q ) 0 1 f ( y 2 ) d q y 2 + = 1 0 q f ( y 2 ) d q y 2 = ( 1 q ) = 0 0 q f ( y 2 ) d q y 2 = ( 1 q ) = 0 1 q 0 q f ( y 2 ) d q ( y 2 ) q = 0 1 1 y 1 0 y 1 f ( y 2 ) d q y 2 d q y 1 = 0 1 0 y 1 f ( y 2 ) y 1 d q y 2 d q y 1 .
 □

References

  1. Euler, L. Introductio in Analysin Infinitorum. Tomus Primus; Sociedad Andaluza de Educación Matemática “Thales”, Seville; Real Sociedad Matemática Española: Madrid, Spain, 2000. [Google Scholar]
  2. Ernst, T. A Comprehensive Treatment of q-Calculus; Birkhäuser/Springer Basel AG: Basel, Switzerland, 2012. [Google Scholar]
  3. Chandler, K.N. The distribution and frequency of record values. J. R. Stat. Soc. Ser. B 1952, 14, 220–228. [Google Scholar] [CrossRef]
  4. Arnold, B.C.; Balakrishnan, N.; Nagaraja, H.N. Records; Wiley Series in Probability and Statistics; John Wiley & Sons, Inc.: New York, NY, USA, 1998. [Google Scholar]
  5. Nevzorov, V.B. Records: Mathematical Theory; Translations of Mathematical Monographs; American Mathematical Society: Providence, RI, USA, 2001; Volume 194. [Google Scholar]
  6. Neuts, M.F. Waitingtimes between record observations. J. Appl. Probab. 1967, 4, 206–208. [Google Scholar] [CrossRef]
  7. Kirschenhofer, P.; Prodinger, H. The path length of random skip lists. Acta Inf. 1994, 31, 775–792. [Google Scholar] [CrossRef] [Green Version]
  8. Kirschenhofer, P.; Martínez, C.; Prodinger, H. Analysis of an optimized search algorithm for skip lists. Theor. Comput. Sci. 1995, 144, 199–220. [Google Scholar] [CrossRef] [Green Version]
  9. Prodinger, H. Combinatorics of geometrically distributed random variables: Left-to-right maxima. Discret. Math. 1996, 153, 253–270. [Google Scholar] [CrossRef] [Green Version]
  10. Bai, Z.D.; Hwang, H.K.; Liang, W.Q. Normal approximations of the number of records in geometrically distributed random variables. Random Struct. Algorithms 1998, 13, 319–334. [Google Scholar] [CrossRef]
  11. Gouet, R.; López, F.; Sanz, G. Central limit theorems for the number of records in discrete models. Adv. Appl. Probab. 2005, 37, 781–800. [Google Scholar] [CrossRef]
  12. Gouet, R.; López, F.; San Miguel, M. A martingale approach to strong convergence of the number of records. Adv. Appl. Probab. 2001, 33, 864–873. [Google Scholar] [CrossRef]
  13. Key, E. On the number of records in an iid discrete sequence. J. Theor. Probab. 2005, 18, 99–107. [Google Scholar] [CrossRef]
  14. Kim, J.; Stanton, D. On q-integrals over order polytopes. Adv. Math. 2017, 308, 1269–1317. [Google Scholar] [CrossRef] [Green Version]
  15. Charalambides, C. Non-central generalized q-factorial coefficients and q-Stirling numbers. Discret. Math. 2004, 275, 67–85. [Google Scholar] [CrossRef] [Green Version]
  16. Kac, V.; Cheung, P. Quantum Calculus; Universitext; Springer: New York, NY, USA, 2002. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

López-Blázquez, F.; Salamanca-Miño, B. Multiple q-Integral and Records from Geometrically Distributed Sequences. Mathematics 2022, 10, 2313. https://doi.org/10.3390/math10132313

AMA Style

López-Blázquez F, Salamanca-Miño B. Multiple q-Integral and Records from Geometrically Distributed Sequences. Mathematics. 2022; 10(13):2313. https://doi.org/10.3390/math10132313

Chicago/Turabian Style

López-Blázquez, Fernando, and Begoña Salamanca-Miño. 2022. "Multiple q-Integral and Records from Geometrically Distributed Sequences" Mathematics 10, no. 13: 2313. https://doi.org/10.3390/math10132313

APA Style

López-Blázquez, F., & Salamanca-Miño, B. (2022). Multiple q-Integral and Records from Geometrically Distributed Sequences. Mathematics, 10(13), 2313. https://doi.org/10.3390/math10132313

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop