Agresti_create_design_matrix function

Creates the design matrix for Agresti's simple diagonal quasi-symmetry model.