Steekproef uit een Poisson-verdeling
We hebben het gemiddelde jaarlijkse aantal mishandelingen, berovingen en slagen in een kleine stad gekregen (respectievelijk 150, 300 en 50). Aangezien de drie typen misdrijven gemodelleerd kunnen worden als een multivariate Poisson-verdeling, maak je 10 verschillende scenario’s voor het aantal misdrijven in het komende jaar.
Deze oefening maakt deel uit van de cursus
Mixture Models in R
Oefeninstructies
- Gebruik de functie
rpoisom 10 observaties te genereren voor elk type misdrijf met de bijbehorende lambda. - Combineer het gesimuleerde aantal misdrijven tot een multivariate Poisson-verdeling.
- Controleer de waarden van je gesimuleerde scenario’s.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
set.seed(1541)
# Create the vector of lambdas
lambda_1 <- c(150, 300, 50)
# Create the sample of each crime
assault_1 <- ___(n = ___, lambda = ___)
robbery_1 <- ___(n = ___, lambda = ___)
battery_1 <- ___(n = ___, lambda = ___)
# Combine the results
cities_1 <- cbind(assault_1, ___, battery_1)
# Check the sample
cities_1