Distribuciones continuas

Download Report

Transcript Distribuciones continuas

8. Distribuciones continuas
1
Transformaciones de variables aleatorias
Densidad
3 / 2 x 2  1  x  1
f ( x)  
 0 en el resto
x  1
0

3
F ( x )   1 / 2 ( x  1)
 1 x 1

Distribución
Transformación o cambio de
variable aleatoria
¿Cuál será la función de densidad de
probabilidad transformada g(y)?
1  x  1
Y  u( X )  2 X
y  u( x)  2 x
1
x  u ( y )  w ( y )  y2 / 2
G ( y )  P (Y  y )  P ( 2 X  y )  P ( X  y / 2 )  F ( y / 2 )
g ( y )  G ' ( y )  F ' ( y / 2)
1
 f ( y / 2)
2
1
2
 3 / 16 y  2  y  2
g ( y)  
 0 en el resto
2
y  2
 0

3
G ( y )  1 / 2 [( y / 2 )  1]

1 y2

2 y2
3
Probemos ahora con una transformación que no sea
biyectiva, como:
Y  u( X )  X
y  u( x)  x
2
2
1
x  u ( y)  w( y)  
G ( y )  P (Y  y )  P ( X
 F (
y )  F (
2
 y )  P (
1
y)
2
f (
y)
2
 f (
y
y  X 
y)
y)
g ( y )  G ' ( y )  F ' (
1
y
 F ' (
y
y)
1
2

y
1
y)
2
y
4
3 / 2 y 0  y  1
g ( y)  
 0 en el resto
y0
0

G ( y)  
y y
 1 y 1

0  y 1
5
Distribución log-normal Log-N(,)
Se trata de la densidad de probabilidad de una variable log x
distribuida según una función normal:
X  N ( , )
G ( y )  P (Y  y )  P ( e
X
g ( y )  G ' ( y )  F ' (log y )
 y )  P ( X  log y )  F (log y )
1
y
g ( y) 
Y e
X
 f (log y )
1
y
2

1
1
(log y   )
exp  
2
y
2

2 


 ;

y0
6
7
8
9
Distribución exponencial Exp ()
La distribución exponencial es el equivalente continuo de
la distribución geométrica discreta. Que recordemos era:
G ( p )  P ( X  x )  1  p  p ,
x
x  0 ,1 , 2 , ...
Describe procesos en los que nos interesa saber el
tiempo hasta que ocurre un determinado evento,
sabiendo que el tiempo que puede transcurrir desde
cualquier instante dado t, hasta que ello ocurra en un
instante tf, no depende del tiempo transcurrido
anteriormente.
10
Distribución exponencial Exp ()
Ejemplos de este tipo de
distribuciones son: el tiempo que
tarda una partícula radiactiva
en desintegrarse (datación de
fósiles o cualquier materia
orgánica mediante la técnica del
carbono 14) o el tiempo que
puede transcurrir en un servicio
de urgencias, para la llegada de
un paciente.
11
Distribución exponencial Exp ()
En un proceso de Poisson donde
se repite sucesivamente un
experimento a intervalos de
tiempo iguales, el tiempo que
transcurre entre la ocurrencia de
dos "sucesos raros" consecutivos
sigue un modelo probabilístico
exponencial. Por ejemplo, el
tiempo que transcurre entre que
sufrimos dos veces una herida
importante (o una coz de burro,
recuerda...)
12
Distribución exponencial Exp ()
f ( x)  e
2.0
 x
para x  0 ,   0




1.8
1.6
1.4
1.2
1.0



f ( x ) dx 
e
 x
0
 e
0.8


 x

1
0
0.6
0.4
0.2
0.0
0
Vida media
1
2
3
 
4


0
5
x e
6
x
dx 
7
8
1

13
dx
Distribución exponencial Exp ()

x
e
t
dt   e
0
t
x
 1 e
 x
0
1  e   x , x  0
F (x)  
x0
 0,
14
15
16
Relación entre la distribución de Poisson y la
exponencial
Entre las distribuciones de Poisson y Exponencial existen
importantes relaciones.
• Distribución de Poisson: Sea Y una P( ) que representa el
número de llegadas en un intervalo de tiempo fijo. Recuerda que
 es la esperanza de esta distribución.
• Distribución exponencial: En este mismo problema
consideramos ahora el tiempo que transcurre entre dos llegadas.
Sea X la v.a que representa dicho tiempo. Se puede demostrar
que entonces X se distrubuye como una Exponencial().
Propiedad de ‘falta de memoria’ de la
distribución exponencial
• Se dice que la distribución exponencial no tiene memoria.
Esto es
P ( X  s  t | X  t )  P ( X  s ).
para todo s, t  0.
• Interpretación: Supongamos que queremos determinar la
probabilidad de que llegue un cliente en la próxima media
hora. Esta propiedad nos dice que nos da igual conocer
cuando llegó el último cliente o calcular directamente cuál
es la prob. De que llegue en los prox. 30 min SIN tener en
cuenta el pasado.
El tiempo en que un producto está de moda en su mercado
se distribuye como una exponencial de parámetro 8 meses.
Si sabemos que ya lleva 5 de moda, ¿cuál es la
probabilidad de que dure 10 más?
• Sea X: tiempo que el producto está de moda. Nos
piden:
P ( X  10  5 | X  5)
• Por la propiedad de ausencia de memoria de las
distribuciones exponenciales sabemos que
P ( X  10  5 | X  5)  P ( X  10).
Por tanto
P ( X  10 )  1  P ( X  10 )  1  (1  e
 10 *8
) e
 10 *8
Tippex de Powerpoint
20
21
23
24
25
26
Fiabilidad
En instalaciones o
aparatos con posibilidad
de accidentes graves:
centrales nucleares,
aviones, coches,... es
imprescindible conocer la
probabilidad de que éstos
acontezcan durante la
vida del sistema.
28
Fiabilidad
Definimos la variable aleatoria:
T = tiempo durante el que el elemento funciona
satisfactoriamente antes de que se produzca un
fallo.
La probabilidad de que el elemento proporcione
unos resultados satisfactorios en el momento t se
puede definir como la fiabilidad o confiabilidad:
R(t) = P(T > t)
29
La infiabilidad Q(t) es la probabilidad de que ocurra un
fallo antes del instante t:
Q(t) = F(t) = 1 - R(t)
Sea λ(t) la tasa de fallos o averías por unidad de tiempo.
Supongamos que un elemento funciona en el instante t.
La probabilidad condicional de que se produzca una
avería entre el momento t y el t + dt puede escribirse:
P (t  T  t   t | T  t ) 
Q (t   t )  Q (t )
R (t )
1
R (t )  R (t   t )
R (t )
t

1
dR ( t )
R (t ) d (t )
  ( t );

R (t )  R (t   t )
  (t )  t
R (t )
  (t )
d  Ln R ( t ) 
   (t )
d (t )
R ( t )  Exp  


t
0
 ( t ) 30
dt 

La curva de la bañera
Curva típica de evolución de la tasa de fallos
Existencia inicial
de dispositivos
defectuosos o
instalados
indebidamente
con una tasa de
fallos superior a
la normal.
Esta tasa de fallos elevada
va disminuyendo
con el tiempo hasta alcanzar
un valor casi constante.
La tercera etapa
de fallos de
desgaste es
debida a la
superación de la
vida prevista del
componente
cuando
empiezan a
aparecer fallos
de degradación
como
consecuencia
del desgaste. Se
caracteriza por
un aumento
rápido de la tasa
de fallos.
Fallos normales o aleatorios. El comportamiento de la
tasa es constante durante esta etapa y los fallos son
debidos a las propias condiciones normales de trabajo
de los dispositivos o a solicitaciones ocasionales
31
superiores a las normales.
R ( t )  Exp  


t
0
 ( t ) dt 

Si la tasa de fallos o averías por unidad de tiempo es
constante: λ(t) = λ, tendremos que la fiabilidad es:
R ( t )  Exp  


t
0
 ( t ) dt   Exp (   t )

una densidad de probabilidad exponencial.
Esta fórmula de fiabilidad se aplica a todos los dispositivos
que han sufrido un rodaje apropiado que permita excluir los
fallos iniciales, y que no estén afectados aún por el
desgaste (la zona plana de la bañera).
32
R ( t )  Exp  


t
0
 ( t ) dt 

En 1951 Weibull propuso que la expresión empírica más
simple capaz de ajustar a una gran variedad de datos reales:

 tt
t
 t  t0 
0




(
t
)
dt


R
(
t
)

Exp





0

  


 tt
0
F ( t )  1  Exp   
  
f (t ) 
  t  t0 

  


 1







 tt
0
Exp   
  














r  

1

 


 y  t  33
t0

Distribución de Weibull W(r, )
X  Exp (  )
Y  X
G Y ( y )  P (Y  y )  P ( X
1/ r
 y )  P ( X  y )  FX ( y )
r
r
r 1
r 1
 y
g Y ( y )  G 'Y ( y )  F ' X ( y ) ry
g ( y )   ry
e
1/ r
 f X ( y ) ry
r
r
;
r
r 1
y0
34
Función generatriz de momentos
g (t )  E [e
Discreta
tX
]

e i P ( X  xi )
tx
i
g (t )  E [e
Continua
tX
]



tx
e f ( x ) dx
2
3


t
t
tX
2
3
g ( t )  E [ e ]    1  tx 
x 
x  ...  f ( x ) dx 

2!
3!



1  tm 1 
t
2
2!
m2 
t
3
3!
m 3  ...
k
mk 
d g (t )
dt
k
t0
35
Función característica
 (t )  E [ e
]
itX

e
itx i
P ( X  xi )
e
itx
f ( x ) dx
i
 (t )  E [ e
]
itX



Observemos que:
f ( x) 
1
2



e
 itx
 ( t ) dt
a partir de la anti-transformada de Fourier de la función
característica obtenemos la densidad de probabilidad.
36
Desarrollando en Taylor la función característica
alrededor de t = 0:
 ' ' (0)
 (t )   ( 0 )   ' ( 0 )t 
t  ... 
2

2!
 (t )  E [ e
itX
 ' ( t )  E [ iXe
(k )
(0)
t  ...
k
k!
]   (0)  1
itX
]   ' ( 0 )  E [ iX ]  im 1
 ' ' (t )  E [i X e
2
2
itX
]   ' ' (0 )  E [i X ]  i m 2
2
2
2
...
d  (t )
d  (0)
k
dt
k
k
 E [i X e
k
k
itX
]
 ( t )  1  im 1t 
dt
i
k
 E [i X ]  i m k
k
2
2!
m 2 t  ... 
2
k
i
k
k
k!
m k t  ...
k
1 d  (t )
k
mk 
i
k
dt
k 37
t0
Calculemos la esperanza y la varianza de la distribución exponencial.
 (t )  E [ e


e
itX
 (   it ) x
]
 ' (0) 
E[ X ] 
i
(   it )
i
e
itx
f ( x ) dx 

  it
e
2


2i

2i
x
dx 

  it
2
(   it )
3
2

 ' ' (0)
2
i


e e
itx
0
E[ X ] 
2
2

 (   it ) x
 ' ' (0) 
1


 ' ' (t ) 
i

 ' (0)

dx  
0
 ' (t ) 


2

2

2
 E [ X ]  ( E [ X ]) 
2
2
2

2

1

38
2

1

2
Sean {X1, X2, ... , XN } n variables aleatorias independientes
con funciones características {1(t), 2(t), ... , N(t) }, e
Y = X1+ X2+ ... + XN.
n
Entonces:
 Y (t ) 

i
(t )
i 1
Ejemplo:
Sean X1 = Exp(), X2 = Exp() ... , Xn = Exp() n variables
aleatorias independientes. ¿Cómo se distribuye Y = X1+
X2+ ... + Xn?
 k (t ) 

  it
n
 Y (t ) 

k 1
; k  1, 2 , ..., n

  


  it    it 
n
39
Distribución de Erlang Er(n,  )

f ( x) 
 (t )  E [ e

n
 (n) 

x
itX
n 1
x
 (n)

]
e
n


e
 (   it ) x
itx
n 1

n,   0 x  0
;



dx 
n
 (n) 

0
1
 ( n ) (   it )
n


0
u
n 1

f ( x ) dx 
0
n
e
x
itx
n
 (n)
 u 


   it 
x
n 1
n 1
e
u
e
x
dx 
1
  it
du 
  
e du 
 (n)  

n
 ( n ) (   it )
   it 
u

e

n
1
n
40
41
42
43
45
46
47
50