LTLA

# Demonstrate that this is chi-squared distributed
# with xyz - (x-1) - (y-1) - (z-1) -1 d.f.
stat <- numeric(1000)
for (i in seq_along(stat)) {
    N <- 1000
    x <- rmultinom(1, size=N, prob=rep(1/27, 27))
    dim(x) <- c(3,3,3)

    p1 <- vapply(1:3, function(i) sum(x[i,,]), 1) / N
    p2 <- vapply(1:3, function(i) sum(x[,i,]), 1) / N

LTLA

Some thoughts on the OSCA books

Organization

The "OSCA" concept should be more of a federated collection of books instead of a centrally managed resource. This will enable greater contributions from the Bioconductor community, rather than being restricted to a handful of maintainers. Most importantly, books can be contributed by the same developers who maintain the package(s) used by the book. This avoids a major problem in book maintenance whereby package updates cause unexpected changes to the results in the book.

The flipside of many separate books is that a greater burden is placed on the reader to assemble the disparate bits of knowledge into a functional workflow.