Too minor for a pull request <g> `X_mean_0 <- sweep(x, 2, colMeas(x))` should be `X_mean_0 <- sweep(x, 2, colMeans(x))`