plot(
T2 <- table(randI(1e6, 2))) ; abline(h = 1e6 / (9*10^(2 - 1)), lty=2, col="gray70")
chisq.test(T2) # typically not at all significant
T3 <- table(randI(1e6, 3))
chisq.test(T3)
stopifnot(exprs = {
identical( 10:99 , as.integer(names(T2)))
identical(100:999, as.integer(names(T3)))
})
Run the code above in your browser using DataLab