yXn function

Prepare design matrix for nested model with groups and subgroups