model.matrix.pim function

Create a model matrix for a probabilistic index model