La última novela de un autor ha tenido un gran éxito, hasta el punto de que el 80% de los lectores ya la han leido. Hallar la probabilidad de que en un grupo de cuatro amigos que son aficionados a la lectura, dos hayan leido la novela.
La probabilidad de que una persona haya leido el libro es de \(p=0.8\), por lo que la probabilidad de que no lo haya leido es de \(q=0.2\)
n <- 4
prob <- 0.80
\(p(x;n;p) = \binom nxp^xq^(n-x); x=0,1,2...n\)
\(\binom nx=\frac{n!}{x!(n-x)!}\)
Para cuando \(x=0\) Se determina la función o probabilidad para cuando \(x=0\) confome a la fórmula y luego conforme a la funcion de R: dbinom(). El resultado es el mismo
x <- 0
n <- 4
p <- prob
q = 1 - p
(factorial(n) / (factorial(x) * factorial(n-x))) * p^x * (1-p)^(n-x)
## [1] 0.0016
#o bien
dbinom(x = 0, size = n, prob = prob)
## [1] 0.0016
x <- 2
dbinom(x,n,prob)
## [1] 0.1536
x<-2
pbinom(x,n,prob)
## [1] 0.1808
x<-1
1 - pbinom(x,n,prob)
## [1] 0.9728
pbinom(3,n,prob) - pbinom(1,n,prob)
## [1] 0.5632
tabla <- data.frame(c(0:4),dbinom(0:4,n,prob), pbinom(0:4,n,prob))
colnames(tabla) <- c("x", "prob.x", "prob.acum.x")
tabla
## x prob.x prob.acum.x
## 1 0 0.0016 0.0016
## 2 1 0.0256 0.0272
## 3 2 0.1536 0.1808
## 4 3 0.4096 0.5904
## 5 4 0.4096 1.0000
barplot(height = tabla$prob.x, names.arg = tabla$x,
xlab = "Valores de x",
ylab = "Probabilidades")
plot(tabla$x, tabla$prob.acum.x, type = "b",
xlab = "Valores de x",
ylab = "Probabilidad acumulada")
\(\mu = np\)
v.e <- n * prob
v.e
## [1] 3.2
\(\sigma^2 = npq\)
vari <- n * p * q
vari
## [1] 0.64
\(\sqrt{\sigma^2}\)
desv.std <- sqrt(vari)
desv.std
## [1] 0.8