Ta stickprov från Poissonfördelningen
Vi har fått information om det genomsnittliga årliga antalet överfall, rån och misshandel i en liten stad (150, 300 respektive 50). Eftersom de tre brottstyperna kan modelleras som en multivariat Poissonfördelning ska du skapa 10 olika scenarier för antalet brott under nästa år.
Den här övningen är en del av kursen
Blandningsmodeller i R
Övningsinstruktioner
- Använd funktionen
rpoisför att skapa 10 observationer för varje brottstyp med dess motsvarande lambda. - Kombinera de simulerade brottsantalen till en multivariat Poissonfördelning.
- Kontrollera värdena i dina simulerade scenarier.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
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