Campionare da una distribuzione di Poisson
Ti sono stati forniti i numeri medi annui di aggressioni, rapine e percosse in una piccola città (150, 300 e 50, rispettivamente). Dato che i tre tipi di reato possono essere modellati con una distribuzione di Poisson multivariata, ti viene chiesto di creare 10 scenari diversi per il numero di reati nel prossimo anno.
Questo esercizio fa parte del corso
Modelli di Mixture in R
Istruzioni dell'esercizio
- Usando la funzione
rpois, crea 10 osservazioni per ciascun tipo di reato utilizzando il rispettivo lambda. - Combina il numero simulato di reati in una distribuzione di Poisson multivariata.
- Controlla i valori dei tuoi scenari simulati.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
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