sample_space <- 1:6 sample_space

event_A <- c(2, 4, 6) event_B <- c(5, 6)

event_A event_B

intersect(event_A, event_B)

P_A <- length(event_A) / length(sample_space) P_B <- length(event_B) / length(sample_space)

P_A

P_A_given_B <- length(intersect(event_A, event_B)) / length(event_B) P_A_given_B

values <- 1:6 probabilities <- rep(1/6, 6)

distribution <- data.frame(X = values, P_X = probabilities) distribution

expected_value <- sum(values * probabilities) expected_value

variance <- sum((values - expected_value)^2 * probabilities) standard_deviation <- sqrt(variance)

variance standard_deviation P_B