Simulation of a single n
With a room of size 50, let us start by finding a solution via simulation.
In this code, n
will represent our room size and match
will be our counter for the number of times at least one matching birthday occurred in the simulation, which starts at 0 and should increment in each iteration where a match occurs.
This exercise is part of the course
Probability Puzzles in R
Exercise instructions
- Fill in the required parameters of the
sample
function, to randomly assign birthdays to each individual. - Increment the
match
variable according to the condition it follows. - Calculate the estimated probability of at least one match.
Hands-on interactive exercise
Have a go at this exercise by completing this sample code.
set.seed(1)
n <- 50
match <- 0
# Simulate 10000 rooms and check for matches in each room
for(i in 1:10000){
birthdays <- sample(___, ___, replace = ___)
if(length(unique(birthdays)) < n){
match <- ___
}
}
# Calculate the estimated probability of a match and print it
p_match <- ___
print(___)