Code
files <- c("M1", "M2", "G1", "G2", "G3")
pairs <- combn(files, 2)
bad <- colSums(matrix(startsWith(pairs, "M"), nrow = 2))
data.frame(
formula = c("choose(5,2)", "choose(2,1)*choose(3,1)", "P(exactly one bad)"),
counted = c(choose(5, 2), choose(2, 1) * choose(3, 1), 0.6),
enumerated = c(ncol(pairs), sum(bad == 1), mean(bad == 1))
) formula counted enumerated
1 choose(5,2) 10.0 10.0
2 choose(2,1)*choose(3,1) 6.0 6.0
3 P(exactly one bad) 0.6 0.6
