Predicted values now calculated across all covariate combos w/ matrix algebra
merged experimental matrix code into main branch