Poisson 분포에서 샘플 추출하기
한 소도시에서 연간 평균 폭행, 강도, 폭력 사건 건수가 각각 150, 300, 50건으로 주어졌어요. 이 세 가지 범죄 유형이 다변량 Poisson 분포로 모형화된다고 가정할 때, 내년도 범죄 발생 건수에 대한 10개의 서로 다른 시나리오를 생성하세요.
이 연습은 강의의 일부입니다
R로 배우는 Mixture Models
연습 안내
rpois함수를 사용해 각 범죄 유형에 대해 해당 람다를 이용한 관측값 10개를 생성하세요.- 시뮬레이션한 범죄 건수를 결합해 다변량 Poisson 분포 형태로 구성하세요.
- 생성한 시나리오의 값을 확인하세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
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