Objetivo. Determinar la probabilidad binomila de lectores de novela

Descripción

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. Un grupo de 4 amigos son aficionados a la lectura:

1 ¿Cuál es la probabilidad de que en el grupo, hayan leido la novela 2 personas?

2 ¿y Cómo máximo 2?. Acumulado hasta 2, incluye: 0, 1 y 2

Valores de los parámetros

n <-  4
p <- 0.80 # EXITO
q <- (1 - p) # 0.2.  # FRACASO
x <-  2 # Variable discdreta a determinar su probabilidad

1 ¿Cuál es la probabilidad de que en el grupo, hayan leido la novela 2 personas?

x = 2

dbinom(x,n, p)
## [1] 0.1536

2 ¿y cómo máximo 2?. Acumulado hasta 2, incluye: 0, 1 y 2

# sumar la probabilidad de 0 mas la prob. de 1 mas la prob de 
# x menor o igual a 2

x <- 2; n <- 4; p <- 0.80  # todo en una misma linea
dbinom(0,n, p) + dbinom(1,n, p) + dbinom(2,n, p)
## [1] 0.1808
# Es lo mismo que si utilizamos pbinom
pbinom(x, n , p)  # Esta es probabilida dacumulada hastya el valor de x
## [1] 0.1808

Plot la probabilidades desde 0 hasta 4

plot (x=c(0:4), y=dbinom(0:4, n, p), type="b", col=terrain.colors(5), pch=20, xlab = 'Variable discreta x de 0 a 4', ylab = 'Probabilidad p(x)' )

legend("topleft", inset=.03, title="Binomial", 
       as.character(paste("p(x) en ", 0:4, " = ",round(dbinom(0:4, n, p),4))), fill=terrain.colors(4), horiz=FALSE)