modelMatrices function

Extract SEM model matrices