Página 68 - ANAlitica6

Versión de HTML Básico

Lilia Quituisaca-Samaniego, Juan Mayorga-Zambrano y Paúl Medina
Analíti a
k
6
Revista de Análisis Estadístico
Journal of Statistical Analysis
B Algoritmo general
Algorithm 1
:
(
T
c
,
T
s
) =
Ponzi
(
h
;
c
0
,
U
,
a
,
b
,
nSim
;
m
,
E
0
;
i
p
,
η
,
σ
1
,
σ
2
;
d
0
,
d
1
,
ω
)
Input
:
h
;
c
0
,
U
,
a
,
b
,
nSim
;
m
,
E
0
;
i
p
,
η
,
σ
1
,
σ
2
;
d
0
,
d
1
,
ω
Output
:
(
T
c
,
T
s
)
// Initialize
t
0
0;
p
0
c
0
;
C
0
c
0
;
1
ω
i
p
1
+
i
p
;
α
d
1
ω
d
0
ω
d
1
d
2
0
d
0
d
2
1
;
β
d
2
0
ω
d
1
ω
d
1
d
2
0
d
0
d
2
1
;
2
L
0
E
0
+
m
;
γ
E
0
m
;
λ
0
L
0
m
;
F
0
E
0
i
p
·
m
·
c
0
;
L
0
L
0
C
0
;
// Process
3
for
sMont
=
0, 1, ...,
nSim
do
4
S
0
1
1
U
;
I
0
1
U
;
S
(
t
0
)
S
0
;
I
(
t
0
)
I
0
;
5
k
1;
Z
L
0
;
paux
(
Math
.
abs
(
a
) +
Math
.
abs
(
b
))
/
periodo
;
haux
a
;
6
while
Z
>
0
do
7
t
k
k
·
h
;
8
haux
haux
+
paux
;
9
S
(
t
k
)
← −
a
·
haux
·
S
(
t
k
1
)
·
I
(
t
k
1
) +
S
(
t
k
1
)
;
10
I
(
t
k
)
haux
·
(
a
·
S
(
t
k
1
)
·
I
(
t
k
1
)
b
·
I
(
t
k
1
)) +
I
(
t
k
1
)
;
11
N
k
I
(
t
k
)
I
(
t
k
1
)
1;
N
k
Normal
N
k
, 1/4 ;
12
c
k
N
k
· C
k
1
;
C
k
← C
k
1
+
c
k
;
13
for
j
=
0, 1, ...,
k
do
14
if
k
j
d
1
then
ω
k
,
j
= (
k
j
)
·
[
α
·
(
k
j
) +
β
]
;
else
ω
k
,
j
=
ω
;
15
ω
k
,
j
Normal
(
ω
k
,
j
,
σ
2
1
)
;
16
if
j
<
k
then
p
k
,
j
p
k
1,
j
(
1
ω
k
,
j
)(
1
+
i
p
)
;
else
p
k
,
j
c
k
;
17
end
18
for
l
=
0, 1, ...,
k
1
do
19
η
l
Normal
(
¯
η
,
σ
2
2
)
20
end
21
P
k
m
·
k
j
=
0
p
k
,
j
;
E
k
E
0
·
k
1
l
=
0
(
1
+
η
l
)
;
22
ˆ
P
k
P
k
+
E
k
;
23
w
k
(
1
+
i
p
)
·
k
1
j
=
0
ω
k
,
j
p
k
1,
j
;
W
k
m
·
w
k
;
24
λ
k
(
1
+
η
k
1
)
·
λ
k
1
+
c
k
w
k
;
L
k
m
·
λ
k
;
25
F
k
F
k
1
P
k
·
i
p
+
η
k
1
·
L
k
1
;
26
for
j
=
0, 1, ...,
k
1
do
27
U
j
m
 
c
j
k
l
=
j
+
1
ω
l
,
j
p
l
,
j
 
;
28
end
29
U
k
m
·
c
k
;
R
k
k
j
=
0
U
j
;
30
L ←
L
k
C
k
;
E
k
λ
k
γ
w
k
;
31
Z
L
k
;
k
k
+
1;
32
end
33
end
34
66
Analítika,
Revista de análisis estadístico
, 3 (2013), Vol. 6(2): 51-66