get_design_matrix_pglmm function

get_design_matrix gets design matrix for gaussian, binomial, and poisson models