1 Probabilidad

<En un colegio hay dos deportes que se practican futbol y baloncesto,un 60% futbol y un 40% baloncesto, el porcentaje de niños lesionados que practican futbol es de un 6% y el de baloncesto es de un 4%,escogemos un niño al azar calcula las siguientes probabilidades.>

l

1.Sabiendo que el niño esta lesionado, ¿cuál es la probabilidad de que practique fútbol?

\(P(JugarF_lesion)\)=\(P(Futbol)\)*\(P(F_lesioando)\)=0.036%

\(P(JugarB_lesion)\)=\(P(Baloncesto)\)*\(P(B_lesioando)\)=0.016%

\(P(Lesionado)\)=\(P(JugarF_lesion)\)+\(P(JugarB_lesion)\)=0.052%

\(P(Futbol|Lesionado)\)=\(P(Futbol\cup F_Lesionado)\)/\(P(Lesionado)\)=0.6923%

2.Sabiendo que el niño no practica fútbol, ¿cuál es la probabilidad de que esté lesionado?

\(P(!futbol\cap Lesionado)\)=\(P(Baloncesto \cap B_lesionado)\)=0.016% \(P(!futbol)\)=1-\(P(Futbol)\)=0.40%

\(P(Lesionado|!futbol)\)=\(P(!futbol\cap Lesionado)\)/\(P(!futbol)\)=0.04%

3.calcula la probabilidad de que juegue a futbol y este lesionado

\(P(Futbol \cap F_lesionado)\)=\(P(Futbol)\)*\(P(F_lesionado)\)=0.036%

4.calcula la probabilidad de que juegue a futbol y este lesionado o que juegue a baloncesto y no este lesionado

\(P(Futbol \cap F_lesionado)\)=0.036% \(P(Baloncesto\cap !B_lesionado)\)=0.384%

\(P(Futbol\cap F_lesioando | Baloncesto\cap !B_lesionado)\)=\(P(Futbol \cap F_lesionado)\)+\(P(Baloncesto\cap !B_lesionado)\)=0.42%

2 Variable aleatoria exponencial

Un fabricante de coches eléctricos produce baterías para un nuevo modelo de coches. La duración de las baterías sigue una distribución exponencial con una vida media de 4 horas. El fabricante produce 1500 baterías para la primera versión del coche que tiene una demanda esperada de 1400 unidades. La probabilidad de que una batería falle es del 5%, y se asume que el comportamiento de las baterías es independiente

  1. ¿Cuál es la distribución de probabilidad para la duración de las baterías?

La función de densidad de probabilidad para esta distribución es:

\[ f(x; \lambda) = \lambda e^{-\lambda x} \] Sabemos que la vida media es de 4 horas, el parámetro \(\lambda\) es \(\frac{1}{4}\).

Entonces la distribucion de probabilidad es \[ f(x) = \frac{1}{4} e^{-\frac{1}{4} x} \]

  1. Calcula la media de la variable aleatoria que representa la duración de las baterías

Como la media (\(\mu\)) de la distribución es igual a la inversa de \(\lambda\).

\[ \mu = \frac{1}{\lambda} = 4 \]

  1. Calcula la mediana de la distribución de duraciones de las baterías.

vamos a calcular la mediana de la distribución asi \(\mu \ln(2)\), que en este caso es \(4 \ln(2) \approx 2.773\) horas.

  1. Supongamos que se seleccionan al azar dos baterías independientes. ¿Cuál es la probabilidad de que ambas duren más de 30 horas?

\[ P(X > 30) = \int_{30}^{\infty} \frac{1}{4} e^{-\frac{1}{4} x} \,dx \]

lo que es aproximadamente 0.0183.

Vamos a calcular la probabilidad del conjunto multiplicando las probabilidades individuales

\[ P(X_1 > 30 \cap X_2 > 30) = P(X_1 > 30) \times P(X_2 > 30) \]

\[ P(X_1 > 30 \cap X_2 > 30) = 0.0183 \times 0.0183 \approx 0.000335 \]

***

3 Distribuciones de probabilidad

La empresa APPLE ofrece garantía contra que se caiga el móvil al suelo y se rompa. Según las estadísticas, la probabilidad de que a una persona se le caiga el móvil por este tipo de garantía es de una cada dos mil personas. La compañía estima que puede vender en promedio cinco mil garantías al año. Calcula:

1- Que el número de caídas cubiertos por la garantía no pase de cuatro por año?

Primero tenemos que calcular la tasa media de eventos por año:

r lambda <- 5000 * (1/2000)

Ahora ya podemos calcular el número de caídas:

\(P(X\leq 4)\)=

{r} ppois(4,lambda)

También se puede poner así:

\(P(X\leq 4)\)=ppois(2,1/3)=r x <-ppois(4,lambda)

2- Número de caídas esperados por año?

Ya lo hemos calculado anteriormente:

r lambda <- 5000 * (1/2000)

3- Que el número de caídas sea superior a tres por año?

\(1 - P(X\leq 3)\)=r x <- 1 - ppois(3, lambda)

4- Que ocurran nueve caídas por año?

\(P(X = 9)\)=r x <- dpois(9, lambda)

4 Probabilidad

La media de los pesos de 500 estudiantes de un colegio es de 70 kg y la desviación típica 3 kg.

Suponiendo que los pesos se distribuyen normalmente, hallar:

  1. La probabilidad de que pesen entre 60 kg y 75 kg.

\(P(60\leq X \leq75)\)={r} pnorm(75,70,3)-pnorm(60,70,3)

[1] 0.9517806
  1. La probabilidad de que pesen más de 90 kg.

\(P(X> 90)=P(Z>(90-70)/3)=P(Z>6.66)=1-P(Z\leq6.66)\)={r} 1-pnorm(6.66)

[1] 1.369138e-11
  1. La probabilidad de que pesen menos de 64 kg.

\(P(X < 64)\)={r} pnorm(64,70,3)

[1] 0.02275013
  1. La probabilidad de que pesen 64 kg.

\(P(X = 64)\)=0 , ya que cuando la distribución es continua, la probabilidad de que la variable tenga un valor exacto es siempre es nula.

5 Contraste de hipótesis - paramétrico

Una empresa de neumáticos afirma que una nueva gama en promedio dura al menos 28.000 km. Las pruebas con 64 neumáticos dan como resultado una duración media de 27.800 km, con una desviación estándar de 1.000 km.

  1. ¿Cuál es la hipótesis nula y cuál es la hipótesis alternativa?

    Como bien sabemos la hipótesis nula (H₀) esla afirmación que quiere ser constrastada y la hipótesis alternativa (H₁) es la opuesta a la hipótesis nula. Por tanto:

    H₀: La media de la duración de los neumáticos es igual o más de 28.000 km.

    H₁: La media de la duración de los neumáticos es menor de 28.000 km.

  2. Determina la región critica, con un nivel de significación del 5%.

\(Z_{\alpha/2}\)={r} qnorm(1-0.05/2) , es el valor que deja en las colas el % de la distribución normal N(0,1).

[1] 1.959964
  1. Comprobar si hay evidencia suficiente para rechazar la afirmación de la empresa.

\(Z_{s}=(27800-28000)/(1000/sqrt{64})\)=1.6

Como \(Z_{s}>-Z_{\alpha/2}\) , entonces no se rechaza la hipótesis nula. Por tanto, la empresa puede afirmar que una nueva gama en promedio dura al menos 28.000 km. 

 ```{r}  n=64
         s=1000
         z=(abs(27800-28000))/(s/sqrt(n))
         z 
 ```
  1. ¿Cuál es el valor de p-valor?

El valor p-valor es el menor nivel de significación para el que se rechaza la hipótesis nula: 

``{r}   pvalor=2*pnorm(z) 
        pvalor
`` [1] 1.890401

6 Contraste no paramétrico

Se han recopilado datos sobre el tiempo (en segundos) que tarda cada empresa en atender una llamada durante un período de un mes. Los datos son los siguientes:

Empresa A: 20, 25, 30, 22, 28, 35, 27, 26, 32, 29 Empresa B: 18, 21, 28, 20, 24, 29, 23, 22, 27, 26

  1. Haz uso de la prueba de los signos con nivel de satisfaccion al 0.07 para probar que la hipotesis nula de que la empresa A no es mas rapida que la B

``{r}
  # Datos de tiempos de espera de las dos empresas
tiempo_empresa_A <- c(20, 25, 30, 22, 28, 35, 27, 26, 32, 29)
tiempo_empresa_B <- c(18, 21, 28, 20, 24, 29, 23, 22, 27, 26)

# Calcular la diferencia entre los tiempos de espera de las dos empresas
diferencia <- tiempo_empresa_A - tiempo_empresa_B

# Contar el número de veces que la empresa A es más rápida que la empresa B
num_positivos <- sum(diferencia > 0)

# Nivel de significación
nivel_significacion <- 0.07

# Realizar la prueba de los signos
p_valor <- binom.test(num_positivos, num_positivos + sum(diferencia < 0), p = 0.5, alternative = "greater")$p.value

# Comprobar si el p-valor es menor que el nivel de significación
if (p_valor < nivel_significacion) {
  mensaje <- "Se rechaza la hipótesis nula. La empresa A es más rápida que la empresa B en atender llamadas."
} else {
  mensaje <- "No se rechaza la hipótesis nula. No hay suficiente evidencia para afirmar que la empresa A es más rápida que la empresa B en atender llamadas."
}

  1. ¿y si queremos comprobar que son iguales?

``{r}
 # Calcular la diferencia entre los tiempos de respuesta de las dos empresas
diferencia <- tiempo_empresa_A - tiempo_empresa_B

# Contar el número de veces que la diferencia es igual a cero
num_ceros <- sum(diferencia == 0)

# Nivel de significación
nivel_significacion <- 0.07

# Realizar la prueba de los signos
p_valor <- binom.test(num_ceros, length(diferencia), p = 0.5, alternative = "two.sided")$p.value

# Comprobar si el p-valor es menor que el nivel de significación
if (p_valor < nivel_significacion) {
  mensaje <- "Se rechaza la hipótesis nula. Los tiempos de respuesta de las empresas A y B son diferentes."
} else {
  mensaje <- "No se rechaza la hipótesis nula. No hay suficiente evidencia para afirmar que los tiempos de respuesta de las empresas A y B son diferentes."
}

  1. Comprueba la hipotesis principal con el test de Wilcoxon y los siguientes valores con una confianza al 91%

``{r}
# Datos de tiempos de espera de las dos empresas
tiempo_empresa_A <- c(20, 25, 30, 22, 28, 35, 27, 26, 32, 29)
tiempo_empresa_B <- c(18, 21, 28, 20, 24, 29, 23, 22, 27, 26)

# Realizar el test de Wilcoxon
resultado_wilcoxon <- wilcox.test(tiempo_empresa_A, tiempo_empresa_B, alternative = "two.sided", conf.level = 0.91)

# Imprimir resultados
print(resultado_wilcoxon)

7 Contraste de hipótesis

Un grupo de automóviles asegura que la velocidad de sus vehículos es en promedio unos 120 km/h. Las purebas con 50 coches dan como resultado una media de 115 km/h, con una desviación estándar de 5 km/h.

  1. Calcula el valor del estadístico de contraste.
#Valor del estadístico de contraste
n=50
s=2*sqrt(50/5)
z=(abs(115-120)) / (s/sqrt(n))
z
## [1] 5.59017

  1. ¿Cuál es el valor del p-valor?

    El valor p-valor es el menor nivel de significación para el que se rechaza la hipótesis nula:

#P-valor
pvalor=2*pnorm(z) 
pvalor
## [1] 2

  1. Supongamos que queremos saber si hay una diferencia significativa en la velocidad promedio entre este grupo de automóviles anterior y otro grupo. Tenemos muestras independientes de la velocidad promedio de cada grupo:

    Grupo 1: [100, 110, 120, 105, 115] (km/h) Grupo 2: [95, 105, 115, 100, 110] (km/h)

8 Anova

La uclm ha realizado distintas mediciones a varias marcas de automóviles, sobre la potencia del motor(CV) de sus modelos, que se expresan en la siguiente tabla:

Audi: 150, 226, 333, 180, 272, 450, 140

BMW: 500, 186, 235, 290, 380, 155, 420

Mercedes: 233, 160, 265, 475, 141, 320, 190

Volkswagen: 116, 202, 400, 321, 150, 256, 175

  1. Calcular el valor estadístico:

  2. Calcula el p-valor:

``{r}
Audi = c(150, 226, 333, 180, 272, 450, 140)
BMW = c(500, 186, 235, 290, 380, 155, 420)
Mercedes = c(233, 160, 265, 475, 141, 320, 190)
Volkswagen = c(116, 202, 400, 321, 150, 256, 175)
datos = c(Audi, BMW, Mercedes, Volkswagen)

marcas = factor(rep(c("Audi", "BMW", "Mercedes", "Volkswagen"), c(7,7,7,7)))
summary(aov(datos~marcas))
##             Df Sum Sq Mean Sq F value Pr(>F)
## marcas       3  23622    7874   0.602   0.62
## Residuals   24 313975   13082

Como podemos observar, tras compilar el fragmento de arriba obtenemos las soluciones de los enunciados. Con el summary obtenemos un conjunto de datos, entre los cuales se encuentran el valor estadístico, que en este caso es 0.602, y el p-valor: 0.62.

  1. Calcula el límite superior de confianza de la potencia de los automóviles de BMW:
mediaBMW <- mean(BMW)
valor_t <- pt(0.05/2, 28 - 1)
sp <-sqrt(13082)
errorEstimacion <- valor_t * (sp/ sqrt(7))
mediaBMW + errorEstimacion
## [1] 331.4709
  1. Crea un diagrama de cajas con los valores dados:
boxplot(datos~marcas, col=c("yellow","blue","red","green"), ylab="Potencia de los automóviles")

9 Regresión lineal

El enlace a Rpubs: https://rpubs.com/DanielBalan/1178087

Enlaces videos:

Primeros dos ejercicios: https://files.fm/u/d86b3nm7v4?share_email_id=13f184d&share_show_reg_form=true

Tercer Ejercicio: https://files.fm/u/u2ev4b8yns

Video Trabajo 2: https://www.flexclip.com/es/share/5650952c118ebffb4afb53384b64392a9041440.html

Video Trabajo 3: