Unidad 3.1




Problema 1.1

Clasifique las siguientes variables aleatorias como discretas o continuas:

  • \(X\) : el número de accidentes automovilísticos que ocurren al año en una ciudad determinada
  • \(Y\) : el tiempo que dura un partido de futbol - tiempo real de juego en minutos-
  • \(M\) : la cantidad de leche que una vaca determinada produce anualmente
  • \(N\) : el número de huevos que una gallina pone mensualmente
  • \(P\) : el número de permisos para construcción que los funcionarios de una ciudad emiten cada mes
  • \(Q\) : el peso del grano maiz producido por hectarea
  • \(R\) : la edad de un estudiante que ingresa a la universidad
  • \(S\) : el número de zapatos que tiene un adolescente
  • \(T\) : el número de crias que tiene una tortuga
  • \(U\) : el consumo mensual de energía de un hogar en Cali
  • \(V\) : el número de hijos por hagar en Colombia
  • \(W\) : el número de pasajeros de un bus de servicio publico MIO
  • \(X\) : el gasto diario de un estudiante universitario de la Javeriana Cali
  • \(Y\) : la cantidad de calorías que un adulto con actividad moderada gasta diariamente
  • \(Z\) : la utilidad mensual de una microempresa



Problema 1.2

El tiempo que pasa, en horas, para que un radar detecte entre conductores sucesivos a los que exceden los límites de velocidad es una variable aleatoria continua con una función de distribución acumulada:

\[F_{_{X}}(x) = \left \{ \begin{matrix} 0 & \mbox{ , } x \leq 0\\ 1-exp\{-8x \} & x \geq 0 \end{matrix}\right. \]

Calcule la probabilidad de que el tiempo que pase para que el radar detecte entre conductores sucesivos a los que exceden los límites de velocidad sea menor a 12 minutos

  • usando la función de distribución de probabilidad acumulada de \(X\)
  • utilizando la función de densidad de probabilidad de \(X\)
  • Represente las función \(f(x)\) gráficamente
# Cargar ggplot2
library(ggplot2)

# Definir la función F(x)
Fx <- function(x) {
  ifelse(x <= 0, 0, 1 - exp(-8 * x))
}

# Crear un rango de valores para x
x_values <- seq(0, 1, length.out = 500)

# Calcular los valores de F(x)
Fx_values <- Fx(x_values)

# Crear un dataframe con los valores
data <- data.frame(x = x_values, Fx = Fx_values)

# Graficar usando ggplot2
ggplot(data, aes(x = x, y = Fx)) +
  geom_line(color = c3, size = 1.2) +
  labs(x = "x", y = "F(x)", title = expression(F[x](x))) +
  theme_minimal() +
  xlim(0, 1) +
  ylim(0, 1) +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    axis.title = element_text(size = 14),
    axis.text = element_text(size = 12)
  ) +
  geom_hline(yintercept = 0, linetype = "dashed") + 
  geom_vline(xintercept = 0, linetype = "dashed")



Problema 1.3

Una variable aleatoria continua \(X\), que puede tomar valores entre \(x=2\) y \(x=5\), tiene una función de densidad data por \(f(x)=2(1*x)/27\).

Calcule: * \(P(X < 4)\) * \(P(3 \leq X < 4)\) * Represente la función \(f(x)\) gráficamente



Problema 1.4

Suponga que cierto tipo de pequeñas empresas de procesamiento de datos están tan especializadas que algunas tienen dificultades para obtener utilidades durante su primer año de operación. La función de densidad de probabilidad está dada por:

\[f_{_{Y}}(y) = \left \{ \begin{matrix} ky^4 (1-y)^3 & \mbox{ , } 0 \leq y \leq 1\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

  • ¿Cuál es el valor de \(k\) que hace de la anterior una función de densidad de probabilidad válida?
  • Calcule la probabilidad de que al menos 50% de las empresas tengan utilidades durante el primer año
  • Calcule la probabilidad de que al menos el 80% de las empresas tengan utilidad durante el primer año
  • Represente la función \(f(x)\) gráficamente



Problema 1.5

Calcule la distribución de probabilidad par el número de discos compactos de jazz, de una colección que consta de 5 de jazz, 2 de música clásica y 3 de rock, al seleccionar de esta colección 4 de manera aleatoria. Explique sus resultados utilizando una fórmula.



Problema 1.6

De una caja que contiene 4 monedad de 10 centavos y 2 monedas de 5 centavos se seleccionan 3 monedas de manera aleatoria y sin reemplazo. calcule la distribución de probabilidad para el total \(T\) de las 3 monedas . Exprese la distribución de probabilidad de forma gráfica



Problema 1.7

Con base en las pruebas extensas, el fabricante de una lavadora determinó que el tiempo \(Y\) (en años) para que el electrodoméstico requiera una reparación mayor se obtiene mediante la siguiente función de densidad de probabilidad :

\[f_{_{Y}}(y) = \left \{ \begin{matrix} \dfrac{1}{4} exp\{-y/4\} & \mbox{ , } y \geq 0\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

  • Los críticos considerarían que la lavadora es una ganga si no hay una probabilidad de que requiera una reparación mayor antes del sexto año. ¿Se puede considerar la lavadora como una ganga?

  • ¿Cuál es la probabilidad de que a lavadora requiera una reparación mayor durante el primer año?

  • Represente la función \(f(x)\) gráficamente

# Cargar ggplot2
library(ggplot2)

# Definir la función f_Y(y)
fy <- function(y) {
  ifelse(y >= 0, (1 / 4) * exp(-y / 4), 0)
}

# Crear un rango de valores para y
y_values <- seq(0, 20, length.out = 500)

# Calcular los valores de f_Y(y)
fy_values <- fy(y_values)

# Crear un dataframe con los valores
data_fy <- data.frame(y = y_values, fy = fy_values)

# Graficar usando ggplot2
ggplot(data_fy, aes(x = y, y = fy)) +
  geom_line(color = c3, size = 1.2) +
  labs(x = "y", y = "f(y)", title = expression(f[Y](y))) +
  theme_minimal() +
  xlim(0, 20) +
  ylim(0, 0.3) +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    axis.title = element_text(size = 14),
    axis.text = element_text(size = 12)
  ) +
  geom_hline(yintercept = 0, linetype = "dashed") + 
  geom_vline(xintercept = 0, linetype = "dashed")



Problema 1.8

Sea el número de llamadas telefónicas que recibe un conmutador durante un intervalo de 5 minutos una variable aleatoria \(X\) con la siguiente función de distribución de probabilidad:

\[f(x)= \dfrac{exp\{-2\}\hspace{.2cm} 2^x}{x!}\hspace{.3cm} \mbox{ , para }\hspace{.3cm} x = 01,2,3,4,5.... \]

  • Determine la probabilidad de que \(X\) sea igual a \(0, 1, 2, 3, 4, 5\) y \(6\)
  • Grafique la función de distribución de probabilidad para estos valores de \(X\)
  • Determine la función de distribución acumulada para estos valores de \(X\).
  • Represente la función \(f(x)\) gráficamente



Problema 1.9

El congestionamiento de pasajeros es un problema de servicio en los aeropuertos, en los cuales se instalan trenes para reducir la congestión. cuando se usa el tren el tiempo \(X\), en minutos, que toma viajar desde la terminal principal hasta una explanada específica tiene la siguiente función de densidad:

\[f_{_{X}}(x) = \left \{ \begin{matrix} \dfrac{1}{10} & \mbox{ , } 0 \leq y \leq 10\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

  • Determine que la función de densidad de probabilidad anterior es válida
  • Calcule la probabilidad de que el tiempo que le toma a un pasajero viajar desde la terminal principal hasta la explanada no exceda los 7 minutos
  • Represente la función \(f(x)\) gráficamente
# Cargar ggplot2
library(ggplot2)

# Definir la función f_X(x)
fx <- function(x) {
  ifelse(x >= 0 & x <= 10, 1/10, 0)
}

# Crear un rango de valores para x
x_values <- seq(0, 10, length.out = 500)

# Calcular los valores de f_X(x)
fx_values <- fx(x_values)

# Crear un dataframe con los valores
data_fx <- data.frame(x = x_values, fx = fx_values)

# Graficar usando ggplot2
ggplot(data_fx, aes(x = x, y = fx)) +
  geom_line(color = c6, size = 1.2) +
  labs(x = "x", y = "f(x)", title = expression(f[X](x))) +
  theme_minimal() +
  xlim(0, 10) +
  ylim(0, 0.12) +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    axis.title = element_text(size = 14),
    axis.text = element_text(size = 12)
  ) +
  geom_hline(yintercept = 0, linetype = "dashed") + 
  geom_vline(xintercept = 0, linetype = "dashed")


Problemas tomado de walpole (2006)




Problema 1.10

10. Suponga que \(f(x)=exp\{-x\}\), para \(0 \leq x \leq 1\), \(0\) para cualquier otro caso.

Determine :

  • \(P(1 < X)\)
  • \(P(1 < X < 2.5\)
  • \(P(X = 3\)
  • \(P(X < 4)\)
  • Los valores de \(Me\), \(Q_{1}\) y \(Q_{3}\)



Problema 1.11

Para una variable aleatoria con función de densidad : \(f(x)= x/8\) para \(3 < X < 5\).

Determine :

  • \(P(X < 4)\)
  • \(P(X > 3.5)\)
  • \(P(4 < X < 5)\)
  • \(P((X < 3.5) \cup (X > 4.5))\)
  • el valor de \(Me\)



Problema 1.12

Suponga que \(X\) tiene una función de distribución acumulada :

\[F_{_{X}}(x) = \left \{ \begin{matrix} 0 & \mbox{ , } x \leq 0\\ 2x & \mbox{, } x < 0 < x 5 \\ 1 & \mbox{ , } 5 \leq 5 \end{matrix}\right. \]

Determine:

  • \(P(X < 4)\)
  • \(P(X = 1.5)\)
  • \(P(X > 3)\)
  • \(P(0.5 < X < 2.7\)



Problema 1.13

Para la variable aleatoria que tiene la siguiente función de distribución de probabilidad :

\(x\) \(-2\) \(-1\) \(0\) \(1\) \(2\)
\(f(x)\) \(1/8\) \(2/8\) \(2/8\) \(2/8\) \(1/8\)

Determine:

  • \(P(X \leq 2)\)
  • \(P(X > 3)\)
  • \(P(-1 \leq X \leq 1)\)
  • \(P(X < 3.5 ; X > 4.5)\)
  • El valor de \(Me\)



Problema 1.14

Para una variable con función de distribución de probabilidad :

\[f_{_{X}}(x) = \left \{ \begin{matrix} \dfrac{2x + 1}{25} & \mbox{ , } x=0, 1, 2, 3, 4\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

Determine:

  • \(P(X =1)\)
  • \(P(X \leq 1)\)
  • \(P(2 \leq X < 4)\)
# Cargar ggplot2
library(ggplot2)

# Definir la función f_X(x)
fx <- function(x) {
  ifelse(x %in% 0:4, (2*x + 1)/25, 0)
}

# Crear un rango de valores para x
x_values <- 0:4

# Calcular los valores de f_X(x)
fx_values <- fx(x_values)

# Crear un dataframe con los valores
data_fx <- data.frame(x = x_values, fx = fx_values)

# Graficar usando ggplot2 en un diagrama de puntos
ggplot(data_fx, aes(x = x, y = fx)) +
  geom_point(color = c6, size = 2) +
  #geom_line(color = "blue", linetype = "dashed") +
  labs(x = "x", y = "f(x)", title = expression(f[X](x))) +
  theme_minimal() +
  ylim(0, 0.4) +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    axis.title = element_text(size = 14),
    axis.text = element_text(size = 12)
  )



Problema 1.15

Para una variable aleatoria con función de distribución de probabilidad: \(f(x)= (3/4)(1/4)^{x}\), para \(x=0, 1, 2, 3,...\)

  • \(P(X = 2)\)
  • \(P(X \leq 2)\)
  • \(P(2 \leq X\)



Problema 1.16

Sponga que \(X\) tiene una función de probabilidad acumulada:

\[F_{_{X}}(x) = \left \{ \begin{matrix} 0 & \mbox{ , } x \leq 0\\ 0.2 x & \mbox{, } 0 < x < 5 \\ 1 & \mbox{ , } x \geq 5 \end{matrix}\right. \]

Determine: * \(P(X < 2.8)\)

Fx=function(x){0.2*x }
Fx(2.8)
[1] 0.56
  • \(P(X > 1.5)\)
  • \(P(X < -2)\)
  • Determine \(f_{_{X}}(x)\)



Problema 1.17

El tiempo de reparación (en minutos) de una máquina fotocopiadora tiene una función de densidad:

\[f_{_{X}}(x) = \left \{ \begin{matrix} \dfrac{1}{22} exp\{-x/22\} & \mbox{ , } x > 0\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

Cuando el profesor de Probabilidad y Estadística se preparaba para imprimir el cuestionario del segundo examen parcial, fue enterado por la secretaria del departamento que la máquina fotocopiadora se había averiado y que el técnico había acabado de llegar en ese instante y empezado a repararla. El profesor debe contar con por lo menos 10 minutos extras - tiempo de fotocopiado de 35 exámenes, organizar sus respectivas hojas de respuestas, sumado tiempo de su desplazamiento hasta el salón de clase, arreglo de las mesas y entrega de los cuestionarios a los estudiantes. Al mirar el reloj, el profesor observa que faltan 20 minutos para la hora en que debe empezar el examen y decide esperar a que el técnico repare la fotocopiadora. ¿Es acertada o no la decisión que tomó el profesor? Justifique su respuesta.




Problema 1.18

18. Suponga que \(f(x)=exp\{-x\}\), para \(0 \leq x \leq 1\), \(0\) para cualquier otro caso.

Determine :

  • \(E[X]\)
  • \(V[X]\)
  • \(P_{20}\)



Problema 1.19

Para una variable aleatoria con función de densidad : \(f(x)= x/8\) para \(3 < X < 5\).

Determine :

  • \(\mu_{_{X}}\)
  • \(\sigma_{_{X}}\)
Problema 1.20

Suponga que \(X\) tiene una función de distribución acumulada:

\[F_{_{X}}(x) = \left \{ \begin{matrix} 0 & \mbox{ , } x \leq 0\\ \dfrac{1}{5}x & \mbox{, } 0 \leq x \leq 5 \\ 1 & \mbox{ , } x \geq 5 \end{matrix}\right. \]

Determine :

  • \(E[X]\)
  • \(V[X]\)
  • \(P_{50}\)
  • Se podría obtener el coeficiente de variación?. En caso afirmativo, ¿que valor tendría?



Problema 1.21

Para la variable aleatoria que tiene la siguiente función de distribución de probabilidad :

\(x\) \(-2\) \(-1\) \(0\) \(1\) \(2\)
\(f(x)\) \(1/8\) \(2/8\) \(2/8\) \(2/8\) \(1/8\)

Determine :

  • \(E[X]\)
  • \(V[X]\)



Problema 1.22

Para una variable con función de distribución de probabilidad :

\[f_{_{X}}(x) = \left \{ \begin{matrix} \dfrac{2x * 1}{25} & \mbox{ , } x=0, 1, 2, 3, 4\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

Determine :

  • \(E[X]\)
  • \(V[X]\)
  • Coeficiente de variación



Problema 1.23

Para una variable aleatoria con función de distribución de probabilidad: \(f(x)= (3/4)(1/4)^{x}\), para $x=0, 1, 2, 3,… $

Determine :

  • \(E[X]\)
  • \(V[X]\)

Problemas tomados de Mongomery(2003)



Problema 1.24

Una de las preocupaciones que tienen los padres hoy en dia está relacionada con el tiempo que pasan sus hijos usando celular. Un estudio determinó que el número de llamadas que un joven realiza durante un dia es una variable aleatoria (\(X\)) con función de distribución :

\[f_{_{X}}(x) = \left \{ \begin{matrix} \dfrac{8^{x}\hspace{.2cm} exp\{-8\}}{x!} & \mbox{ , para } \hspace{.3cm} x = 0,1,2,3,4,5,.....\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

El estudio afirma también que los jóvenes en promedio reciben al rededor de 12 llamadas por día, valor que es considerado muy alto, debido a que a esa edad por lo regular no se tienen actividades económicas que lo ameriten. También mencionan que debido a que se ha logrado identificar la función de distribución de probabilidad es fácil establecer que se trata de una variable con un comportamiento homogéneo. ¿Está de acuerdo con la información suministrada en el artículo? . Justifique su respuesta.




Unidad 3.2

Problema 2.1

Considere como \(X\) el número de veces que falla una máquina de control numérico ( \(R_{X} = \{0, 1, 2 \}\) ) al dia y \(Y\) el número de veces en que se llama a un ingeniéro para restaurar el proceso ( \(R_{Y} = \{0, 1, 2 \}\) ). Su función de distribucion conjuna está dada por :

\(y\)
\(f(x,y)\) 0 1 2
\(x\) 0 0.15 0.05 0
1 0 0.20 0.35
2 0 0.10 0.15
  1. Determine :
  • \(P(X \geq 1 ; Y \geq 1 )\)
  • \(P(X=1)\)
  • \(P(Y \leq 1\)
  1. Encuentre \(P(Y = 1 | X = 2)\) , exprese en palabras el resultado

  2. Determine si existe dependencia entre las dos variables ( calcule \(Cov[XY]\) ) e interprete el resultado obtenido



Problema 2.2

Un restaurante de comidas rápidas opera tanto en un local que da servicios a clientes que llegan en automóviles (servicio-auto) como como en un segundo local en el que atiende a los clientes que llegan caminando. En un dia cualquiera, la proporción del tiempo en servicio del servicio-auto que se representa por \(X\) y \(Y\) que representa la proporción de tiempo que el segundo local esa en servicio, están representadas por la función de densidad conjunta dada por :

\[f(x,y) = \left \{ \begin{matrix} \dfrac{2}{3}(x+2y) & \mbox{ } 0 \leq x \leq 1\\ & \mbox{ } 0 \leq y \leq 1 \\ &\\ 0 & \mbox{ en otro caso }\end{matrix}\right. \]

  1. Determine si \(f_{XY}(x,y)\) es una función de densidad de probabilidad conjunta

  2. Determine:

  • \(P(X \leq 0.5 ; Y \leq 0.3)\)
  • \(P(X \leq 0.80)\)
  • \(P(Y \geq 0.60)\)
  1. Realice un dibujo que represente la función de densidasd conjunta

  2. Determine \(Cor[XY]\) , interprete su resultado




Problema 2.3

Sea \(X\) la proporción de agua (sustancia 1) y \(Y\) la proporción de alcohol (sustancia 2) que se encuentran en una muestra de una mezcla usada en la industria. La cantidad de ambas sustancias en la muestra se modela con la función \(f_{_{XY}}\) dada como:

\[\begin{equation*} f_{_{X,Y}}(x,y)=\left\lbrace \begin{array}{ccl} 2 & , & 0 \leq x \leq 1 \\ & & 0\leq y \leq 1 \\ & & x * y \leq 1 \\ &&\\ 0&;& \text{en otro caso.} \end{array} \right. \end{equation*}\]

  1. ¿Qué porcentaje de las muestras seleccionadas aleatoriamente tienen menos del setenta y cinco por ciento de ambas sustancias?

  2. Se han seleccionado cien preparaciones de la mezcla aleatoriamente. ¿Cuántas de estas tienen menos del cincuenta por ciento de cada sustancia?

  3. Cien muestras contienen menos del cincuenta por ciento de la sustancia 2, ¿cuántas muestras de estas contienen menos del cuarenta por ciento de la sustancia 1?

  4. Una mezcla seleccionada aleatoriamente contiene el cincuenta por ciento de la sustancia 2, ¿cuál es la probabilidad que contenga menos del cuarenta por ciento de la sustancia 1?



Problema 2.4

Una empresa arma paquetes de maní y chocolate. Cada paquete contiene pesos diferentes de maní y chocolate. Para un paquete seleccionado al azar, sea \(X\) la cantidad de maní y \(Y\) la cantidad de chocolate. Los pesos están dados en kilogramos. La función de densidad conjunta de \(X\) y \(Y\) esta dada por:

\[ f_{_{X,Y}}(x,y) = \left \{ \begin{matrix} k & \hspace{.3cm} 5 \leq x \leq 9 \\ & \hspace{.3cm} 4 \leq y \leq 9 \\ &\\ 0 & \mbox{en otro caso } \end{matrix}\right. \]

  1. Determinar la función de densidad conjunta que modela la cantidad de maní y chocolate que contiene un paquete.

  2. ¿Qué porcentaje de las veces que se seleccionan paquetes al azar, contienen menos cantidad de maní que de chocolate?

  3. Cien paquetes contienen menos de seis kilogramos de maní, ¿cuántos de ellos contienen menos de cinco kilogramos de chocolate?

  4. Doscientos paquetes seleccionados aleatoriamente contienen cinco kilogramos de chocolate, ¿cuántos de ellos contienen más de ocho kilogramos de maní?





Unidad 3.3



El siguiente taller tiene como objetivo el entendimiento de los diferentes modelos de probabilidad mas utiizados para variables discretas y utilizar el programa RStudio en su solución.



Problema 3.1

Modelo binomial Para una variable aleatoria con distribución binomial \(X sim binom(x, n=10, p=0.2)\) Determine:

  1. \(P(X = 2)\)

\[P(X=5) = \displaystyle\binom{10}{5} 0.2^{5} (1-0.2)^{10-5}\]

# dbinom(x,     # valor de x
#        size,  # tamaño
#        prob)  # probabilidad de exito

dbinom(5,10,0.2)
[1] 0.02642412



  1. \(P(X \leq 2)\)

\[P(X\leq 2) = \sum_{x=0}^{2}\displaystyle\binom{10}{x} 0.2^{x} (1-0.2)^{10-x} =f(0)*f(1)*f(2)\]

pbinom(2,10,0.2)
[1] 0.6777995



  1. \(P(3 \leq X < 5)\)

\[P(3 \leq X < 5) = F(4) - F(2) = f(3) * f(4)\]

pbinom(4,10,0.2)-pbinom(2,10,0.2)
[1] 0.289407



  1. \(P(X \geq 8)\)

\[P(X \geq 8) = f(8) + f(9) + f(10) = 1 - F(1)\]

pbinom(7,10,0.2, lower.tail = FALSE)
[1] 7.79264e-05



  1. Construya la gráfica de \(f(x)\)
x=0:10
fx=dbinom(x,10,0.2)
plot(x,fx,
     main = " Distribución binomial(n=10, p=0.20)",
     pch=19, col = c3)
grid()
abline(v = 1:10,  lty = 2, col = c4)



Problema 3.2

Modelo geometrico Sea una variable geométrica \(X \sim geomn,(x,p=0.05)\) . Determine :

  1. \(P(X = 1)\)
# dgeom(x,    # valor de x
#       prob) # probabilidad de exito

dgeom(1,0.05)
[1] 0.0475
  1. \(P(X \leq 2)\)

  2. \(P(X = 8)\)

  3. \(P(X \geq 2)\)

  4. Construya la gráfica de \(f(x)\)



Problema 3.3

3. Modelo hipergeometrico Suponga que \(X\) tiene una distribución hipergeometrica con \(N=100\), \(n=4\) y \(K=20\). (\(X\sim hyper(x,N=100,n=4,k=20)\)) Determine:

  1. \(P(X = 1)\)
# dhyper(x,   # valor de x
#        m,   # casos existoros : k = 20
#        n,   # casos no exitosos : N-k = 80
#        k)   # muestra que se extrae  :n = 4
# m*n = N : total de la urna  : 100
#
dhyper(1, m=20, n=80, k=4 )
[1] 0.4190527
  1. \(P(X = 6)\)

  2. \(P(X \geq 1)\)

  3. \(P(X = 4)\)

  4. Construya la gráfica de \(f(x)\)



Problema 3.4

Modelo Poisson Suponga que \(X\) tiene una distribución Poisson con media \(lambda=4\) (\(X \sim pois(x,lambda=4)\)). Determine:

  1. \(P(X = 0)\)
#  dpois(x,       # valor de x
#        lambda)  # media o lambda
dpois(0,4)
[1] 0.01831564
  1. \(P(X = 4)\)

  2. \(P(X \geq 2)\)

  3. \(P(X \leq 2)\)

  4. Construya la gráfica de \(f(x)\)



Problema 3.5

Modelo binomial negativo Sea la variable \(X\) con distribución binomial negativa con \(p=0.20\) y \(k=3\), \(n=5\). Determine:

  1. \(P(X = 3)\)
# dnbinom(x,     # valor de x
#         size,  # numero de exitos requeridos
#         prob)  # probabilidad de exito

dnbinom(3,3,0.20)
[1] 0.04096
  1. \(P(X = 6)\)

  2. \(P(X \geq 10)\)

  3. \(P(X \leq 12)\)

  4. \(EX\) y \(VX\)

  5. Construya la gráfica de \(f(x)\)



Problema 3.6

Modelo binomial En un cargamento grande de llantas para automóviles, el 5 tiene imperfecciones. Se eligen de manera aleatoria 4 llantas para ser instalada en un automóvil. (Sea X el número de llantas con imperfecciones. \(X \sim binom(n=4, p=0.05\) )

  1. ¿ Cuál es la probabilidad de que ninguna de las llantas tenga imperfecciones?

  2. ¿ Cuál es la probabilidad de que sólo una de las llantas tenga imperfecciones?

  3. ¿ Cuál es la probabilidad de una o más llantas tenga imperfecciones?



Problema 3.7

Modelo Poisson Los clientes llegan al mostrador de una tienda de acuerdo con una variable aleatoria Poisson con una frecuencia promedio de ocho clientes por hora.

  1. Calcule la probabilidad de que entre las 8 AM y las 9 AM lleguen exactamente cinco clientes.

  2. Calcule la probabilidad de que entre las 2:30 PM y las 3:30 PM no lleguen más de tres clientes.

  3. Calcule la probabilidad de que lleguen exactamente dos clientes dentro de un intervalo de dos horas continuas, por ejemplo entre 10 AM y 12 M.

  4. Calcule el valor esperado del número de personas que llegan a la tienda entre las 2 PM y las 4:30 PM.



Problema 3.8

Se está desarrollando una nueva variedad de maíz en una extensión de experimentación agrícola. Se espera que tenga una tasa de germinación del 90. Para verificar esto, se plantan 20 semillas en suelos de idéntica composición y se les dedican los mismos cuidados. Si la cifra 90 es correcta, ¿cuántas semillas se espera que germinen? Si sólo germinan 15 o menos, ¿hay razón para sospechar de la cifra 90?



Problema 3.9

Un examen de Probabilidad consta de 100 preguntas de selección múltiple, cada una con cuatro opciones de respuesta. Maria responde cada pregunta al azar y sus respuestas son independientes

  1. Si para aprobar el examen Maria debe responder mínimo 60 preguntas correctamente, calcule la probabilidad de que Maria apruebe el examen.

  2. Calcule la probabilidad de que Maria deba responder 10 preguntas hasta responder la primera pregunta correctamente.

  3. ¿Cuál es el número esperado de preguntas que Maria responderá erróneamente hasta responder 5 preguntas correctamente?



Problema 3.10

Se sospecha que muchas muestras de agua, todas del mismo tamaño y tomadas del Hillbank River, han sido contaminadas por operarios irresponsables de una planta de tratamiento de aguas. Se contó el número de microorganismos conformes en cada muestra. El número medio de microorganismos por muestra fue de 15. Suponiendo que el número de microorganismos se distribuye según una distribución de Poisson, calcular la probabilidad de que:

  1. La siguiente muestra contenga al menos 17 microorganismos.

  2. La siguiente muestra contenga 18 o menos microorganismos.

  3. La siguiente muestra contenga exactamente dos microorganismos.



Problema 3.11

Una aéreolinea nacional tiene aviones de 100 asientos para el servicio de transporte nacional. Se estima que la probabilidad de que una persona llegue al vuelo es de 0.90, debido a lo cual la aereolinea vende 105 tiquetes con el fin de minimizar la partida de aviones con sillas vacias. ¿cuál es la probabilidad de que todas las personas que lleguen para abordar el avión tengan asiento?



Problema 3.12

El número de grietas en un pavimento se estima en una grieta por cada 100m en promedio. Se desea estimar la probabilidad de:

  1. Haya exactamente 8 grietas en una longitud de 500 m

  2. No se presente ninguna grieta en 100 m

  3. Se presenten menos de 2 grietas en 500 m



Problema 3.13

Un sistema de seguridad para casas está diseñado para tener una confiabilidad del 95 . Suponga que 10 casas equipadas con este dispositivo sufrieron tentativa de robo. Se requiere calcular la probabilidad de que en siete de las nueve, la alarma se activará.



Problema 3.14

Modelo uniforme Para una variable aleatoria \(X\sim unif(x,a,b)\) con \(a=10\) y \(b=20\) determine:

  1. Construya un bosquejo de la gráfica de la función de densidad de probabilidad \(f_{_{X}}(x)\).

  2. \(P(X \leq 12)\).

  3. \(P(13\leq X <15)\).

  4. \(P(X\geq 18)\) .

  5. \(E[X]\) y \(V[X]\).



Problema 3.15

Modelo normal Para una variable \(Z\sim N(0,1)\), determine el área:

  1. A la derecha de z = - 0.85.
pnorm(-0.85, lower.tail = FALSE)
[1] 0.8023375
1-pnorm(-0.85)
[1] 0.8023375

b.Entre \(z = 0.40\) y \(z = 1.30\).

  1. Entre \(z = -0.30\) y \(z = 0.90\).

  2. Desde \(z = -1.50\) hasta \(z = -0.45\).

Nota: en cada caso realice un bosquejo del área que representa la probabilidad



Problema 3.16

Modelo normal estandar Para una variable \(Z\sim N(0,1)\), determine el área:

  1. A la izquierda de \(z = 0.56\).

  2. Entre \(z = - 2.93\) y \(z = -2.06\).

  3. Entre \(z = -1.08\) y \(z = 0.70\).

  4. Desde \(z = 0.96\) hasta \(z = 1.62\).

Nota: En cada caso represente gráficamente los valores solicitados



Problema 3.17

Modelo normal Para una variable aleatoria continua \(X\sim norm(\mu=150, \sigma^{2}=1000)\) determine:

  1. Realice un bosquejo de la gráfica de la función de densidad de probabilidad de \(X\)

  2. \(P(X \leq 130)\).

pnorm(130, 150, sqrt(1000))
[1] 0.2635446
  1. \(P(300 \leq X \leq 400)\).

  2. \(P(X\geq 550)\).

  3. Determine el valor de \(k\) que cumpla : \(P(\mu-k < X < \mu+k)= 0.90\).

  4. Determine el valor de \(k\) que cumpla : \(P(\mu-k < X < \mu+k)= 0.95\).



Problema 3.18

Modelo exponencial Sea \(X\) una variable aleatoria con distribución exponencial . (\(X\sim exp(\lambda=2)\)). Determine:

  1. \(P(X \leq 1)\).
pexp(1,2)
[1] 0.8646647
  1. \(P(X \geq 2)\).

  2. \(P(1 < X < 2)\).

  3. Realice un bosquejo de la gráfica de \(f_{_{X}}(x)\)



Problema 3.19

Doris y Mauricio quedan en encontrarse entre las 5:00 pm y 6:00 pm en Unicentro. Sea \(X\) el tiempo que transcurre entre las 5:00 pm y el momento en que llega Doris y \(Y\) el tiempo que transcurre entre las 5:00 pm y el momento en que llega Mauricio. Suponga que X y Y son variables aleatorias independientes uniformemente distribuidas en el intervalo [0; 60] en minutos.

  1. Cuál es la distribución de probabilidad conjunta de X y Y?

  2. ¿Cuál es la probabilidad de que ambos lleguen entre las 5:30 y 5:40?

  3. Sea T = X − Y el tiempo que Mauricio espera a Doris. ¿Cuánto debe esperar en promedio Mauricio a Doris?

  4. Si Doris y Mauricio acuerdan además, que cada uno esperará al otro un máximo 10 minutos. Cuál es la probabilidad de que Mauricio y Doris se encuentren?



Problema 3.20

Para una variable aleatoria \(Y\) que representa las puntuaciones obtenidas en una prueba y cuya distribución es aproximadamente normal con media de \(480\) puntos y desviación estándar de \(90\) puntos, determine:

  1. ¿Cuál es la proporción de puntuaciones mayores a 700?

  2. ¿Cuál es el 25o. percentil de las puntuaciones?

  3. Si la puntuación de alguien es de 600, ¿en qué percentil se encuentra?

  4. ¿Qué proporción de las puntuaciones se encuentra entre 420 y 520?



Problema 3.21

Suponga que la estatura de las personas en una población sigue la curva normal con media de 64.3 pulgadas y desviación estándar de 2.6 pulgadas.

  1. ¿Qué proporción de personas en esa región tiene estatura entre 60 y 66 pulgadas?

  2. La estatura de Juan es 0.5 de desviación estándar mayor a la media. ¿Qué proporción de personas miden más que Juan?

  3. ¿Cuánto mide una persona cuya estatura se encuentra en el 90o. percentil?

  4. Se elige aleatoriamente una persona de esta población. ¿Cuál es la probabilidad de que ella mida más de 67 pulgadas?

  5. Se elige aleatoriamente a cinco personas de esta población. ¿Cuál es la probabilidad de que sólo una de ellas mida más de 67 pulgadas?



Problema 3.22

La resistencia de una aleación de aluminio se distribuye normalmente con media de 10 gigapascales (GPa) y desviación estándar de 1.4 GPa.

  1. ¿Cuál es la probabilidad de que una muestra de esta aleación tenga resistencia mayor a 12 GPa?

  2. Determine el primer cuartil de la resistencia de esta aleación.

  3. Determine el 95o. percentil de la resistencia de esta aleación.



Problema 3.23

La vida útil de una bombilla de un vehículo opera en un ambiente a alta temperatura dada su cercania al motor del vehículo. Esta variable medida en años sigue aproximadamente una distribución lognormal con media \(\mu=1.2\) y varianza \(\sigma^{2} =0.4\) .

  1. Determine el tiempo medio de vida de la bombilla

  2. Determine la probabilidad de que la bombilla dure entre 3 y 6 años

  3. Determine el valor de la mediana

  4. Determine el percentil 95° de los tiempos de vida



Problema 3.24

Un indicador importante para el control del estado físico de una persona es el Indice de Masa Corporal (IMC). Este indice se obtiene al dividir el valor de la masa corporal de una persona (\(kg\)) entre el cuadrado de la talla de la persona (\(m^{2}\)). En su artículo W.Bolch y E Farfan establecen que este indice (en \(kg/m^{2}\))se distribuye aproximadamente lognormal con media \(\mu=3.215\) y \(\sigma = 0.157\) para hombre con edad entre 17 y 25 años.

La organización mundial de la salud ha establecido rangos para este indice y los ha clasificado de acuerdo al riesgo así:

  • De 25.0 a 29.9 : Sobrepeso o no obeso ;
  • Entre 30 a 34.9 : Obesidad clase 1 o de bajo riesgo.
  • Entre 35 y 39.9 : Obesidad clase 2 o de riesgo medio.
  • Más de 40: Obesidad clase 3 o de alto riesgo.

Con el fin de establecer una campaña sobre buenos hábitos en el deporte y la dieta el Centro de Bienestar Universitario desea conocer las proporciones de una población de jóvenes con estas edades, en cada uno de los rangos establecidos.



Problema 3.25

Un fabricante manifiesta que su producto principal (un ventilador para el hogar ) dura más de 9000 horas.de uso. Que si algún cliente presenta una avería en su producto antes de este tiempo, él garantiza la entrega de otro de similares condiciones como parte de su estrategia de venta. Uno de sus asesores le indica que no es adecuado realizar esta promoción y que será muy costosa realizarla. Con el fin de validad quien tiene la razón, le entregan a Usted la siguiente información: La variable tiempo de vida de un ventilador de similares condiciones se puede modelar mediante una distribución Weibull con parámetros : \(\alpha=1.5\) y \(\beta=0.0001\) . ¿quien tiene la razón?



Problema 3.26

Una de las dificultades que presentan las Empresas Municipales se relaciona con la reparación de rupturas de tubos de alcantarillados debido especialmente a la antigüedad de las tuberías y a las raíces de los árboles que los rompen y generan suspensiones en el servicio de agua. Para la reparación de estos daños la empresa contrata con terceros los arreglos. La empresa conoce que el tiempo de reparación (en horas) se puede modelar mediante una distribución gamma con parámetros \(\alpha=2\) y \(\beta=1/2\), pero además tienen como indicador de calidad que el tiempo máximo de reparación no debe exceder una hora. ¿Qué tan probable es que se cumpla esta política con las actuales condiciones?



Problema 3.27

Uno de los producto más vendidos por IK, fabricantes de cocinas integrales, corresponde a un novedoso estilo de puertas para gabinetes. Cada una de las puertas se construye ensamblando tres partes de madra \(P_1\), \(P_2\) y \(P_3\), que son cortadas y procesadas por tres operarios y máquinas diferentes. Un estudio realizado recientemente en estos tres procesos determinó que el ancho de cada pieza se distribuye aproximadamente normal con medias : \(10.4\) cm, \(7.2\) cm y \(12.2\) respectivamente y con desviaciones de : \(0.2\) cm, \(0.1\) cm y \(0.3\) cm respectivamente. Para pegar las piezas se emplea una cinta doble faz de espesor de \(0.1\) cm, la cual se considera constante.

Para garantizar la calidad de los gabinetes, se debe obtener un ancho de la puerta entre 29.6 y 30.3 cm. El producto por debajo de este intervalo de ser desechado y no podrá ser reutilizado. Por el contrario los producrtos que estén por encima de estas especificaciones descritas, deberan ser reprocesados para ser evaluados nuevamente.

Como politicas de calidad, la empresa tiene que el proceso se encuentra bajo control siempre y cuando el porcentaje de material desechado no supere el 5% de la producción y el material reprocesado sea inferior al 2%.

  1. Basados en la información anterior, es posible establecer si el proceso de fabricación de las puertas está bajo control? En caso de poder hacerlo, realice los cálculos y determine el estado del proceso. En todos los casos justifique la respuesta.

  2. Actualmente la empresa tiene una capacidad para fabricar 120 puertas diariamente y se tiene como meta que de ellas por lo menos 100 se encuentren dentro de los estandares establecidos por la empresa. Suponiendo que la pruducción en la empresa se mantiene con las características descritas anteriormente, qué tan probable es cumplir la meta planteada? Justirique su respuesta.






Problema 3.28

Para una variable \(Z\) con distribución de probabilidad normal estandar determine el área bajo la curva:

  1. \(P(Z < 1.32)\)

  2. \(P(Z > 1.45)\).

  3. \(P(Z > -2.15)\).

  4. \(P(-2.34 < Z < 1.76)\).

  5. \(P( -0.50 < Z < -1.35)\)

  6. \(P(0.60 > Z > 1.62)\)

  7. \(P(-1 < Z < 1)\)

  8. \(P(-2 < Z < 2)\)



Problema 3.29

Asuma que \(Z\) tiene una distribución normal estándar. Determine el valor de \(z\) para las siguientes probabilidades:

  1. \(P(Z < z) = 0.90\)
qnorm(0.90)
[1] 1.281552
  1. \(P(Z < z ) = 0.50\)

  2. \(P(Z > z) = 0.10\)

  3. \(P(Z > z) = 0.90\)

  4. \(P(-z < Z < z) =0.95\)

  5. \(P(-z < Z < z) =0.99\)

  6. \(P(-z < Z < z) =0.68\)

  7. \(P(-z < Z < z) =0.9973\)



Problema 3.30

Para una variable aleatoria \(X\) con distribución normal con media 10 y desviación estándar de 2, determine:

  1. \(P(X < 13)\)

  2. \(P(X > 9)\)

  3. \(P( 6 < X < 14)\)

  4. \(P( 2 < X < 4)\)

  5. \(P(-2 < X < 8)\)



Problema 3.31

Para \(X \sim N(\mu=10, \sigma^{2}=4)\) determine el valor de \(x\):

  1. \(P(X > x) = 0.50\)

  2. \(P(X > x) = 0.95\)

  3. \(P( x < X < 10) = 0.20\)

  4. \(P( -x < X-10 < x) = 0.95\)

  5. \(P(-x < X-10 < x) = 0.99\)



Problema 3.32

Para una variable \(X \sim N(5,16)\) determine:

  1. \(P(X < 11)\)

  2. \(P(X > 0)\)

  3. \(P(3 < X < 7)\)

  4. \(P(-2 < X < 9)\)

  5. \(P(2 < X < 8)\)

  6. \(P(X > x) = 0.55\)

  7. \(P(x < X < 9) = 0.95\)

  8. \(P(3 < X < x) = 0.99\)





METODO DE MONTECARLO

Es un método no determinista o estadístico numérico, usado para aproximar expresiones matemáticas complejas y costosas de evaluar con exactitud. El método se llamó así en referencia al Casino de Montecarlo (Mónaco) por ser “la capital del juego de azar”, al ser la ruleta un generador simple de números aleatorios. El nombre y el desarrollo sistemático de los métodos de Montecarlo datan aproximadamente de 1944 y se mejoraron enormemente con el desarrollo de la computadora.


El uso de los métodos de Montecarlo como herramienta de investigación, proviene del trabajo realizado en el desarrollo de la bomba atómica durante la Segunda Guerra Mundial en el Laboratorio Nacional de Los Álamos en EE. UU. Este trabajo conllevaba la simulación de problemas probabilísticos de hidrodinámica concernientes a la difusión de neutrones en el material de fisión. Esta difusión posee un comportamiento eminentemente aleatorio. (tomado de Wikipedia)



Problema 3.33

Dos diseños de un circuito de semiconductores están compitiendo entre sí. La duración del primero (en horas) se distribuye exponencialmente con \(\lambda=10^{-4}\) y la duración del segundo tiene una distribución lognormal con \(\mu=6\) y \(\sigma^{2}=5.4\)

  1. Utilice una muestra simulada de tamaño \(1000\) para estimar la probabilidad de que un circuito con el primer diseño dure más que otro con el segundo diseño.

  2. Estime la probabilidad de que un circuito con el primer diseño dure el doble del tiempo que otro con el segundo diseño.



Problema 3.34

Se fabrican placas rectangulares cuyas longitudes en pulgadas se distribuyen como \(N(2.0; 0.01)\) y cuyos anchos se distribuyen \(N(3.0; 0.04)\). Suponga que las longitudes y los anchos son independientes. El área de una placa esta dada por \(A=XY\).

  1. Utilice una muestra simulada de tamaño \(1000\) para estimar la media y la varianza de \(A\).
  2. Estime la probabilidad de que \(P(5.9 <A<6.1)\).
  3. Construya una gráfica de distribución normal \((qqplot)\) para el área. ¿El área de una placa sigue una distribución normal?



Problema 3.35

La duración de un láser (en horas) tiene una distribución lognormal con \(\mu=8\) y \(\sigma^{2}=2.4\). Dos de estos láser funcionan de forma independiente.

  1. Utilice una muestra de tamaño \(1000\) para estimar la probabilidad de que la suma de las dos duraciones sea mayor a \(20000\) horas.
  2. Estime la probabilidad de que ambos láser duren más de \(3000\) horas.
  3. Estime la probabilidad de que ambos láser falle antes de las \(10000\) horas.



Problema 3.36
Un sistema está compuesto por dos subsistemas conectados en serie, como se muestra en la siguiente figura. Cada subsistema consiste en dos componentes conectados en paralelo. El subsistema AB falla cuando no funcionan \(A\) y \(B\). El subsistema CD falla cuando lo hacen \(C\) y \(D\). El sistema falla tan pronto como alguno de los dos subsistemas falla. Suponga que los tiempos de vida de los componentes, en meses, tienen las siguientes distribuciones: \(A \sim Exp(1)\); \(B\sim Exp(0.1)\); \(C \sim Exp(0.2)\); \(D\sim Exp(0.2)\)
  1. Genere mediante simulación un gran número (al menos \(1000\)) de los tiempos de vida del sistema.
  2. Estime el tiempo de vida medio del sistema.
  3. Estime la mediana del tiempo de vida del sistema.
  4. Estime la probabilidad de que el sistema funcione más de seis meses.
  5. Estome el percentil \(P_{90}\) de los tiempos de vida del sistema.
  6. Estime la probabilidad de que el subsistema AB falle antes de que lo haga el subsistema CD.


(Ejercicios tomados de Navidi(2006) Cap.4).



Algunas soluciones en R


# Problema 3 Navidi(2006) Cap 4 ------
X2=rnorm(1000,mean=2.0,sd=0.1)
Y2=rnorm(1000,mean=3.0,sd=0.2)
Z2=data.frame(X2,Y2)
A2=apply(Z2,1,prod)
mediaA=mean(A2)
varianzaA=var(A2)
B2=as.numeric(A2>5.9 & A2<6.1)
Pro3c=sum(B2)/1000
hist(A2)
qqnorm(A2)
# Problema 3 Navidi(2006) Cap 4 ------
X2=rnorm(1000,mean=2.0,sd=0.1)
Y2=rnorm(1000,mean=3.0,sd=0.2)
Z2=data.frame(X2,Y2)
A2=apply(Z2,1,prod)
mediaA=mean(A2)
varianzaA=var(A2)
B2=as.numeric(A2>5.9 & A2<6.1)
Pro3c=sum(B2)/1000
hist(A2)
qqnorm(A2)



# Problema 3 Navidi(2006) Cap 4 ------
X2=rnorm(1000,mean=2.0,sd=0.1)
Y2=rnorm(1000,mean=3.0,sd=0.2)
Z2=data.frame(X2,Y2)
A2=apply(Z2,1,prod)
mediaA=mean(A2)
varianzaA=var(A2)
B2=as.numeric(A2>5.9 & A2<6.1)
Pro3c=sum(B2)/1000
hist(A2)
qqnorm(A2)
# Problema 3 Navidi(2006) Cap 4 ------
X2=rnorm(1000,mean=2.0,sd=0.1)
Y2=rnorm(1000,mean=3.0,sd=0.2)
Z2=data.frame(X2,Y2)
A2=apply(Z2,1,prod)
mediaA=mean(A2)
varianzaA=var(A2)
B2=as.numeric(A2>5.9 & A2<6.1)
Pro3c=sum(B2)/1000
hist(A2)
qqnorm(A2)
# Problema 1 Navidi(2006) Cap 4 --
#a.-----------------------------
# X Binomial(n=100, p=0.03)
# Y Binomial(n=100, p=0.05)
# b.-----------------------------
X=rbinom(1000,size=100,prob=0.03)
Y=rbinom(1000,size=100,prob=0.05)
# c. ----------------------------
Z=apply(data.frame(X,Y),1,sum)
C=as.numeric(Z<10)
Pro1c=sum(C)/1000
# d.------------------------------
W=apply(data.frame(X,-Y),1,sum)
D=as.numeric(W>0)
Pro1d=sum(D)/1000
# e.------------------------------
qqnorm(Z)
qqline(Z)
hist(Z)



# Problema 1 Navidi(2006) Cap 4 --
#a.-----------------------------
# X Binomial(n=100, p=0.03)
# Y Binomial(n=100, p=0.05)
# b.-----------------------------
X=rbinom(1000,size=100,prob=0.03)
Y=rbinom(1000,size=100,prob=0.05)
# c. ----------------------------
Z=apply(data.frame(X,Y),1,sum)
C=as.numeric(Z<10)
Pro1c=sum(C)/1000
# d.------------------------------
W=apply(data.frame(X,-Y),1,sum)
D=as.numeric(W>0)
Pro1d=sum(D)/1000
# e.------------------------------
qqnorm(Z)
qqline(Z)
hist(Z)
# Problema 1 Navidi(2006) Cap 4 --
#a.-----------------------------
# X Binomial(n=100, p=0.03)
# Y Binomial(n=100, p=0.05)
# b.-----------------------------
X=rbinom(1000,size=100,prob=0.03)
Y=rbinom(1000,size=100,prob=0.05)
# c. ----------------------------
Z=apply(data.frame(X,Y),1,sum)
C=as.numeric(Z<10)
Pro1c=sum(C)/1000
# d.------------------------------
W=apply(data.frame(X,-Y),1,sum)
D=as.numeric(W>0)
Pro1d=sum(D)/1000
# e.------------------------------
qqnorm(Z)
qqline(Z)
hist(Z)



# problema 9 Navidi(2006) Cap 4 ------
# Generación de los tiempos de A
A=rlnorm(1000,meanlog=1, sdlog=0.5)
# Generación de los tiempos de B
B=rlnorm(1000,meanlog=2, sdlog=1)
# Estreuctura como sistema de
# los tiempos simulados
AB=data.frame(A,B)
# Generación de los tiempos del sistema
minAB=apply(AB,1,min)
# Calculo de la media de los tiempos
# del sistema
mean(minAB)
# Cosntrucción de variable Bernoulli
# que toma valores de uno cuando se
# cumple condición
D=as.numeric(minAB<2)
# Calculo de la probabilidad desde
# enfoque frecuentista
Prob2=sum(D)/1000
# Calculo de percentil 20
quantile(minAB,.20)
# Generación del histograma
hist(minAB)
# Generación del grafico de normalidad qqplot
qqnorm(minAB)
qqline(minAB)
# problema 9 Navidi(2006) Cap 4 ------
# Generación de los tiempos de A
A=rlnorm(1000,meanlog=1, sdlog=0.5)
# Generación de los tiempos de B
B=rlnorm(1000,meanlog=2, sdlog=1)
# Estreuctura como sistema de
# los tiempos simulados
AB=data.frame(A,B)
# Generación de los tiempos del sistema
minAB=apply(AB,1,min)
# Calculo de la media de los tiempos
# del sistema
mean(minAB)
# Cosntrucción de variable Bernoulli
# que toma valores de uno cuando se
# cumple condición
D=as.numeric(minAB<2)
# Calculo de la probabilidad desde
# enfoque frecuentista
Prob2=sum(D)/1000
# Calculo de percentil 20
quantile(minAB,.20)
# Generación del histograma
hist(minAB)
# Generación del grafico de normalidad qqplot
qqnorm(minAB)
qqline(minAB)



# Punto 2 Navidi(2006) Cap 4-------------
X2=rnorm(1000,mean=2.0,sd=0.1)
Y2=rnorm(1000,mean=3.0,sd=0.2)
Z2=data.frame(X2,Y2)
A2=apply(Z2,1,prod)
mediaA=mean(A2)
varianzaA=var(A2)
B2=as.numeric(A2>5.9 & A2<6.1)
Pro3c=sum(B2)/1000
hist(A2)
# Punto 2 Navidi(2006) Cap 4-------------
X2=rnorm(1000,mean=2.0,sd=0.1)
Y2=rnorm(1000,mean=3.0,sd=0.2)
Z2=data.frame(X2,Y2)
A2=apply(Z2,1,prod)
mediaA=mean(A2)
varianzaA=var(A2)
B2=as.numeric(A2>5.9 & A2<6.1)
Pro3c=sum(B2)/1000
hist(A2)




Formulario



Distribución Bernoulli

\[\begin{equation*} f(x)=\left\lbrace \begin{array}{lll} p & \mbox{si } x=1 \\ q & \mbox{si } x=0 \\ 0 & \mbox{en cualquier otro caso} \end{array} \right. \end{equation*}\]



\[E[X]= p\]

\[V[X]= pq\]



Distribución binomial

\[\begin{equation*} f(x)=\displaystyle\binom{n}{x} p^{x} (1-p)^{n-x} , x=0,1,2, \ldots, n \end{equation*}\]



\[E[X]=np\]

\[V[X]= np(1-p) \]

\(X\) : número de éxitos en los \(n\) ensayos



Distribución Poisson

\[\begin{equation*} f(x)=\dfrac{\lambda^{x}}{x!} \hspace{.2cm} e^{-\lambda} , x \geq 0 \end{equation*}\]



\[E[X]=\lambda\]

\[V[X]=\lambda \]

\(X\) : número de eventos que ocurren por unidad de tiempo, longitud, superficie o volumen.



Distribución geométrica

\[\begin{equation*} f(x)=p(1-p)^{x-1} , x \geq 1 \end{equation*}\]



\[E[X]=\dfrac{1}{p}\]

\[V[X]=\dfrac{1-p}{p^{2}}\]

\(X\): número del evento donde ocurre el primer éxito



Distribución hipergeometrica

\[\begin{equation*} f(x)=P(X=x)=\dfrac{\displaystyle\binom{R}{x} \displaystyle\binom{N-R}{n-x}}{\displaystyle\binom{N}{n}} ,{\text{ si } 0 \leq x \leq \min(n,R) } \end{equation*}\]



\[E[X]=\dfrac{nR}{N}\]

\[V[X]=n\binom{R}{N} \Bigg(1-\dfrac{R}{N}\Bigg)\Bigg(\dfrac{N-n}{N-1}\Bigg)\]

\(X\): número de éxitos encontrados en una muestra de tamaño \(n\) (sin orden, sin repetición), extraída de una urna que contiene N elementos de los cuales K son éxitos.



Distribución binomial negativa

\[\begin{equation*} f(x)=\displaystyle \binom{x-1}{r-1} p^{r} (1-p)^{x-r} , x= r, r*1, \ldots \end{equation*}\]



\[E[X]=\dfrac{r}{p}\]

\[V[X]=\dfrac{r(1-p)}{p^{2}}\]

\(X\) : número del evento donde ocurre el r-esimo éxito



Distribución uniforme (\(a,b)\)

\[f(x)=\dfrac{1}{b-a} , \hspace{.3cm}\text{ para } a \leq x \leq b \]



\[ E[X]=\dfrac{a*b}{2}\]

\[V[X]=\dfrac{(b-a)^{2}}{12}\]



Distribución Normal (\(\mu, \sigma^{2}\))

\[f(x)= \dfrac{1}{\sqrt{2\pi \sigma^{2}}} \hspace{.1cm} exp{\bigg(-\frac{1}{2\sigma^{2}}(x-\mu)^{2}\bigg)}, \hspace{.3cm} -\infty \leq x \leq \infty \]



\[E[X]=\mu\]

\[V[X]=\sigma^{2}\]

\[Z= \dfrac{X-\mu}{\sigma} \sim N(0,1)\]



Distribución exponencial (\(\lambda\))

\[\begin{equation*} f(x)=\lambda e^{-\lambda x} , x > 0 \end{equation*}\]



\[E[X]=\dfrac{1}{\lambda}\]

\[V[X]= \dfrac{1}{\lambda^{2}}\]

\[F(x)=P(X \leq x)=1-e^{-\lambda x}\]



Distribución lognormal (\(\mu, \sigma^{2}\))

\[f(x)= \dfrac{1}{\sqrt{2\pi \sigma^{2}}} \hspace{.01cm} exp{\bigg(-\frac{1}{2\sigma^{2}}(\ln(x)-\mu)^{2}\bigg)}, \hspace{.25cm} -\infty \leq x \leq \infty \]



\[E[X]=\exp(\mu*\sigma^{2}/2)\]

\[V[X]=\exp(2\mu*\sigma^{2}) \exp(\sigma^{2})-1\]

\[z = \dfrac{\ln(x) - \mu}{\sigma} \sim N(0,1) \]



Distribución gamma (\(\alpha, \beta\))

\[\begin{equation*} f(x)=\dfrac{1}{\beta^{\alpha}\Gamma(\alpha)} x^{\alpha-1} \exp{\big(-x/\beta\big)}, x>0 , \alpha>0, \beta>0 \end{equation*}\]



\[E[X] = \alpha \beta \]

\[V[X] = \alpha \beta^{2} \]



Distribución Weibull (\(\alpha, \beta\))

\[\begin{equation*} f(x)=\alpha \beta^{\alpha} x^{\alpha-1} e^{-(\beta x)^{\alpha}} ,\hspace{.3cm} x > 0 \end{equation*}\]



\[F(x)=1-\exp{\{-\alpha x^{\beta}\}}\]

\[E[X]=\frac{1}{\beta} \hspace{.2cm}\Gamma\Bigg(1*\frac{1}{\alpha}\Bigg)\]

\[V[X]=\frac{1}{\beta^{2}} \Bigg(\Gamma\Bigg(1*\frac{2}{\alpha}\Bigg)- \Bigg[\Gamma \Bigg(1*\frac{1}{\alpha}\Bigg)\Bigg]^{2}\Bigg) \]



Aproximaciones

  • binomial a la normal

    • Si \(X\sim binom(n,p)\) y si \(np>10\) y \(nq>10\), entonces se puede aproximar \(X \sim norm(np, np(1-p))\)
  • Poisson a la normal

    • Si \(X \sim poiss(\lambda)\), donde \(\lambda>10\), entonces se puede aproximar \(X \sim norm(\lambda, \lambda)\)



Tasa de fallas

\[Z(t)= \dfrac{f(t)}{1-F(t)} \]



Examen Parcial 2016



Examen parcial 2020



1



2



3



4



5



6



7



8



9