# Mathematical Model of Stock Prices via a Fractional Brownian Motion Model with Adaptive Parameters

Authors: Tidarut Areerak
Journal: ISRN Applied Mathematics
 (2014)

---

## Abstract

The paper presents a mathematical model of stock prices using a fractional Brownian motion model with adaptive parameters (FBMAP). The accuracy index of the proposed model is compared with the Brownian motion model with adaptive parameters (BMAP). The parameters in both models are adapted at any time. The ADVANC Info Service Public Company Limited (ADVANC) and Land and Houses Public Company Limited (LH) closed prices are concerned in the paper. The Brownian motion model with adaptive parameters (BMAP) and fractional Brownian motion model with adaptive parameters (FBMAP) are applied to identify ADVANC and LH closed prices. The simulation results show that the FBMAP is more suitable for forecasting the ADVANC and LH closed price than the BMAP.


---

## Body

## 1. Introduction

The ideas of using a Brownian motion process to explain the behavior of the risky asset prices were presented by Black et al. [1–3]. The stock prices presented in the paper are also the type in the risky asset prices. Therefore, the Brownian motion is usually used to model a stock price. However, Brownian motion process has the independent increments property. This means that the present price must not affect the future price. In fact, the present stock price may influence the price at some time in the future. Hence, Brownian motion process is not suitable to explain the stock price. Another process, a fractional Brownian motion process, exhibits a long range dependent property. Therefore, a fractional Brownian motion process can be used to describe the behavior of stock price instead of Brownian motion process.The rate of return and volatility in general asset pricing model are usually the constant parameters. Actually, the rate of return and volatility in the model are not constant at any time. In the paper, these parameters are updated depending on time by using the new information.The ADVANC Info Service Public Company Limited (ADVANC) and Land and Houses Public Company Limited (LH) stock prices are considered in the paper. These two stocks are chosen from different stock exchange of Thailand (SET) industry groups. The ADVANC and LH prices are selected from technology group (TECH) and property and construction group (PROPCON), respectively.The ADVANC and LH stock price models are studied by using fractional Brownian motion process to explain uncertainly behavior instead of Brownian motion process. The Brownian motion model with adaptive parameters (BMAP) and the fractional Brownian motion model with adaptive parameters (FBMAP) are presented to model the ADVANC and LH stock prices.The paper is organized as follows. Preliminaries on a fractional Brownian motion are given in Section2. The estimation of the rate of return and volatility is shown in Section 3. In Section 4, the BMAP and the FBMAP are explained. Finally, Section 5 concludes the work in the paper.


## 2. Preliminaries on a Fractional Brownian Motion Process

In the general accepted model, the randomness of stock price is modelled by Brownian motion process. A stock price process(
S
t
,
t
≥
0
) is represented by the stochastic differential equation (SDE) as shown in
(1)
d
S
t
=
S
t
(
μ
d
t
+
σ
d
W
t
)
.
Note that the parameters μ and σ are the rate of return and the volatility, respectively. The process (
W
t
,
t
≥
0
) in (1) is a standard Brownian motion process. The stochastic differential equation (1) is driven by the Brownian motion process (
W
t
,
t
≥
0
). In the real world, μ and σ in (1) are not constant at any time. Hence these parameters in the paper are the adaptable parameters based on time. In the paper, model (1) is called a Brownian motion model with adaptive parameters (BMAP).In practice, the dynamics of stock price have a long memory (long range dependence). The BMAP model in (1) is not suitable to describe the dynamics of stock price. Therefore, the fractional Brownian motion process is considered in the paper. The fractional Brownian motion process (
B
t
H
,
t
≥
0
) with Hurst index H is a centered Gaussian process. If H
=
0.5, then (
B
t
H
,
t
≥
0
) is a standard Brownian motion process. If H
≠
0.5, then (
B
t
H
,
t
≥
0
) is neither a semimartingale nor a Markov process. For H
≠
0.5 case, the (
B
t
H
,
t
≥
0
) is the long memory process. The (
B
t
H
,
t
≥
0
) is represented in (2) by Mandelbrot and Van Ness [4]. Consider the following:
(2)
B
t
H
=
1
Γ
(
1
+
α
)
[
Z
t
+
B
t
]
.
The function Γ
(
·
) is the gamma function. The process (
Z
t
,
t
≥
0
) is defined by Z
t
=
∫
-
∞
0
[
(
t
-
s
)
α
-
(
-
s
)
α
]
d
W
s. The process (
B
t
,
t
≥
0
) is described by B
t
=
∫
0
t
(
t
-
s
)
α
d
W
s. The parameter α
=
H
-
1
/
2, where H
∈
(
0,1
) and (
W
t
,
t
≥
0
) is a standard Brownian motion process.The rate of return and volatility are not constant at any time. Hence, the paper also proposes the new approach of the asset pricing model. In this case, the driving process of model (1) is replaced by a fractional Brownian motion process. The rate of return and volatility are adaptive parameters. In this case, the model can be represented by the stochastic differential equation (SDE) as shown in
(3)
d
S
t
=
S
t
(
μ
d
t
+
σ
d
B
t
H
)
.
The parameters μ and σ in (3) are the rate of return and the volatility, respectively. The μ and σ are adaptive parameters the same as the previous model. The (
B
t
H
,
t
≥
0
) is a fractional Brownian motion process. In the paper, model (3) is called a fractional Brownian motion model with adaptive parameters (FBMAP).Alòs et al. [5] have proposed to use the process (
B
t
,
t
≥
0
) instead of (
B
t
H
,
t
≥
0
), since (
Z
t
,
t
≥
0
) has absolutely continuous trajectory. So the process (
B
t
,
t
≥
0
) has long range dependence. Hence, the model (3) can be considered as shown in
(4)
d
S
t
=
S
t
(
μ
d
t
+
σ
d
B
t
)
.An approximate approach to stochastic differential equation perturbed by fractional Brownian motion was proposed by Thao [6]. The process (
B
t
ε
,
t
≥
0
) is introduced. For every ε
>
0, the process (
B
t
ε
,
t
≥
0
) is defined by
(5)
B
t
ε
=
∫
0
t
(
t
-
s
+
ε
)
α
d
W
s
.
The process (
B
t
ε
,
t
≥
0
) is a semimartingale. Therefore, this process can be written as in
(6)
B
t
ε
=
α
∫
0
t
φ
s
ε
d
s
+
ε
α
W
t
,
where φ
t
ε
=
∫
0
t
(
t
-
s
+
ε
)
α
-
1
d
W
s.The process(
B
t
ε
,
t
≥
0
) converges to (
B
t
,
t
≥
0
) in L
2
(
Ω
) when ε approaches to 0. This convergence is uniform with respect to t
∈
[
0
,
T
]. Hence, the model (4) can be considered as shown in
(7)
d
S
t
=
S
t
(
μ
d
t
+
σ
d
B
t
ε
)
.


## 3. The Estimation of the Rate of Return and Volatility

In this paper, the ADVANC and LH closed prices are identified by two asset pricing models. In the BMAP, the driving process is Brownian motion. On the other hand, the driving process is fractional Brownian motion in the FBMAP. The parametersμ and σ in both models are adaptive parameters at any time. The ADVANC and LH simulated stock prices are compared with these empirical prices. The ADVANC and LH empirical prices can be obtained from http://www.set.or.th/th/index.html. The data of ADVANC and LH empirical prices from July 9, 2010, to July 8, 2013, are used in the paper. These data are divided into two joint sets for two purposes. The first set (July 9, 2010–July 5, 2013) is used to estimate the drift rate and volatility. The second set (July 8, 2011–July 8, 2013) is used for model validations.The rate of return and volatility contained in the BMAP and the FBMAP are adaptive parameters based on time. Therefore, these parameters are not constant. In this section, the rate of return and volatility of ADVANC and LH stock prices are estimated.The ADVANC and LH closed prices from July 9, 2010, to July 5, 2013, are used to estimateμ
j and σ
j by using (8) and (9), respectively [7]. Consider the following:
(8)
μ
j
=
252
M
∑
i
=
1
M
R
i
,
(9)
σ
j
=
252
M
-
1
∑
i
=
1
M
(
R
i
-
R
-
)
2
,
where R
i is the return of stock price which can be computed by R
i
=
(
S
i
+
1
-
S
i
)
/
S
i
,
R
- is the average of return R
i, and M is the number of returns. The parameters μ
j and σ
j are estimated by the set of data as shown in Figure 1. In this figure, the data from July 9, 2010, to July 8, 2011, are used to estimate the initial μ
0 and σ
0. The data from July 9, 2010, to July 11, 2011, are used to estimate μ
1 and σ
1. The stock market is closed on the weekend. Therefore, the closed prices on July 9, 2011, and July 10, 2011, are not available. The data from July 9, 2010, to July 12, 2011, are used to estimate μ
2 and σ
2, and so on. Using the same procedure, the data from July 9, 2010, to July 5, 2013, are used to estimate μ
482 and σ
482.Figure 1
The flowchart to estimate the parametersμ
j and σ
j.From the estimation results, the parametersμ
j and σ
j of ADVANC and LH closed prices are shown in Figures 2 and 3, respectively.The historical drift rateμ
j for closed prices prediction from July 11, 2011, to July 8, 2013.
(a)
ADVANC
(b)
LHThe historical volatilityσ
j for closed prices prediction from July 11, 2011, to July 8, 2013.
(a)
ADVANC
(b)
LH


## 4. Stock Prices Mathematical Models

### 4.1. Brownian Motion Model with Adaptive Parameters (BMAP)

The BMAP can be considered by the SDE as shown in (1). The rate of return μ and the volatility σ are adaptive parameters and can be estimated using the flowchart in Figure 1. In the paper, the Euler discretization method is applied to solve the SDE. The solution of discretized form of the SDE (1) is denoted by S
j. Therefore, the Euler discretization form of (1) can be written in
(10)
S
j
+
1
=
S
j
+
μ
j
S
j
Δ
t
+
σ
j
S
j
Δ
W
j
,
where j is time index (
j
=
0
,
…
,
N
), N is the number of datasets, Δ
t is a sampling time, and μ
j and σ
j are estimated in the previous section. For the paper, N is equal to 484 and Δ
t is set to 1/252. The initial value S
0 is equal to stock price at July 8, 2011. The term Δ
W
j can be approximated by
(11)
Δ
W
j
=
Z
j
Δ
t
.
The random variable Z
j is the standard normally distributed random variable with mean = 0 and variance = 1. It is generated by method of Box and Muller [8].The ADVANC and LH stock prices calculated from the BMAP are simulated by MATLAB programming. These simulated data are compared with the second set data of empirical prices for a model validation. The average relative percentage error (ARPE) as given in (12) is the accuracy index in the paper. Consider the following:
(12)
ARPE
=
1
N
∑
i
=
1
N
|
X
i
-
Y
i
|
X
i
×
100
,
where N is the number of datasets, X
i is the empirical price (market price), and Y
i is the model price (simulated price).For the simulation results by the BMAP, Figure4 shows the empirical prices compared with the prices simulated by the BMAP for a given path of Brownian motion process. In the paper, the date period for simulation is between July 8, 2011, and July 8, 2013.The simulation results using the BMAP.
(a)
For ADVANC closed prices (ARPE = 23.268%)
(b)
For LH closed prices (ARPE = 18.483%)


### 4.2. Fractional Brownian Motion Model with Adaptive Parameters (FBMAP)

The FBMAP can be described by SDE (7). The rate of return μ and the volatility σ in this model are variable parameters depending on time. The μ and σ can be estimated using the block diagram as shown in Figure 1. The SDE (7) is solved by using the Euler discretization method. Therefore, the Euler discretization form of (7) can be written in
(13)
S
j
+
1
ε
=
S
j
ε
+
μ
j
S
j
ε
Δ
t
+
σ
j
S
j
ε
[
α
φ
j
Δ
t
+
Δ
W
j
ε
α
]
.
In (13), the S
j
ε is the discretized solution of the SDE (7). j, Δ
t, and N have the same meaning as the BMAP case. The parameter α
=
H
-
0.5, where H is Hurst index and H
∈
(
0,1
). The estimation of this parameter is shown in Section 4.2.1. The parameters μ
j and σ
j are calculated the same as those of the BMAP case. In the paper, N
,
Δ
t, and ε are set equal to 484, 0.005, and 1/252, respectively. The initial value S
0
ε is equal to stock price at July 8, 2011. The term Δ
W
j can be generated by (11). The term φ
j in (13) can be calculated by [9]
(14)
φ
j
=
j
Δ
t
N
∑
k
=
0
N
-
1
(
t
-
k
j
Δ
t
N
+
ε
)
α
-
1
Z
k
.
The random variable Z
k in (14) is the standard normally distributed random variable with mean = 0 and variance = 1. It is generated by Box and Muller method.The MATLAB programming is also used to calculate the ADVANC and LH closed prices in FBMAP. For model validation, these simulated data are compared with the empirical prices on July 8, 2011–July 8, 2013. The accuracy index in this case uses the average relative percentage error (ARPE) as calculated by (12).


#### 4.2.1. Parameter Estimation

The parametersα for ADVANC and LH stock prices are the unknown values. Therefore, the estimation of these parameters is proposed in this section. The parameter α is calculated by α
=
H
-
0.5. In this equation, H is Hurst index and 0
<
H
<
1. Hence -
0.5
<
α
<
0.5. Firstly, α is varied from -
0.5 to 0.5 with step size equal to 0.1. However, the parameter α cannot be equal to -
0.5 or 0.5. Therefore, -
0.49 and 0.49 are used instead of -
0.5 and 0.5, respectively. The 10,000 sample paths of Brownian motion process are considered. In each path, the average relative percentage error (ARPE) is computed using every value α. The simulation results for ADVANC and LH stock prices with α varied from -
0.49 to 0.49 are addressed in Tables 1 and 2, respectively.Table 1
ADVANC.
α
−0.49
−0.4
−0.3
−0.2
−0.1
0
0.1
0.2
0.3
0.4
0.49
Average of ARPE
12.383
14.713
21.128
32.555
42.932
21.496
114.55
3.8315
×
10
14
3.7474
×
10
210
Inf
Inf
SD of ARPE
3.5199
5.4755
9.3025
18.166
30.492
10.474
744.64
2.8687
×
10
16
Inf
Inf
InfTable 2
LH.
α
−0.49
−0.4
−0.3
−0.2
−0.1
0
0.1
0.2
0.3
0.4
0.49
Average of ARPE
13.362
17.722
28.059
44.745
59.309
29.370
118.19
5.0349
×
10
41
8.0212
×
10
268
Inf
Inf
SD of ARPE
3.3415
6.2432
13.703
31.977
59.243
15.935
749.36
4.0179
×
10
43
Inf
Inf
InfIt can be seen that the average and standard deviation of ARPE in case ofα
=
-
0.49 are the minimum in both stock prices (ADVANC and LH). Therefore, α
=
-
0.49 is chosen for ADVANC and LH cases.


#### 4.2.2. Model Validation

For the simulation results using the FBMAP, Figure5 shows the empirical prices compared with the prices simulated by the FBMAP (7) with the same scenario of Brownian motion of Figure 4. In Figure 5, the date period to simulate the stock prices using FBMAP is between July 8, 2011, and July 8, 2013.The simulation results using the FBMAP.
(a)
For ADVANC closed prices (ARPE = 6.9208%)
(b)
For LH closed prices (ARPE = 6.8788%)


### 4.3. Comparison of Accuracy Index between BMAP and FBMAP

For a given standard Brownian motion sample path, Figures4 and 5 show that the average relative percentage error (ARPE) of the FBMAP is smaller than those of the BMAP in case of ADVANC and LH.In general, the 10,000 scenarios or sample paths of Brownian motion process are considered. In each path, the ARPE is computed from both models. The comparison results between the BMAP and the FBMAP can be seen from Table3. It can be seen that the average ARPE of the FBMAP is less than the average ARPE of the BMAP. Moreover, the standard deviation of ARPE from the FBMAP is smaller compared with the BMAP. The simulation results show that the FBMAP can provide the small ARPE compared with the BMAP in case of ADVANC and LH.Table 3
The average and standard deviation of ARPE using the BMAP and FBMAP.
Stock name
Model
Average of ARPE
Standard deviation of ARPE
ADVANC
BMAP
21.496%
10.474%
FBMAP
12.383%
3.5199%
LH
BMAP
29.370%
15.935%
FBMAP
13.362%
3.3415%


## 4.1. Brownian Motion Model with Adaptive Parameters (BMAP)

The BMAP can be considered by the SDE as shown in (1). The rate of return μ and the volatility σ are adaptive parameters and can be estimated using the flowchart in Figure 1. In the paper, the Euler discretization method is applied to solve the SDE. The solution of discretized form of the SDE (1) is denoted by S
j. Therefore, the Euler discretization form of (1) can be written in
(10)
S
j
+
1
=
S
j
+
μ
j
S
j
Δ
t
+
σ
j
S
j
Δ
W
j
,
where j is time index (
j
=
0
,
…
,
N
), N is the number of datasets, Δ
t is a sampling time, and μ
j and σ
j are estimated in the previous section. For the paper, N is equal to 484 and Δ
t is set to 1/252. The initial value S
0 is equal to stock price at July 8, 2011. The term Δ
W
j can be approximated by
(11)
Δ
W
j
=
Z
j
Δ
t
.
The random variable Z
j is the standard normally distributed random variable with mean = 0 and variance = 1. It is generated by method of Box and Muller [8].The ADVANC and LH stock prices calculated from the BMAP are simulated by MATLAB programming. These simulated data are compared with the second set data of empirical prices for a model validation. The average relative percentage error (ARPE) as given in (12) is the accuracy index in the paper. Consider the following:
(12)
ARPE
=
1
N
∑
i
=
1
N
|
X
i
-
Y
i
|
X
i
×
100
,
where N is the number of datasets, X
i is the empirical price (market price), and Y
i is the model price (simulated price).For the simulation results by the BMAP, Figure4 shows the empirical prices compared with the prices simulated by the BMAP for a given path of Brownian motion process. In the paper, the date period for simulation is between July 8, 2011, and July 8, 2013.The simulation results using the BMAP.
(a)
For ADVANC closed prices (ARPE = 23.268%)
(b)
For LH closed prices (ARPE = 18.483%)


## 4.2. Fractional Brownian Motion Model with Adaptive Parameters (FBMAP)

The FBMAP can be described by SDE (7). The rate of return μ and the volatility σ in this model are variable parameters depending on time. The μ and σ can be estimated using the block diagram as shown in Figure 1. The SDE (7) is solved by using the Euler discretization method. Therefore, the Euler discretization form of (7) can be written in
(13)
S
j
+
1
ε
=
S
j
ε
+
μ
j
S
j
ε
Δ
t
+
σ
j
S
j
ε
[
α
φ
j
Δ
t
+
Δ
W
j
ε
α
]
.
In (13), the S
j
ε is the discretized solution of the SDE (7). j, Δ
t, and N have the same meaning as the BMAP case. The parameter α
=
H
-
0.5, where H is Hurst index and H
∈
(
0,1
). The estimation of this parameter is shown in Section 4.2.1. The parameters μ
j and σ
j are calculated the same as those of the BMAP case. In the paper, N
,
Δ
t, and ε are set equal to 484, 0.005, and 1/252, respectively. The initial value S
0
ε is equal to stock price at July 8, 2011. The term Δ
W
j can be generated by (11). The term φ
j in (13) can be calculated by [9]
(14)
φ
j
=
j
Δ
t
N
∑
k
=
0
N
-
1
(
t
-
k
j
Δ
t
N
+
ε
)
α
-
1
Z
k
.
The random variable Z
k in (14) is the standard normally distributed random variable with mean = 0 and variance = 1. It is generated by Box and Muller method.The MATLAB programming is also used to calculate the ADVANC and LH closed prices in FBMAP. For model validation, these simulated data are compared with the empirical prices on July 8, 2011–July 8, 2013. The accuracy index in this case uses the average relative percentage error (ARPE) as calculated by (12).


### 4.2.1. Parameter Estimation

The parametersα for ADVANC and LH stock prices are the unknown values. Therefore, the estimation of these parameters is proposed in this section. The parameter α is calculated by α
=
H
-
0.5. In this equation, H is Hurst index and 0
<
H
<
1. Hence -
0.5
<
α
<
0.5. Firstly, α is varied from -
0.5 to 0.5 with step size equal to 0.1. However, the parameter α cannot be equal to -
0.5 or 0.5. Therefore, -
0.49 and 0.49 are used instead of -
0.5 and 0.5, respectively. The 10,000 sample paths of Brownian motion process are considered. In each path, the average relative percentage error (ARPE) is computed using every value α. The simulation results for ADVANC and LH stock prices with α varied from -
0.49 to 0.49 are addressed in Tables 1 and 2, respectively.Table 1
ADVANC.
α
−0.49
−0.4
−0.3
−0.2
−0.1
0
0.1
0.2
0.3
0.4
0.49
Average of ARPE
12.383
14.713
21.128
32.555
42.932
21.496
114.55
3.8315
×
10
14
3.7474
×
10
210
Inf
Inf
SD of ARPE
3.5199
5.4755
9.3025
18.166
30.492
10.474
744.64
2.8687
×
10
16
Inf
Inf
InfTable 2
LH.
α
−0.49
−0.4
−0.3
−0.2
−0.1
0
0.1
0.2
0.3
0.4
0.49
Average of ARPE
13.362
17.722
28.059
44.745
59.309
29.370
118.19
5.0349
×
10
41
8.0212
×
10
268
Inf
Inf
SD of ARPE
3.3415
6.2432
13.703
31.977
59.243
15.935
749.36
4.0179
×
10
43
Inf
Inf
InfIt can be seen that the average and standard deviation of ARPE in case ofα
=
-
0.49 are the minimum in both stock prices (ADVANC and LH). Therefore, α
=
-
0.49 is chosen for ADVANC and LH cases.


### 4.2.2. Model Validation

For the simulation results using the FBMAP, Figure5 shows the empirical prices compared with the prices simulated by the FBMAP (7) with the same scenario of Brownian motion of Figure 4. In Figure 5, the date period to simulate the stock prices using FBMAP is between July 8, 2011, and July 8, 2013.The simulation results using the FBMAP.
(a)
For ADVANC closed prices (ARPE = 6.9208%)
(b)
For LH closed prices (ARPE = 6.8788%)


## 4.2.1. Parameter Estimation

The parametersα for ADVANC and LH stock prices are the unknown values. Therefore, the estimation of these parameters is proposed in this section. The parameter α is calculated by α
=
H
-
0.5. In this equation, H is Hurst index and 0
<
H
<
1. Hence -
0.5
<
α
<
0.5. Firstly, α is varied from -
0.5 to 0.5 with step size equal to 0.1. However, the parameter α cannot be equal to -
0.5 or 0.5. Therefore, -
0.49 and 0.49 are used instead of -
0.5 and 0.5, respectively. The 10,000 sample paths of Brownian motion process are considered. In each path, the average relative percentage error (ARPE) is computed using every value α. The simulation results for ADVANC and LH stock prices with α varied from -
0.49 to 0.49 are addressed in Tables 1 and 2, respectively.Table 1
ADVANC.
α
−0.49
−0.4
−0.3
−0.2
−0.1
0
0.1
0.2
0.3
0.4
0.49
Average of ARPE
12.383
14.713
21.128
32.555
42.932
21.496
114.55
3.8315
×
10
14
3.7474
×
10
210
Inf
Inf
SD of ARPE
3.5199
5.4755
9.3025
18.166
30.492
10.474
744.64
2.8687
×
10
16
Inf
Inf
InfTable 2
LH.
α
−0.49
−0.4
−0.3
−0.2
−0.1
0
0.1
0.2
0.3
0.4
0.49
Average of ARPE
13.362
17.722
28.059
44.745
59.309
29.370
118.19
5.0349
×
10
41
8.0212
×
10
268
Inf
Inf
SD of ARPE
3.3415
6.2432
13.703
31.977
59.243
15.935
749.36
4.0179
×
10
43
Inf
Inf
InfIt can be seen that the average and standard deviation of ARPE in case ofα
=
-
0.49 are the minimum in both stock prices (ADVANC and LH). Therefore, α
=
-
0.49 is chosen for ADVANC and LH cases.


## 4.2.2. Model Validation

For the simulation results using the FBMAP, Figure5 shows the empirical prices compared with the prices simulated by the FBMAP (7) with the same scenario of Brownian motion of Figure 4. In Figure 5, the date period to simulate the stock prices using FBMAP is between July 8, 2011, and July 8, 2013.The simulation results using the FBMAP.
(a)
For ADVANC closed prices (ARPE = 6.9208%)
(b)
For LH closed prices (ARPE = 6.8788%)


## 4.3. Comparison of Accuracy Index between BMAP and FBMAP

For a given standard Brownian motion sample path, Figures4 and 5 show that the average relative percentage error (ARPE) of the FBMAP is smaller than those of the BMAP in case of ADVANC and LH.In general, the 10,000 scenarios or sample paths of Brownian motion process are considered. In each path, the ARPE is computed from both models. The comparison results between the BMAP and the FBMAP can be seen from Table3. It can be seen that the average ARPE of the FBMAP is less than the average ARPE of the BMAP. Moreover, the standard deviation of ARPE from the FBMAP is smaller compared with the BMAP. The simulation results show that the FBMAP can provide the small ARPE compared with the BMAP in case of ADVANC and LH.Table 3
The average and standard deviation of ARPE using the BMAP and FBMAP.
Stock name
Model
Average of ARPE
Standard deviation of ARPE
ADVANC
BMAP
21.496%
10.474%
FBMAP
12.383%
3.5199%
LH
BMAP
29.370%
15.935%
FBMAP
13.362%
3.3415%


## 5. Conclusion

Two asset pricing models are presented in the paper. One is the Brownian motion model with adaptive parameters called BMAP and another one is the fractional Brownian motion model with adaptive parameters called FBMAP. The rate of return and volatility in both models are adaptive at any time. The driven process in the BMAP is Brownian motion, while the driven process in the FBMAP is a fractional Brownian motion. The BMAP and the FBMAP are applied to simulate the ADVANC and LH stock prices. The simulated prices from both models are compared with the empirical prices. The accuracy index ARPE is used in the paper. From the 10,000 scenarios of simulated prices of each model, the average and standard deviation of ARPE from both models show that the FBMAP provides a better appropriateness with the dataset than the BMAP in case of ADVANC and LH. Therefore, the FBMAP is suitable to predict the ADVANC and LH closed prices in the future.



---

*Source: 791418-2014-04-07.xml*