Capítulo 2 (3).doc

CAPÍTULO II
2. MODELOS ESTOCÁSTICOS A UTILIZARSE PARA
IMPUTACIÓN DE DATOS
2.1 Introducción
Este segundo capítulo aborda el tema de las técnicas y principios
científicos que permiten la generación de números aleatorios, los mismos
que son necesarios para la simulación de sistemas que se explican
estocásticamente. Para esto, la sección 2.2 trata acerca de la Distribución
Uniforme; la siguiente sección detalla los Métodos para la Generación de
números aleatorios; después se describe un Método de Generación de
Variables Aleatorias no Uniformes y por último se muestra las pruebas de
hipótesis sobre los parámetros o la distribución de una población.

26
Elaborado por: G. Cuenca
2.2 Distribución Uniforme
Una variable aleatoria X tiene una distribución uniforme continua con
parámetros  y  si y solo sí su función de densidad f , está dada por
 








X
de
resto
el
para
,
0
β
α,
X
,
α
β
1
x
f )
(
Los parámetros  y  de esta variable son constantes reales con 
  ,
véase Gráfico 2.1.
(2.1)
 
x

 
1
f(x)
GRÀFICO 2.1
Efectos de la Imputación en el análisis de datos multivariados
Densidad de la Distribución Uniforme

27
La media y la varianza de la distribución uniforme están dadas por:
μ
2
α
β
α)
2(β
α)
α)(β
(β
α)
2(β
α
β
α)
2(β
α
α)
2(β
β
α)
2(β
x
dx
α
β
1
x
xf(x)dx
X
E
2
2
2
2
β
α
2
β
α


































 




)
(
GRÁFICO 2.2
Efectos de la Imputación en análisis de datos multivariados
Media de la distribución uniforme
Por lo tanto,  se ubica en el punto medio entre  y  , véase Gráfico 2.2
(2.2)
 
2

 
)
(x
f
x

28
Además si  ,
β
α,
U
~
X
  
 
   
 
12
α
β
12
α
2α
β
12
3β
6α
3α
4α
4α
4β
β
2α
α
3
α
αβ
β
4
4
β
2α
α
3
α
αβ
β
2
β
α
3
α
αβ
β
Var(X)
α
β
3
α
αβ
β
α
β
3
α
β
α
β
1
3
α
3
β
α
β
1
3
x
α
β
1
dx
x
α
β
1
dx
α
β
1
x
f(x)dx
x
)
E(X
μ
)
E(X
μ)
E(X
Var(X)
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
3
3
3
3
β
α
3
β
α
2
2
2
2
2
2
2




























 
















 






















































Para comprobar que la probabilidad de que un valor ocurra en un intervalo,
solo depende de la longitud de dicho intervalo, efectuamos lo siguiente:
(2.3)

29
GRÁFICO 2.3
Efectos de la Imputación en análisis de datos multivariados
Números en el intervalo X Є )
,
( 

Sea X ~ )
,
( 

U , tal que:
 








X
de
resto
el
para
,
0
β
α,
X
,
α
β
1
x
f )
(
Además, el intervalo (a, b) está incluido en  
β
α, al igual que (c, d), esto es:
       
β
α,
d
c,
y
β
α,
b
a,
Si 

Supongamos además que:
x
a
b 

 )
( y que igualmente x
d
c 

 )
( , esto significa que:
 
   
 
   
































x
c
d
dx
d
c
X
x
a
b
dx
b
a
X
d
c
b
a
1
1
,
P
1
1
,
P
b
a c d β

← Δ → ← Δ →
)
(x
f

30
Por tanto  
   
 

 





x
d
c
X
P
b
a
X
P ,
,
Caso particular: Si 1

 y 10

 ; X  )
10
,
1
(
U
En este caso la densidad f es tal que:
 
 
   
 
    3
x
entonces
9
3
1
4
9
1
9
1
4
,
1
P
4
d
,
1
c
si
9
3
2
5
9
1
9
1
5
,
2
P
5
b
,
2
si
de
resto
0
10
,
1
9
1
)
(
4
1
5
2

























dx
X
dx
X
a
X
X
x
f
En este caso particular se puede comprobar que escogiendo números que se
encuentran en intervalos de igual longitud 3

x , la probabilidad de que se
efectúe una lectura en dicho intervalo es la misma, cuando la variable
aleatoria es uniforme.
En cambio si escogemos intervalos de longitud diferente a tres, la
probabilidad de que algo ocurra en dicho intervalo, no es
9
3
.
 
   
9
5
1
6
9
1
d
9
1
6
,
1
P
5
1
6
entonces
6
,
1
si
6
1











 x
X
x
i
h
(2.4)

31
La Distribución Acumulada de X  )
,
( 

U está dada por:
















β
X
si
1
β
X
α
si
α
β
α
x
α
X
si
0
F(x)
x)
P(X
Si X  )
1
,
0
(
U ;









1
si
1
)
1
,
0
(
si
0
si
0
X
X
x
X
F(x)
esto es   x
x
F  ,  
1
,
0

x lo cual significa que:
.
,
99
.
0
)
99
.
0
(
)
99
.
0
(
15
.
0
)
15
.
0
(
)
15
.
0
(
10
.
0
)
10
.
0
(
)
10
.
0
(
etc
X
P
F
X
P
F
X
P
F









Por definición la función generadora de momentos de una variable aleatoria
continua X es:
  dx
x
f
e
e
E
(t)
M Xt
X
X )
(
t






la variable independiente es t , y por lo general estamos interesados en los
valores de t en una vecindad de cero, por ejemplo |t |< h
Ahora se calcula la Función Generadora de Momentos de la Distribución
Uniforme
(2.5)
(2.6)

32
   
 
0
t
α
β
t
e
e
e
e
e
t
1
dx
e
dx
e
(t)
M
,
f(x)dx
e
)
E(e
(t)
M
tα
tβ
tα
tβ
β
α
tX
tX
tX
X
tX
tX
X





























,
t
1
*
1
*
1
1
1












Esto es, (t)
M X no está definida en t=0.
2.3 Prueba de Bondad de Ajuste Ji Cuadrada 2

La prueba de bondad de ajuste 2
 se aplica a situaciones en las que
queremos determinar si un conjunto de datos se puede considerar como
una muestra aleatoria de una población que tiene una distribución dada.
El contraste de hipótesis y el estadístico de prueba utilizados para éste
análisis, se presentan en el Cuadro 2.1
Cuadro 2.1
Efectos de la Imputación en el Análisis de Datos Multivariados
Contraste de Hipótesis de la Prueba de Bondad de Ajuste
H0: La distribución de la población donde se obtuvo la muestra es F0 (x)
vs.
H1: No es verdad H0
Estadístico de Prueba es : 


m
1
i i
2
i
i
e
)
e
f
(
que sigue una distribución
2
 y con (m-p-1) grados de libertad
(2.7)

33
Donde m es el número de términos en la suma y p es el número de
parámetros que se estiman en el modelo con base en los datos
muestrales.
Con el propósito de ilustrar esta prueba de hipótesis, considere si los
siguientes números aleatorios provienen de una distribución Poisson con
parámetro conocido 3


Tabla 2.1
Prueba de Bondad de Ajuste
i
Frecuencia
Observada
i
f
Probabilidad
Poisson
3


i
e
0 18 0.050 22.000
1 53 0.149 65.700
2 103 0.224 98.600
3 107 0.224 98.600
4 82 0.168 73.900
5 46 0.101 44.440
6 18 0.050 22.200
7 10 0.022 9.680
8 3 0.012 5.280
Cuadro 2.2
Prueba de Bondad de Ajuste
H0: La distribución de la población donde se obtuvo la muestra es Poisson 3


vs.
H1: No es verdad H0
Estadístico de Prueba es : 


m
1
i i
2
i
i
e
)
e
f
(
= 6.828
Valor p = 0,998

34
De acuerdo a los resultados obtenidos mediante la prueba de bondad de
ajuste, el valor p es 0.998, por lo tanto no existe evidencia estadística
suficiente para rechazar la hipótesis nula, es decir, la distribución de la
población donde se obtuvo la muestra es Poisson 3

 .
2.4 Prueba de Kolmogorov-Smirnov
La prueba de bondad de ajuste KS es una alternativa a la prueba 2
 que
permite comprobar si una muestra aleatoria proviene de una población
con una distribución dada, pero se prefiere el uso de la prueba KS en el
caso de distribuciones continuas ya que esta prueba trabaja
directamente sobre las observaciones y en cambio la prueba 2
 trabaja
sobre los datos agrupados.
Recordemos que dada una muestra aleatoria n
X
X
X ,...,
, 2
1 y la muestra
ordenada )
(
)
2
(
)
1
( ,...,
, n
X
X
X , la distribución empírica de la muestra es:













 
(n)
1)
(k
(k)
(1)
n
X
X
si
1
1
n
1,2,...,
k
si
;
X
X
X
si
n
k
X
X
si
0
(x)
F̂

35
La prueba KS consiste en verificar el contraste de hipótesis:
Con el propósito de ilustrar esta prueba de hipótesis, se tiene una
matriz de datos cuyas columnas son muestras tomadas de cinco
poblaciones todas ellas Normal, independientes e idénticamente
distribuidas, con parámetros  =0 y 2
 =1, 5
x
4
M

Χ , i= 1,2,3 y
j= 1,2,3,4,5
Tabla 2.2
Matriz de Datos de variables aleatorias independientes
con distribución Normal (0, 1)
Tamaño de muestra n=4
0.464 0.137 2.455 -0.323 -0.068
0.906 -0.513 -0.525 0.595 0.881
-0.482 1.678 -0.057 -1.229 -0.486
-1.787 -0.261 1.237 1.046 -0.508
Cuadro 2.3
Efectos de la Imputación en el Análisis de Datos Multivariados
Contraste de Hipótesis de la Prueba de Kolmogorov-Smirnov
H0: La distribución de la población donde se obtuvo la muestra es F0 (x)
vs.
H1: No es verdad F0 (x)
Estadístico de Prueba es : (x)
F
(x)
F
max 0
n 
ˆ
que sigue una distribución D y con (n, p) grados de libertad

36
Tabla 2.3
Prueba de Kolmogorov-Smirnov
Xn (x)
F20
ˆ F0(x) (x)
F
(x)
F
max 0
n 
ˆ
-1.787 1/20 0.037 0.013
-1.229 2/20 0.109 0.009
-0.525 3/20 0.299 0.149
-0.513 4/20 0.305 0.105
-0.508 5/20 0.305 0.056
-0.486 6/20 0.312 0.012
-0.482 7/20 0.316 0.034
-0.323 8/20 0.375 0.026
-0.261 9/20 0.397 0.053
-0.068 10/20 0.472 0.028
-0.057 11/20 0.476 0.074
0.137 12/20 0.556 0.044
0.464 13/20 0.677 0.027
0.595 14/20 0.726 0.026
0.881 15/20 0.811 0.061
0.906 16/20 0.819 0.014
1.046 17/20 0.853 0.003
1.237 18/20 0.893 0.008
1.678 19/20 0.954 0.004
2.455 20/20 0.993 0.006
De acuerdo a los resultados obtenidos mediante la prueba de
kolmogorov-Smirnov, el valor p es 0.766, por lo tanto no existe
evidencia estadística suficiente para rechazar la hipótesis nula, es
decir, la distribución de la población donde se obtuvo la muestra es
Normal (0, 1).
Cuadro 2.4
Prueba de Kolmogorov-Smirnov
H0: La distribución de la población donde se obtuvo la muestra es Normal(0,1)
vs.
H1: No es verdad H0
Estadístico de Prueba es : (x)
F
(x)
F
max 0
n 
ˆ = 0.149
Valor p = 0.766

37
2.5 Generación de Números Pseudo Aleatorios U(0,1)
Los números “pseudo aleatorios” son la base en la construcción de los
modelos de simulación donde hay presencia de variables estocásticas,
ya que estos permiten el funcionamiento de las abstracciones con los
que un fenómeno que no se puede construir físicamente, sea
numéricamente construìdo o recreado.
Existe un gran número de métodos que permiten la generación de
números aleatorios entre 0 y 1, la importancia del método a utilizar radica
en los números que genera, ya que estos números deben cumplir ciertas
características para que sean válidos. Dichas características son:
- Ser uniformemente distribuidos.
- Ser estocàsticamente independientes lo cual significa que si 1
X y 2
X
son dos variables aleatorias, 1
X y 2
X , son independientes si y sólo
si )
(
)
(
)
,
( 2
2
1
1
2
1
12 X
f
X
f
X
X
f  ; siendo 12
f la distribución conjunta de
1
X y 2
X y además 1
f y 2
f las marginales de 1
X y 2
X
respectivamente.
- Además es recomendable que los períodos del generador sean
“largos” es decir sin repetición dentro de una longitud determinada de
la sucesión de valores generados. [2]

38
2.5.1 Generadores Congruenciales Lineales
La generación de números pseudos aleatorios se realiza a través
de una “relación de recurrencia”, es decir para una sucesión
n
1
0 ,...,
, X
X
X , es una expresión que define a cada término n
X , en
función de uno o más de los términos que le preceden. Los valores
de los términos necesarios para empezar a calcular se llaman
condiciones iniciales. Se han propuesto varios esquemas como los
métodos congruenciales: congruencial mixto y congruencial
multiplicativo. [2]
2.5.1.1 Método Congruencial Mixto
El Método Congruencial Mixto genera una sucesión de números
pseudo aleatorios en la cual el sucesor 1

n
X del número pseudo
aleatorio n
X es determinado justo a partir de n
X . Particularmente
para el caso del generador congruencial mixto la relación de
recurrencia es la siguiente:
m
c
aX
X n
n mod
)
(
1 

 ,
Donde
0
X > 0: representa la semilla y es un valor que elige el
investigador;
a > 0 : se denomina multiplicador;
c > 0 : es una constante aditiva la que se denomina incremento;
(2.8)

39
m es el “módulo”, siendo; 0
X
m  , a
m  y además c
m 
Esta “relación de recurrencia” nos dice que 1

n
X es el residuo de
dividir c
aXn  para el módulo. Es decir que los valores posibles de
1

n
X son 1
,...,
3
,
2
,
1
,
0 
m , tal que, m representa el número posible
de valores diferentes que pueden ser generados. [2]
Para ilustrar la generación de números pseudoaleatorios por medio
de este método, suponga que se tiene un generador en el cual los
valores de sus parámetros son: 5

a , 7

c , 4
0 
X y 8

m .
Como se puede apreciar en la Tabla 2.4 el “período del generador”
es ocho, esto es la sucesión se repite una vez que se obtuvo el
octavo número generado

40
TABLA 2.4
Método Congruencial Mixto
Números pseudos aleatorios del generador 8
mod
)
7
5
( n
1
n 

 X
X
n Xn (5Xn+7)/8 Xn+1
Números
Uniformes
0 4 3+3/8 3 0.375
1 3 2+6/8 6 0.750
2 6 4+5/8 5 0.625
3 5 4+0/8 0 0.000
4 0 0+7/8 7 0.875
5 7 5+2/8 2 0.250
6 2 2+1/8 1 0.125
7 1 1+4/8 4 0.500
8 4 3+3/8 3 0.375
9 3 2+6/8 6 0.750
10 6 4+5/8 5 0.625
11 5 4+0/8 0 0.000
12 0 0+7/8 7 0.875
El procedimiento para obtener los números pseudo aleatorios se
realiza de la siguiente forma, donde la semilla es 4:
8
mod
)
7
5
(
1 

 n
n X
X
si 0

n
375
.
3
8
3
3
8
27
8
7
)
4
(
5
8
mod
)
7
5
( 0
1







 X
X
donde
8
3
es el residuo y al dividir 3 para 8, el resultado es el
número uniforme 0.375.

41
si 1

n
75
.
2
8
6
2
8
22
8
7
)
3
(
5
8
mod
)
7
5
( 1
2







 X
X
donde el residuo es
8
6
y al dividir 6 para 8, el resultado es el
si 2

n
62
.
4
8
5
4
8
37
8
7
)
6
(
5
8
mod
)
7
5
( 2
3







 X
X
donde el residuo es
8
5
número uniforme 0.625
y así sucesivamente se calculan los restantes cinco números
uniformes de la sucesión, los cuales son: 0.000, 0.875, 0.250,
0.125 y 0.500.
Al analizar este ejemplo se podría pensar que el período de
todo generador es siempre igual a m . Sin embargo, esto no es
verdad ya que el período del generador depende de los valores
asignados a los parámetros 0
,
, X
c
a y m , es decir, se requiere

42
seleccionar valores adecuados para estos parámetros con el fin
de que el generador tenga “período largo”.
Con el fin de ilustrar el caso que se presenta cuando el período
del generador es menor que m , suponga que se tiene un caso
en el cual los valores de los parámetros son: 7
0 

 c
X
a y
10

m . Para estos valores, la sucesión de números pseudo
aleatorios y uniformes son mostrados en la Tabla 2.5. Se puede
apreciar que el período del generador es cuatro, lo cual deja
claro que una selección inadecuada de los valores de los
parámetros del generador, puede conducirnos a obtener
períodos indeseables.
TABLA 2.5
Método Congruencial Mixto
Números pseudoaleatorios del generador   10
mod
7
7
1 

 n
n X
X
n Xn (7Xn+7)/10 Xn+1
Números
Uniformes
0 7 5+6/10 6 0.600
1 6 4+9/10 9 0.900
2 9 7+0/10 0 0.000
3 0 0+7/10 7 0.700
4 7 5+6/10 6 0.600
5 6 4+9/10 9 0.900
6 9 7+0/10 0 0.000

43
De la misma forma que el ejemplo anterior los números pseudo
aleatorios se obtienen de la siguiente manera, donde la semilla
es 7:
  10
mod
7
7
1 

 n
n X
X
si 0

n
600
.
5
10
6
5
10
56
10
7
)
7
(
7
10
mod
)
7
7
(







 0
1 X
X
donde
10
6
si 1

n
900
.
4
10
9
4
10
49
10
7
)
6
(
7
10
mod
)
7
7
(







 1
2 X
X
donde
10
9
si 2

n
000
.
7
10
0
7
10
70
10
7
)
9
(
7
10
mod
)
7
7
(







 2
3 X
X

44
donde
10
0
si 3

n
700
.
0
10
7
0
10
7
10
7
)
0
(
7
10
mod
)
7
7
(







 3
4 X
X
donde
10
7
Azarang y García expresan:
“Se advierte la necesidad de establecer algunas reglas que
pueden ser utilizadas en la selección de los valores de los
parámetros, para que el generador resultante tenga período
completo”. [1]
El valor apropiado del módulo m debe ser el número entero
más grande que la computadora acepte, el multiplicador a debe
ser un entero impar no divisible para 3 ó 5, la constante aditiva
c , puede ser cualquier constante y el valor de la semilla 0
X , es
irrelevante, para el generador congruencial mixto, es decir, el
valor de este parámetro resulta tener poca o ninguna influencia
sobre las propiedades estadísticas de las sucesiones.

45
2.5.1.2 Método Congruencial Multiplicativo
El Método Congruencial Multiplicativo al igual que el congruencial
mixto genera una sucesión de números pseudos aleatorios en la
cual el sucesor 1

n
X del número pseudo aleatorio n
X es
determinado justo a partir de n
X , de acuerdo a la siguiente
relación de recurrencia:
m
aX
X n
n mod
1 
 ,
Al igual que el generador anterior, en éste también se debe
seleccionar adecuadamente los valores de los parámetros 0
, X
a y
m , con el fin de asegurar un período largo para las sucesiones
generadas por este método.
Para ilustrar la obtención del período de un generador utilizando el
Método Congruencial Multiplicativo, suponga que se tiene un
generador con los siguientes parámetros: 5

a , 5
0 
X y m=32.
Estos valores se muestran en la Tabla 2.6.
(2.9)

46
TABLA 2.6
Efectos de la imputación en el análisis de datos multivariados
Método Congruencial Multiplicativo
Números pseudoaleatorios del generador 32
mod
5
1 n
n X
X 

n Xn 5Xn / 32 Xn+1
Números
Uniformes
0 5 0+25/32 25 0.781
1 25 3+29/32 29 0.906
2 29 4+17/32 17 0.531
3 17 2+21/32 21 0.656
4 21 3+9/32 9 0.281
5 9 1+13/32 13 0.406
6 13 2+1/32 1 0.031
7 1 0+5/32 5 0.156
8 5 0+25/32 25 0.781
9 25 3+29/32 29 0.906
10 29 4+17/32 17 0.531
11 17 2+21/32 21 0.656
Se puede apreciar en Tabla 2.6 que el período del generador es
ocho, esto es la sucesión se repite una vez que se obtuvo el octavo
número generado.
El procedimiento para obtener los números pseudo aleatorios se
realiza de la siguiente forma, donde la semilla es 5.
32
mod
5
1 n
n X
X 

si 0

n
781
.
0
32
25
0
32
25
32
)
5
(
5
32
mod





 0
1 5X
X

47
donde
32
25
si 1

n
906
.
3
32
29
3
32
125
32
)
25
(
5
32
mod





 1
2 5X
X
donde el residuo es
32
29
si 2

n
531
.
4
32
17
4
32
145
32
)
29
(
5
32
mod





 2
3 5X
X
donde el residuo es
32
17
y así sucesivamente se calculan los restantes cinco números de la
uniformes de la sucesión, los cuales son: 0.656, 0.281, 0.406,
0.031 y 0.156.
2.6 Métodos de Generación de Variables Aleatorias No Uniformes
Generalmente en las simulaciones de sistemas estocásticos existen una
o varias variables aleatorias interactuando, las cuales siguen

48
distribuciones diferentes a la distribución uniforme. Por consiguiente,
para simular este tipo de variables, es necesario contar con un
generador de números uniformes y una función que transforme estos
números uniformes en valores de la distribución de probabilidad
deseada. Para esto, se suele utilizar el método de la transformada
inversa. [1]
2.6.1 Método de la Transformada Inversa
El “Método de la Transformada Inversa” utiliza la distribución
acumulada F(x) de una variable aleatoria X que se va a simular.
Puesto que F(x) está definida en el intervalo (0,1), y que además
F(x)=x para  
1
,
0

x se puede generar un número aleatorio
uniforme y y tratar de determinar el valor de la variable aleatoria
para la cual su distribución acumulada es igual a y. Recordemos
que F es una función sobreyectiva e inyectiva y por tanto un
isomorfismo, además 0
)
(
lim
x



x
F y 1
)
(
lim
x



x
F .
Para convertir a un valor x, tomado de una distribución específica,
a partir de un valor uniforme, se deberá encontrar y en términos de
x, a partir de:
(y)
F
x
(y)
F
x
F
F
y
F(x)
1
1
1
-





))
(
(
(2.10)

49
Este método tiene la dificultad principal de que en algunas
ocasiones es difícil encontrar la transformada inversa. Sin
embargo, si esta función inversa ya ha sido establecida, generando
números aleatorios uniformes se podrán obtener valores de la
variable aleatoria que sigan la distribución de probabilidad
deseada.
2.6.1.1 Distribución exponencial
Utilizando el “Método de la Transformada Inversa” a continuación
se desarrolla un generador de variables aleatorias con distribución
exponencial. La función de densidad de probabilidad de una
variable aleatoria exponencial con parámetro  es:









0
si
0
0
si
e
1 /
x
x
f(x)
x 

Su función acumulada F es:










0
si
0
0
si
e
1 /
x
x
F(x)
x)
P(X
x 
Aplicando el método de la transformada inversa, se tiene que si
X es exponencial con parámetro  y y uniforme con parámetro
0 y 1:
y
x

  
/
e
1
y
x



1
e / 
(2.11)
(2.12)

50
y)
x



1
ln(
e
ln / 
y)
x 

 1
ln(

Para ilustrar este generador de variables aleatorias con distribución
exponencial, sea  ,
β
1,
E
~
X tal que 2

 y  ,
0,1
U
~
Y si por
ejemplo y=0.25.
0.575
x
0.25)
x



 1
ln(
2
Este es un valor tomado de una población exponencial.
2.6.2 Procedimientos Especiales
Existen algunas distribuciones como las distribuciones normal,
binomial, poisson, etc. cuya simulación a través del método de
transformada inversa resulta complicada. Para estas distribuciones
es posible utilizar algunas de sus propiedades para facilitar y
agilizar el proceso de generación de números aleatorios.
2.6.2.1 Variables que siguen una Distribución Normal
Puesto que no es posible expresar la distribución acumulada de la
distribución normal en forma explícita, por ende no se puede
utilizar el método de la transformada inversa.
2
2
1
2
1
)
(





 

 



x
e
x
f 


 x
(2.13)
(2.14)

51
Entonces si se desea generar números aleatorios que sigan una
Distribución Normal con parámetros conocidos  y 2
 , se puede
hacer uso del Teorema del Límite Central el cual establece que la
suma de n variables aleatorias independientes se aproxima a una
distribución Normal con media  y varianza 2
 , a medida que n
se aproxima al infinito.
Si por ejemplo, n
X
X
X ,....,
, 2
1 es una sucesión de n variables
aleatorias independientes, tal que i
X ~ U(0, 1) , i=1,2,…n
Para la distribución uniforme con parámetros 0

 y 1
β  se
conoce que:
2
1

 y
12
1
2


Aplicando el Teorema del Límite Central a los i
X tenemos que:
12
1
i
i
n
n
n
X
Z





Donde Z es un número aleatorio tomado de una población que
sigue una distribución Normal con parámetros 0

 y 1
2

 .
Para 12

n se obtienen buenas aproximaciones, entonces se tiene
que:
(2.15)

52
6
12
1
.
12
)
2
/
1
(
12
12
1
.
12
12 12
1
12
1
12
1











 





i
i
i
i
i
i
X
X
X
Z

Recuérdese que i
X ~ U(0, 1)
Si se desea obtener números aleatorios X que sigan una
distribución Normal con media  y varianza 2
 , se parte del
siguiente resultado, que se relaciona con la distribución normal
estándar esto es, X ~ )
,
(
N 2

 , entonces




X
Z es normal
estándar.
2.6.2.2 Variables que siguen una Distribución Poisson
Una variable aleatoria X tiene una distribución poisson si y sólo si
su distribución de probabilidades está dada por:
x!
e
λ
x)
P(X
λ
x 

 para x=0, 1, 2, …
Entonces la generación de números al azar que sigan una
distribución poisson, se lo puede hacer aplicando el método de la
transformada inversa.
0
;
1





 i
p
i
i)
P(X
p i
1
i

(2.16)
(2.17)
(2.18)

Capítulo 2 (3).doc

Más contenido relacionado

La actualidad más candente (19)

Similar a Capítulo 2 (3).doc (20)

Capítulo 2 (3).doc