KRIGING SCENARIO FOR CAPITAL MARKETS

arXiv:cs/0509012v5 [cs.CE] 1 Aug 2009

T. SUSLO
Abstract. An introduction to numerical statistics.

1. Introduction
Let us consider an unknown constant mean m and a variance σ 2 of the random
variable V and the linear unbiased estimation statistics (weighted variable with
sum of its weights equal to one) at j ≥ n + 1
V̂j =

n
X

ωji Vi = ωji Vi

i=1

with minimized (the best statistics)
E{[(Vj − m) − (ωji Vi − m)]2 }

(1)

E{[ωji Vi − m]2 }

(2)

and minimized

our aim is to satisfy for finite n the asymptotic property of (1)
lim E{[(Vj − m) − (ωji Vi − m)]2 } = σ 2

n→∞

and the asymptotic property of (2)
lim E{[ωji Vi − m]2 } = 0

n→∞

equivalent to
lim ωji Vi = lim ωji vi = m

n→∞

n→∞

to replace (conditionally) weighted variable ωji Vi by weighted average ωji vi .
2. The numerical statistics
The minimized variance of the estimation statistics (2)
!
n
X

i
i
2
2
ωj ρij − µj = ∓σ 2 ωji ρij − µj
E{[ωj Vi − m] } = ∓σ

(3)

i=1

and (1)

E{[(Vj − m) − (ωji Vi − m)]2 } = σ 2 1 ± ωji ρij + µj



(4)

2

T. SUSLO

in terms of correlation function ρij is given by the kriging system of equations


 1 


ωj
ρ11 . . . ρ1n 1
ρ1j
 .. 
 ..
 .. 
.. 
..
..



 .

.
. 
.
 ·  .  =  . 

 ωjn 
 ρn1 . . . ρnn 1 
 ρnj  .
(5)
1 ...
1
0
1
µj
|
| {z }
{z
}
| {z }
(n+1)×(n+1)

(n+1)×1

(n+1)×1

To replace (conditionally) weighted variable by weighted average in (3)

E{[ωji Vi − m]2 } = min[ωji vi − m]2 = ∓σ 2 ωji ρij − µj

(6)

we have to find (on computer) a numerical disjunction of the best linear estimation
statistics which satisfies the asymptotic property of (4) constrained by numerical
approximation to root of the equation
ωji ρij + µj = 0

(7)

and fulfills the asymptotic property of (3) for j → ∞.
Proof. Since for j → ∞ the correlation has to vanish
 




ρ1j
ρ(|1 − j|)
1

 ..  
 . 
..
ξ → 0− (or ξ → 0+ )
 = ξ  .. 
 . =
.
ρnj
ρ(|n − j|)
1
| {z } |
| {z }
{z
}
n×1

n×1

n×1

the constraint (7) can be easy fulfilled

ωji ρij + µj = ξ (ωj1 + . . . + ωjn ) +µj = ξ + µj = 0
|
{z
}
1

the kriging system that does not

ρ11 . . . ρ1n
 ..
..
..
 .
.
.

 ρn1 . . . ρnn
1 ...
1
{z
|
(n+1)×(n+1)

or (if rewritten)

depend on j now
 1 

1
ω
 .. 
.. 


. 
 ·  . 
 ωn 
1 
−ξ
0
}
| {z }

=

(n+1)×1




ξ
 .. 
 . 


 ξ 
1
| {z }

(8)

(n+1)×1

Λω − ξF = ξF
and
F ′ω = 1 ,
where
ω

=




ω1
 .. 
 . , F
ωn
| {z }
n×1

=




1
 .. 
 . , Λ
1
| {z }
n×1



ρ11

=  ...
ρn1
|

...
..
.
...
{z

n×n


ρ1n
..  =
. 
ρnn
}

Λ′ ,

KRIGING SCENARIO FOR CAPITAL MARKETS

3

gives the classic least-squares weights for an estimate ω ′ v
−1 −1
Λ−1 F
ω = 2ξΛ−1 F = F ′ Λ−1 F
Λ F = ′ −1 ,
F Λ F
where from (6) constrained by (7) holds
min[ω ′ v − m]2 = ±σ 2 2µj = ∓σ 2 2ξ
based on

ξ → 0− (or ξ → 0+ )




v1


v =  ...  .
vn
| {z }
n×1

Example. The numerical least-squares estimator constrained by (7) which also
fulfills the asymptotic property of (3) is derived by frozen model (see Tab. 1).
Fig.
Index
n
j
ωji vi min[ωji vi − m]2 ωji ρij + µj
1
FTSE 100 132 426 8463
0.003150
-0.000058
2
DJI
113 356 13604
0.002569
0.000007
3
S&P 500 102 270 1789
0.033930
0.000252
i
Table 1. The numerical least-squares estimator ωj vi of final mean of
frozen model tries at some j for input sample of size n to fulfill the
asymptotic property of the best linear unbiased statistics.

FTSE 100
8000
7000
6000
5000
4000
3000
2000

100

200

300

400

Figure 1. FTSE 100 from 1 September 1989 up to 1 May 2009 (237
monthly close quotes). The classic (grey line) least-squares estimator of
mean compared for frozen model with the numerical (black dots) leastsquares estimator of mean at its co-ordinate. The assumption j → ∞ is
fulfilled at j = 426. Dashed vertical line denotes n = 132.

j

4

T. SUSLO

DJI
14000
12000
10000
8000
6000
4000
2000

50

100

150

200

250

300

350

j

Figure 2. DJI from 1 April 1991 up to 1 May 2009 (218 monthly close
quotes). The classic (grey line) and numerical (black dots) least-squares
estimator of mean for frozen model. The assumption j → ∞ is fulfilled
at j = 356. Dashed vertical line denotes n = 113.

Remark. Crisis emerges always when someone’s wish is to hide the cash in pocket
instead of spending it. This is also the nature of present crisis. The root of capitalist

S&P500
1800
1600
1400
1200
1000
800
600
400

50

100

150

200

250

Figure 3. S&P 500 from 3 April 1992 up to 1 May 2009 (206 monthly
close quotes). The classic (grey line) and numerical (black dots) leastsquares estimator of mean for frozen model. The assumption j → ∞ is
fulfilled at j = 270. Dashed vertical line denotes n = 102.

j

KRIGING SCENARIO FOR CAPITAL MARKETS

5

production consists in industrial production of high quality components and semimanufactured goods for the needs of world economy. Each dollar spent in any
place of the world will sooner or later reinforce the roots of capitalism. In situation
of the lack of such sound root, the effects of positive balance of payments can be
deposited only in foreign exchange reserves, to avoid immediate consumption. The
money that cannot be spent can be called dead money. Global financial crisis has
its cause in the outflow of money from the international financial network, towards
the monetary reserves of the systems which, having no sound capitalist roots, while
enjoying high share in global trade, froze the assets in reserves. Under the impact
of the crisis, the foreign exchange reserves held by above mentioned systems have
started to shrink. Sooner or later, they will feed the roots of capitalism. The
revitalized and released dead money will result in growth of inflation in industrial
countries. This development will herald the end of financial crisis and the beginning
of a new economic boom, since investors will start to buy up the shares in order to
gain protection against inflation.
References
[1] E. H. Isaaks and R. M. Srivastava, An Introduction to Applied Geostatistics, New York:
Oxford Univ. Press (1989).
[2] T. Suslo, Modern Statistics by Kriging, arXiv:cs/0609079 (2009).
[3] T. Suslo, Mean and Variance Estimation by Kriging, arXiv:cs/0409033 (2009).

