scratch - HenrikBengtsson/matrixStats GitHub Wiki
mu <- matrixStats::colMedians(X, rows=33:158, cols=1001:3000)
is much faster and uses way less memory than
mu <- apply(X[33:158, 1001:3000], MARGIN=2, FUN=median)
mu <- matrixStats::colMedians(X, rows=33:158, cols=1001:3000)
is much faster and uses way less memory than
mu <- apply(X[33:158, 1001:3000], MARGIN=2, FUN=median)