Probabilitati laborator repartitii

#Uniforma pe multimea (1,2,3,4)
x = sample(c(1,2,3,4), 100, replace = TRUE,  prob = c(0.25, 0.25 0.25, 0.25))
#Bernoulli de parametru p
p <- 0.3
x1 = sample(c(0,1), 10, replace = TRUE,  prob = c(1-p,p))

#binomiala
s <- sum(x1)

#repartitia binomiala de parametrii n si p
p <- 0.6
n <- 20

par(mfrow = c(1,2)) #pentru a afisa 2 grafice simultan
plot(0:20,dbinom(0:20,n,p),col="magenta",main = "Functia de masa a binomialei",xlab = "nr. de reusite", ylab="probabilitate")

#functia de repartitie
x2 <- seq(-5,30,0.001)
#(0:20,pbinom(0:20,n,p),col="green",main="Functia de repartitie a binomialei")

#plot(0:20,pbinom(0:20,n,p),col="green",main="Functia de repartitie a binomialei")

plot(x2,pbinom(x2,n,p),col="green",main="Functia de repartitie a binomialei",lwd=0.5,lty=1)
y <- rbinom(10,n,p)