> dbinom(0, 10, 0.3)
[1] 0.02824752
> dbinom(1, 10, 0.3)
[1] 0.1210608
> pbinom(1, 10, 0.3)
[1] 0.1493083
> dbinom(2, 10, 0.3)
[1] 0.2334744
> pbinom(2, 10, 0.3)
[1] 0.3827828
> #Testiküsimusel on 4 varianti, neist 1 õige
> #Testis on 5 küsimust
> #Millise tõenäosusega saab teemaga tundmatu inimene 3 või rohkem õiget vastust
> dbinom(3, 5, 0.25)
[1] 0.08789063
> dbinom(3:5, 5, 0.25)
[1] 0.0878906250 0.0146484375 0.0009765625
> sum(dbinom(3:5, 5, 0.25))
[1] 0.1035156
> pbinom(2, 5, 0.25)
[1] 0.8964844
> 1-pbinom(2, 5, 0.25)
[1] 0.1035156
> 

toenaosus=0.3
seeriapikkus=5
uuritavpihtas=3
choose(seeriapikkus, uuritavpihtas)*toenaosus**uuritavpihtas*(1-toenaosus)**(seeriapikkus-uuritavpihtas)
dbinom(uuritavpihtas, seeriapikkus, toenaosus)

binoomtoenaosus <- function(k, n, p){
   return (factorial(n)/(factorial(k)*factorial(n-k))*(p**k)*((1-p)**(n-k)))
}

binoomtoenaosus(3, 5, 0.3) 



> x=0:5
> y=sapply(x, function(x) dbinom(x, 5, 0.3))
> plot(x, y)
> boxplot(x, y)
> barplot(x, y)
> barplot(y)
> help(barplot)
> barplot(y, legend.text=x)
> barplot(y, names.arg=x)
> barplot(y, names.arg=x, col="red")
> barplot(y, names.arg=x, col="red", main="binoomjaotus")
> barplot(y, names.arg=x, col="red", main="binoomjaotus", ylab="osakaal", xlab="tabamuste arv")
> barplot(y, names.arg=x, col=c("red", "blue"), main="binoomjaotus", ylab="osakaal", xlab="tabamuste arv")
> #Koostage funktsioon, mille parameetriteks on binoomjaotuse pikkus ja katse õnnestumise 
> # #tõenäosus ning mis joonistab vastava binoomjaotuse kohta tulpdiagrammi

binoomjoonis <- function(n, p){
   x=0:n
   y=sapply(x, function(x) dbinom(x, n, p))
   barplot(y, names.arg=x, ylab="osakaal")
}

binoomjoonis(10, 0.4)

> par(mfrow=c(1, 2))
> binoomjoonis(5, 0.2)
> binoomjoonis(5, 0.6)
> par(mfrow=c(1, 1))

joonisedFaili <- function(n, failinimi){
   pdf(failinimi, width=10, height=10)
   par(mfrow=c(2, 2))
   binoomjoonis(n, 0.2)
   binoomjoonis(n, 0.4)
   binoomjoonis(n, 0.6)
   binoomjoonis(n, 0.8)
   par(mfrow=c(1, 1))
   dev.off()
}


joonisedFaili(5, "c:\\jaagup\\viiene.pdf")