1. Introduction
We develop a two-parameter version of an arc-sine law for a last hitting time, and this is achieved in the two-parameter Wiener space. If a random variable
X has its cumulative distribution
, then
X is called arcsine-distributed. It is well known that Lévy proved an arc-sine theorem for Brownian motion [
1] in 1939 and Erdös and Kac extended Levy’s result to general random variables [
2] in 1947. Since then a number of similar laws for various processes have been introduced by researchers, most of which are studied in the Wiener space, for example, Reference [
3]. One recent work [
4] studies the arc-sine law in the analogue of the Wiener space; the analogue of the Wiener space is introduced as ’a kind of the generalisation of the Wiener space’ in Reference [
5].
Some phenomena in many fields of science have been explained by arc-sine distributions; Barato et al. applied to thermodynamics [
6]. Arc-sine laws play crucial roles in problems arising in dynamic systems using random walks. Due to the aspect of the U-shaped density
of an arc-sine distribution, arc-sine laws notably have been applied in the fields of mathematical finance. Dale and Workman discuss the law in relation to the treasury bill future market [
7].
Among many arc-sine laws, we are interested in the process called a last hitting time in this paper, which is also called the second arc-sine law (In some literature, such as Reference [
8], the name may indicate a slightly different event). In order to state it we assume that
is the Wiener process defined on the Wiener space
. Let
be the last time the process hits zero. Then the probability
[
9].
The existing arc-sine laws are about a stochastic process with one parameter t, which is a time parameter in most applications. If there is another varying key factor of an event described by the process, then we need to consider another parameter besides t. That is, we need a system of random variables with two parameters, say . In this paper we challenge to develop such an idea. A two-parameter process in relation to our interest is very rare in literature. To the best of the author’s knowledge, arc-sine laws have not been studied yet in two-parameter processes. We devise a two-parameter process analogous to a last hitting time and develop the two parameter version of the arc-sine law using the properties of the two-parameter Wiener process. Though the probability of an event described by is much more complicated than that of , a new arc-sine law has been accomplished based on our development.
In the next section we briefly introduce the two-parameter Wiener space for background knowledge of our development.
Section 3 provides the procedure of how to construct a probability density function (pdf) for the stochastic process of our interest. Using the pdf, the arc-sine law for last hitting points in the the two-parameter Wiener space is calculated in
Section 4. Finally conclusions are given in
Section 5.
2. Background: The Two-Parameter Wiener Process
In this section we briefly introduce the two-parameter Wiener process for background knowledge of our development as it is far less accessible than the Wiener process. This section is mostly based on Reference [
10]. For readers who are not familiar with the two-parameter Wiener process, a more detailed digest of it can be found in the author’s recent work [
11].
Given a set in , will denote the space of real-valued continuous functions on R. We are interested in the subset of defined below.
Definition 1. Let in . We denote the set of continuous functions by Definition 2 (Yeh)
. Let in , and , ,where , , and , .- (1)
I is called a strict interval in .
- (2)
If D is a measurable set in , then I is called a cylinder.
- (3)
, ..., are called restriction points of I.
- (4)
= the collection of all strict intervals I. (Then is a semi-algebra.)
Throughout the paper, we denote an exponential function by , usually for a random variable; we use it and interchangeably.
Definition 3 (Yeh)
. For , the two-parameter Wiener measure of I is defined bywhereHere, for , .
The completed measure space of the two-parameter Wiener measure
is denoted by
and called the two-parameter Wiener space. The space is also called the Yeh-Wiener space in some literature because the space was first introduced by Yeh in Reference [
10] and further work [
12]. One can see the measure space is a probability space since
.
Proposition 1 (Yeh)
. Let . Let be a Lebesgue measurable function on and , then F is two-parameter Wiener measurable andThe equality means that if one side exists then the other side also exists and both sides coincide. It tells us that the integral over the infinite dimensional function space on the left is reduced to an ordinary Lebesgue integral over the finite dimensional space on the right [10]. Example 1. - (a)
- (b)
Here, denotes the expectation of a random variable.
Now we introduce the two-parameter Wiener process defined on to be used in the next sections.
Definition 4 (Yeh).
Let and is defined on the two-parameter Wiener space byThen is called the two-parameter Wiener process.
From Example 1, we can see that follows a normal distribution . Obviously, for all , for all t and for all s as .
The function
is a random variable and the collection
was originally called a two-time parameter Wiener process by Yeh [
10]; later, it is called the Yeh-Wiener process in the literature, such as References [
13,
14,
15]. We keep the name, the two-parameter Wiener process, throughout this paper.
By the definition of in Definition 3
Since
is a probability measure, we can express it by
P:
In the later sections we will use the following well-known integral in our calculations.
Proposition 2. For a positive real number a,
3. Construction of a PDF for , a Two-Parameter Version of a First Hitting Time
In this section, we aim to build a probability density function (pdf) of an event , where c is a real number.
Let x be in , that is, x is continuous on and for all and . A one-parameter process is called the first point that hits c or a first hitting time to c, and its distribution is well-known.
We define the two-parameter version of denoted by for below. For this we adapt the inequality < for one-dimensional to two-dimensional .
Definition 5. Let and are two points in the rectangle R.
- (1)
We define if and .
We also define if either or [ and ].
- (2)
Similarly, we define if either or .
We also define if either or [ and ].
In Definition 5, if R is divided into four subrectangles with common vertex , then the symbol < is for the description of the bottom left subrectangle and the symbol is used for the description of the complement of the top right subrectangle.
Definition 6. Let c be a real number. is the set of points satisfying:
- (i)
and
- (ii)
for any with , but .
For one-variable functions , the first hitting point is a single point in the domain as introduced at the beginning of this section. On the other hand, for two-variable functions , the set is a curve in general. If and , then never hits c. Hence, can be called a least point hitting c in some sense or the first point in a radial direction from .
As introduced in
Section 2,
is the two-parameter Wiener measure defined on
, where
R is a rectangle
. When a specific point
in
R is chosen, our interest is the region beyond the point
; that is, we are interested in the subrectangle
of
R. If we regard the set of
x’s (surfaces) restricted on this subrectangle as a new space,
is insufficient to configure the set of
x’s defined on the subrectangle; we cannot apply the two-parameter Wiener measure
directly. The reason is that
or
on the boundary of the subrectangle may not be zero as required for
. This is one of the greatest difficulties in the development of a process in the two-parameter Wiener space compared to the Wiener space (For the Wiener space, if a specific point
is chosen, then curves
x restricted on its subinterval
can be regarded as a new Wiener space by translating
x by
so that we can make the translated
x has value zero at
t. Then the Wiener measure can be applied directly to a set of those
x’s defined on the subinterval).
Remark 1. Let and . Then we have the following relationships of the two-parameter Wiener measure (see Definition 3). Being motivated by the relationships in the remark above, we introduce a restriction of on a subrectangle of R as follows.
Definition 7. Let , then is a subrectangle of R. We denote the subrectangle by ; . Let for .
Define . Let b be a real number and define .
Definition 8. Let us define a measure w on , satisfying the followings. For a real number α and a point , We will calculate the measure w when , a specific value of at , that is, and denote it by in the next theorem.
Theorem 1. If is the measure when in Definition 8 (here, b is the value of at , that is, ), then Proof. Each exponent of exponential functions
in the integrand of
is a quadratic function of an integral variable. Using completing the square on a quadratic function with respect to
, the exponents of
and
can be expressed by
where
. Then the inner integral has the form of
for some constants a and d. Since
as in Proposition 2, the inner integral with respect to
is calculated; then a, d and the front
are functions of two variables
and
. Following the same technique for
as done for
, the middle integral is calculated too and becomes a function of
, which is the variable of the outer integral.
Though the process of these calculations are tremendously lengthy, the resulting integral has been compactly expressed as a single exponential function (In the author’s calculations, five sheets of A4 paper are filled with equations written in small-sized letters; the width of a page was not long enough for a line even in landscape orientation). Therefore, the triple integral can be reduced to a single integral with respect to
:
This is the Equation (
5). □
Lemma 1. For a real number and , Proof. We use a trivial set relationship of an event, , Since the last set among the three sets on the right side has two-parameter Wiener measure (or probability) zero, we do not consider it. The rest of the two sets have an equal probability. Thus we examine only one of them.
Obviously, and we need to show the converse inclusion. If is contained in the set on the right side, that is, , then cannot be and . This means that is contained in the set on the left side. Therefore, the sets on the both sides are equivalent and they have the same measure as desired. □
For , we have a similar result to Lemma 1 and state in the next corollary.
Corollary 1. For and , Using Lemma 1, we can express the measure of an event of as a single integral.
Theorem 2. Let . For and ,For and , . Proof. We know that
by Lemma 1. Using Equation (
5) for the set on the right, we have
Here, the last equation is derived by the change of variables . Therefore, . The second statement for follows from the exactly same manner as the case and Corollary 1. □
Remark 2. Let and .
- (1)
Suppose that .
Let , , and . (See Definition 5 for the symbol )
Then and - (2)
By (1), is a constant function of associated with b; so and . Therefore, - (3)
In the case of or , we can obtain the equality (6) using the same way as (1) and (2) above.
Theorem 3. For , Proof. Let
as in Remark 2. First we will show that
From Theorem 2,
. In order to get
, we use the Fundamental Theorem of Calculus; the integral is a function of
which appears only on an integral boundary. Then the same technique is applied to the second partial derivative
.
Next we will use Equation (
7) with
and
for both cases
and
. At the beginning of the proof we regarded
as a cumulative distribution function of the event
; that is,
is known to be the pdf of F. We can compute
using the pdf as follows:
and we obtain the desired result. □
The pdf in the proof of Theorem 3 will play a key role in our main theorem.
4. The Two-Parameter Version of an Arc-Sine Law for a Hitting Time
Introducing the one-parameter case for the Wiener process
, a last hitting time
has the probability
([
9]), which is known to be an arc-sine law.
We now define a process L in the two-parameter Wiener space that is the two parameter version of a last hitting time and aim to develop for our process L a law similar to the existing arc-sine law for . We keep the notations and . Differently from the one-dimensional feature of the parameter t, the word “last” for a two-dimensional parameter will have a weakened meaning. We define an analogous concept to the last hitting time in the one-dimension as follows.
Definition 9. Let be the set of points in a rectangle R with the following properties:
- (i)
and for any , ;
- (ii)
on the boundary or of R, if for any and , there exists a point hitting zero in the region , then we define .
Here, the inequality > used in (i) has been defined in Definition 5.
Remark 3. (1) If , then x never hits zero beyond the point . Hence, in some sense, the point can be called the last point in a radial direction from the origin . Therefore, we can call L the set of last hitting points of radial directions.
(2) Considering a degenerate case of the definition of L, if is an increasing function of t and is an increasing function of s, then we define .
Lemma 2. For , the two-parameter Wiener measure of L can be expressed by Proof. Let . Then the surface of x over R does not hit zero at , that is, , say . Also for any with , . For all , x does not hit zero on .
The subrectangle has been denoted by (Definition 7) and a new measure denoted by on (Definition 8) has been defined and used in Theorems 2 and 3. In the event discussed here, b in the theorems is replaced by . This means that x may hit zero after S for s and T for t; that is, . The event can be written as . Therefore, □
We provide the main theorem of our work.
Theorem 4. Proof. From Lemma 2,
Using Equation (
6) and Theorem 3 for
,
and so
is a triple integral of the form of
. The integrand of the triple integral consists of three terms (fractions are involved each); we calculate each term of the triple integral one by one.
here, the Fubini theorem is used in the first equality. The second equality comes from
for a positive real number a.
Here, the second equality comes from for a positive real number a.
Adding the three terms (i), (ii) and (iii) yields
Substituting
and
, so
and
, the integrand is replaced by
and the integral is simplified to
The inner integral of Equation (
8) (
will be multiplied later) with respect to u becomes
Set
for the outer integral of (
8), then
,
and
Finally, we have obtained that the probability of follows an arc-sine distribution as desired. □
5. Conclusions
We have developed a two-parameter process analogous to a last hitting time for the Wiener process and provided the probability distribution of the process. This can be called the two-parameter version of the arc-sine law for a last hitting time. Differently from the one-dimensional feature of the parameter t, we were cautious to say the word “the last time” for a two-dimensional parameter . We devised an analogous concept for the two-dimensional parameter, which can be called last points hitting b in radial directions, so our word “last” has a weakened meaning.
We have calculated the probability of the event of last points associated with the process and confirmed that it follows an arc-sine distribution. The resulting probability calculated in Theorem 4 shows a similar form to the arc-sine law for the last hitting time of the Wiener process. However, the derivation of our result is not as simple as the derivation for a one-parameter process. The probability of our two-parameter process is naturally led to a multiple integral and required some device on measures as well as very lengthy calculations and a series of substitutions. Nevertheless, we have obtained that it follows an arc-sine distribution.
Our result can be added to many existing arc-sine laws, in particular, uniquely for two-dimensional parameters. One can use our results to predict an event happened in a system of random variables with two parameters, which is not available among existing laws for one parameter processes. The result of our work can be adapted to other acrsine laws, such as the first or third laws, which is our ongoing work.