Updated beta.R to vectorize arguments, still need to find a block design