masa = dbinom(3, size = 3,prob = 0.8)
masa
## [1] 0.512
library("Rlab")
## Rlab 2.15.1 attached.
##
## Attaching package: 'Rlab'
## The following objects are masked from 'package:stats':
##
## dexp, dgamma, dweibull, pexp, pgamma, pweibull, qexp, qgamma,
## qweibull, rexp, rgamma, rweibull
## The following object is masked from 'package:datasets':
##
## precip
rango = seq(0,3)
distribucion = dbinom(rango, size = 3,prob = 0.8)
datos=data.frame(rango,distribucion)
#Gráfico
library("ggplot2")
grafico = ggplot(data=datos,aes(x=rango,y=distribucion))
grafico = grafico + geom_bar(stat="identity",fill="black")
grafico = grafico + theme_bw() + ggtitle("probabilidad de telefonos")
grafico = grafico + xlab("telefonos") + ylab("Probabilidad")
plot(grafico)
print((1-dnbinom(x=2, size= 2, prob=0.1)))
## [1] 0.9757
#valor esperado en esta distrución es r/p
valor_esperado = 2/0.1
valor_esperado
## [1] 20
personas = seq(0,10)
positivos=2
distribucion = dnbinom(x=personas, size=positivos, prob=0.1)
datos=data.frame(personas,distribucion)
#Gráfico
library("ggplot2")
grafico = ggplot(data=datos,aes(x=personas,y=distribucion))
grafico = grafico + geom_bar(stat="identity",fill="red")
grafico = grafico + theme_bw() + ggtitle("Probabilidad para 2 personas")
grafico = grafico + xlab("personas") + ylab("Probabilidad")
plot(grafico)
print(dhyper(x=1, m=240, k=10, n=560))
## [1] 0.1200794
print(1-dhyper(x=1, m=240, k=10, n=560))
## [1] 0.8799206
library("Rlab")
poblacion=240 #Exitos
total=560 #Fracasos
hombres=10 #Experimentos
exitos=seq(0,10) #Exitos
distribucion = dhyper(x=exitos, m=poblacion, k=hombres, n=total)
datos=data.frame(exitos,distribucion)
#Gráfico
library("ggplot2")
grafico = ggplot(data=datos,aes(x=exitos,y=distribucion))
grafico = grafico + geom_bar(stat="identity",fill="green")
grafico = grafico + theme_bw() + ggtitle("Distribución de
probabilidades")
grafico = grafico + xlab("hombres") + ylab("Probabilidad")
plot(grafico)
lambda= 8
print(dpois(5,lambda))
## [1] 0.09160366
print(dpois(3,lambda) + dpois(2,lambda) + dpois(1,lambda))
## [1] 0.04204465
media= 3 * 0.8
desviacion= (3 * 0.8 * 0.2)**(1/2)
aproximacion= pnorm((3.5 - media)/desviacion)
aproximacion
## [1] 0.9438244
k = 5
lambda= 8
X = (k-lambda)/(sqrt(lambda))
aproximacion= pnorm(X, 0, 1)
aproximacion
## [1] 0.1444222
k = (seq (0,3))
lambda= 8
X = (1-lambda)/(sqrt(lambda))
X2 = (2-lambda)/(sqrt(lambda))
X3 = (3-lambda)/(sqrt(lambda))
aproximacion= pnorm(X, 0, 1) + pnorm(X2, 0, 1) + pnorm(X3, 0, 1)
aproximacion
## [1] 0.06216153