zaaltje

inputs:

populatiegrootte = pop

besmettingsgraad = b

zaalgrootte = z

aantal trekkingen van een nieuw zaaltje = n

output:

frequentie waarbij een zaaltje minstens 1 besmet persoon bevatte

pop<- rep(0, 17000000)             # zeg, populatie is 100000000
b <-0.0001                          # stel besmettingsgraag 1:10.000
a<- length(pop)*b 
pop[1:a]<-1                         # de eerste maak ik 'besmet'

sum(pop)                            # absoluut aantal besmette personen

## [1] 1700

length(pop)                        # absoluut populatie grootte

## [1] 17000000

z <- 1000
zaaltje <-vector("numeric", z)    # maak de lege vector aan
n=1000                              # trek 1000 keer een zaaltje van 100 uit pop
for (i in 1:n)                      
  zaaltje[i] <- sum(sample(pop,100))  #besmet kijkt of er >1 'besmet persoon zit'
ja <- sum(zaaltje)
prob <- ja/1000
###########################################################################

Plots

plot(zaaltje)