Eksempel på Binomiske Sannsynligheter

require(pacman)
p_load(mosaic)

Beregne Binomiske Sannsynligheter (fra Kristensen og Wikan), eksempel s. 117. For X~B(4, 0.6).

Funksjonen i R er: dbinom(x, size, prob)

2. Finn P(X=0):

dbinom(0, 4, 0.6) 

[1] 0.0256

For alle 5 verdier i utfallsrommet, samtidig:

dbinom(0:4, 4, 0.6)   

[1] 0.0256 0.1536 0.3456 0.3456 0.1296

Et plot av denne fordelingen:

plotDist("binom", size=4, prob=0.6)

Et plot av den kumulative fordelingen:
plotDist("binom", size=4, prob=0.6, kind='cdf')

Gjennomsnitt, formel 3.3, s. 108

sum(0:4*dbinom(0:4, 4, 0.6))

[1] 2.4

Lagre gjennomsnitt i m

m <- sum(0:4*dbinom(0:4, 4, 0.6))

Varians, formel 3.5, s. 109

sum((0:4-m)^2*dbinom(0:4, 4, 0.6))

[1] 0.96

St.dev

sqrt(sum((0:4-m)^2*dbinom(0:4, 4, 0.6)))

[1] 0.9797959