Blanchet All MC2QMC
Blanchet All MC2QMC
Blanchet All MC2QMC
1 Introduction
In the financial industry, there are three main approaches to investment: the
fundamental approach, where strategies are based on fundamental economic
principles, the technical analysis approach, where strategies are based on past
prices behavior, and the mathematical approach where strategies are based on
mathematical models and studies. The main advantage of technical analysis is
that it avoids model specification, and thus calibration problems, misspecification risks, etc. On the other hand, technical analysis techniques have limited
theoretical justifications, and therefore no one can assert that they are riskless, or even efficient (see [LMW00]).
Consider a nonstationary financial economy. It is impossible to specify and
calibrate models which can capture all the sources of instability during a long
time interval. Thus it might be useful to compare the performances obtained
by using erroneously calibrated mathematical models and the performances
obtained by technical analysis techniques.
To our knowledge, this question has not yet been investigated in the literature. The purpose of this paper is to present its mathematical complexity
and preliminary results.
Here we consider the case of an asset whose instantaneous expected rate
of return changes at an unknown random time. We compare the performances
of traders who respectively use:
a strategy which is optimal when the model is perfectly specified and
calibrated,
mathematical strategies for misspecified situations,
a technical analysis technique.
In all this paper, we limit ourselves to the case in which the traders utility function is logarithmic. Of course, it is a severe limitation from a financial point of view. This choice is also questionable from a numerical point
of view because logarithmic utilities tend to smoothen the effects of the different strategies. However, we will see that, even in this case and within a
simplified model, the analytical formulae are rather cumbersome and that our
analysis requires nonelementary mathematical and numerical tools. See also
the Remark 1 below.
Our study is divided into two parts: a mathematical part which, when possible, provides analytical formulae for portfolios managed by means of mathematical and technical analysis strategies; a numerical part which provides
quantitative comparisons between all these various strategies.
(3)
This model was considered by Shiryaev [Shi63] who studied the problem of
detecting the change time as early and reliably as possible when one only
observes the process (St )t0 .
Assumptions and Notation
The algebra generated by the observations at time t is denoted by
FtS := (Su , 0 u t) , t [0, T ].
Note that the Brownian motion (Bt )0tT is not adapted to the filtration
(FtS )t0 .
The Brownian motion (Bt )t0 and the random variable are independent.
The change time follows an exponential law 4 of parameter :
P ( > t) = et ,
t 0.
(4)
2
2
< r < 2
.
2
2
tn = nt.
Su du.
(7)
We suppose that, at time 0, the agent knows the history of the risky asset
prices before time 0 and has enough data to compute M0 .
At each tn , n [1 N ], the agent invests all his/her wealth into the risky
asset if the price Stn is larger than the moving average Mtn . Otherwise, he/she
invests all the wealth into the riskless asset. Consequently,
tn = 1(St M ) .
n
tn
(8)
Denote by x the initial wealth of the trader. The wealth at time tn+1 is
!
St0n+1
Stn+1
tn + 0 (1 tn ) ,
Wtn+1 = Wtn
Stn
Stn
and therefore, since St0n+1 /St0n = exp(rt),
WT = x
N
1
Y
n=0
tn exp(Rtn+1 Rtn ) exp(rt) + exp(rt) .
(9)
2 2 y
(1)
2
e
i2 /2
dzdy,
(10)
p =
2y
2 y
0
y
Z Z
3/2
(2 //2)2 (v) (1+z22 )
z2 2
2 2 y
(2)
)
(
2
2
p =
1
e
z1
2y
4
2
+z2
y2
0
R
y1
1 3/2 (1 //2)2 v (1+z2 )
z2
z1
z1
2 2 y1
2
1 i 2
i2 (v)/2
e
v/2
2 y2
2y1
2 y1
ev dy1 dz1 dy2 dz2 dv,
Z Z 1 3/2 (1 //2)2 (1+z2 )
z
z
2 2 y
(3)
2
e
i2 /2
dzdy,
p =
2y
2 y
0
y
Z
2
ze /4y z cosh(u)u2 /4y
iy (z) =
e
sinh(u) sin(u/2y)du.
y 0
(11)
(12)
(13)
Proof. The tedious calculation involves an explicit formula, due to Yor [Yor01],
Rt
for the density of ( 0 exp(2Bs )ds, Bt ). See [BSDG+ 05] for details.
3.4 Empirical Determination of a Good Windowing
One can optimize the choice of by using Proposition 1 and deterministic
numerical optimization procedures, or by means of Monte Carlo simulations.
In this subsection we present results obtained from Monte Carlo simulations,
which show that bad choices of may weaken the performance of the technical
analyst strategy. For each value of we have simulated 500,000 trajectories
of the asset price and computed the expectation E log(WT ) by a Monte Carlo
method. The parameters used to obtain Figure 1(a) and Figure 1(b) are all
equal but the volatility. It is clear from the figures that the optimal choice of
varies. When the volatility is 5 percent, the optimal choice of is around 0.3
4.8
4.88
4.87
4.79
E(log(W_T))
4.77
4.76
4.75
4.74
4.85
4.84
4.83
4.82
4.81
4.8
4.79
4.73
4.72
4.78
0
0.2
0.4
0.6
0.8
1.2
1.4
4.77
1.6
0.1
0.2
0.3
(b)
5.07
5.06
5.05
5.04
5.03
5.02
5.01
0.1
0.2
0.3
0.4
0.5
0.6
(a)
E(log(W_T))
E(log(W_T))
4.86
4.78
0.4
0.5
0.6
0.7
0.8
0.9
(c)
Parameter 1 2
Figure 1(a) 0.2 0.2 2
Figure 1(b) 0.2 0.2 2
Figure 1(c) 0.2 0.2 2
0.15
0.05
0.05
r
0.0
0.0
0.0
T
2.0
2.0
3.0
0.7
1=0.12=0.2
4.86
4.82
4.84
E[log(WT)]
4.84
4.8
4.82
4.78
4.8
4.76
4.78
4.74
0.5
1.5
4.76
0.5
1.5
1.5
1=0.12=0.1
1=0.22=0.1
4.7
4.7
4.69
E[log(WT)]
4.69
4.68
4.68
4.67
4.67
4.66
4.65
0.5
1
delta (years)
1.5
4.66
0.5
1
delta (years)
Fig. 2. Volatility and Optimal Moving Average Window Size: Plot of Expected
Value of the Log of Terminal Wealth vs. Window size with T = 2, = 2, - = 0.1,
- - = 0.15, and -. = 0.2
Remark 1. One can observe that the empirical optimal choices of are close
to the classical values used by the technical analysts, that is, around 200 days
or 50 days. One can also observe from Monte Carlo simulations that these
optimal values also hold when the traders utility function belongs to the
HARA family: see [BSDG+ 05].
sup E U (WT ).
A(x)
2
0
where F is the conditional a posteriori probability (5).
Let D the subset of the {FtS } progressively measurable processes :
[0, T ] R such that
Z
E
The bond price process S 0 () and the stock price S() satisfy
Rt
St0 () = 1 + 0 Su0 ()(r + (u))du,
Rt
St () = S0 + 0 Su () (1 + (2 1 )Fu + (u) + (u))du + dB u .
For each auxiliary unconstrained market driven by a process , the value
function is
V (, x) := sup Ex U (WT ()),
A(,x)
where
dWt () = Wt () (r + (t))dt + t (t)dt + (2 1 )Ft dt + (1 r)dt + dB t
Proposition 2. If there exists e such that
V (e
, x) = inf V (, x)
(14)
then there exists an optimal portfolio for which the optimal wealth is
Wt = Wt (e
).
(15)
:=
t
1 r + (2 1 )Ft + e(t)
Rt
+
rt
e (s)ds
0
Ht W t e
!
, (16)
Rt
es L1
s ds
,
Rt
t
s
1 + e Lt 0 e
L1
s ds
et Lt
e
+
Ht = exp 0
dB s
!
2
1 R t 1 r + e(s) (2 1 )Fs
0
+
ds ,
2
Here, v is the Lagrange multiplier which makes the expectation of the left hand
side equal to x for all x.
Proof. See Karatzas-Shreve [KS98, p. 275] to prove (15). We obtain (16) by
solving the classical unconstrained problem for e.
10
(17)
where
(1 r + (2 1 )Ft )
1 r + (2 1 )Ft
< 0,
if
2
1 r + (2 1 )Ft
e(t) =
(18)
[0, 1],
0
if
2
(1 r + (2 1 )Ft )
otherwise,
and, as above,
e (t) = inf (0, e(t)) .
Remark 2. The optimal strategies for the constrained problem are the projections on [0, 1] of the optimal strategies for the unconstrained problem.
Remark 3. In the case of the logarithmic utility function, when t is small and
thus before the change time with high probability, one has Ft close to 0; as,
by hypothesis, one also has 1r
0, the optimal strategy is close to 0 ; after
2
the change time , one has Ft close to 1, and the optimal strategy is close to
min(1, 2r
2 ). In both cases, we approximately recover the optimal strategies
of the constrained Merton problem with drift parameters equal to 1 or 2
respectively.
Using (18) one can obtain an explicit formula for the value function corresponding to the optimal strategy:
E log(WT ) = log(x) + rT
Z T Z "
2
a
1
+
1 r + (2 1 )
2 1 + r
1+a
2
0
0
a>
2 2 + r
#
2
1
a
1 r + (2 1 )
1
1 r
2 1 + r
2
1+a
<a<
r
2 2 + r
2
et (1 + a)g(a, t)dadt.
(19)
Here, g(a, t)da is the density function of
exp
11
2 1
(2 1 )2
s
+
s
ds,
Bs
2 2
1/2
(1 2s)e/s 2
s
ds =
(1 s)2+
with = 2 2 /(2 1 )2 .
1/2
(2s 1)e/s 2
s
ds
(1 s)2+
12
We thus are in a position to compute the wealth of the trader who uses
the strategy consisting in investing all of his/her money in the bond until K
and in the stock after K . The value of the portfolio at maturity T is
WT =
xS0K
ST 1(K T ) + xST0 1(K >T ) .
S K
In the particular case of the logarithmic utility function, one can exhibit an
exact formula for E(log(WT )). Unfortunately this new formula has a complexity similar to (19): see [BSDG+ 05]. However we can numerically compare
the performances of the two change time detection strategies (one based on
technical analysis, the other one based on a mathematical model), and the
optimal portfolio allocation strategy. Figure 3 illustrates, based on the typical
results that we have obtained so far, that the methods using mathematical
models have better performances than the technical analyst method.
4.85
Optimal allocation
Model and detect
Technical Analyst
E[log(Wt )]
4.8
4.75
4.7
4.65
4.6
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
Time
Fig. 3. Comparison
13
)R
)
+
2(
)(
)
t ,
t
2
1
2
1
1
1
2
2 2
2 2
R
1
t
et Lt 0 es Ls ds
Ft =
.
Rt
1
1 + et Lt 0 es Ls ds
6.2 On the Misspecified Optimal Allocation Strategy
Observing the stock price St , the trader computes a pseudo optimal allocation
by using the erroneous parameters 1 , 2 , and . Thus the value of his/her
misspecified optimal allocation strategy is
t = proj[0,1]
(1 r + (2 1 )F t )
,
2
Wt
rt
Z
= e exp
u d(eru Su )
Numerical Example
In this section, we compare numerically the performance of two traders who
respectively use a misspecified model and the true model. We fix the value
of 1 = 0.2, = 0.15, r = 0.0 and = 2.0, and we assume that they
are perfectly known by the trader. A contrario 2 is misspecified. Its true
value is 2 = 0.2. Figure 4 shows the functions t E(log(Wt )) for three
values of 2 . It suggests that it is better to overestimate 2 (2 > 2 ) than
to underestimate it (2 < 2 ).
6.3 On Misspecified Model and Detect Strategies
The erroneous stopping rule is
Z t
K
1
= inf t 0, et Lt
es Ls ds
0
p
1 p
1
where p is the unique solution in ( , 1) of the equation
2
Z 1/2
Z p
(1 2s)e/s 2
(2s 1)e/s 2
s
ds =
s
ds
(1 s)2+
(1 s)2+
0
1/2
14
2 = 0.2
2 = 0.3
2 = 0.1
4.85
E[log(Wt )]
4.8
4.75
4.7
4.65
4.6
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
Time
Fig. 4. Error on 2 for the optimal trader
with = 2 2 /(2 1 )2 .
The value of the corresponding portfolio is
0
W T = xS
K
ST
+ xST0 1(K >T ) .
1 K
SK ( T )
15
4.85
MSP Case II
Technical
Analyst
E[log(Wt )]
4.8
4.75
MSP Case I
4.7
4.65
4.6
0.2
0.4
0.6
0.8
Time
1.2
1.4
1.6
1.8
0.15
0.25
1 -0.3
2 (0,0.13)
0.25
1 -0.3
2 0.1
(0.2,)
1 -0.3
2 0.1
0.25
(0,1.5)
alyst may overperform the misspecified optimal allocation strategy but not
the misspecified model and detect strategy. One can also observe that, when
2 /1 decreases, the performances of well specified and misspecified model
and detect strategies decrease. See [BSDG+ 05].
16
Acknowledgment
This research is part of the Swiss national science foundation research program
NCCR Finrisk which has funded A. Diop during her postdoc studies at INRIA
and University of Z
urich.
References
[Ach00]
[BS02]
[Yor01]
17