log_linear_quasi_symmetry_model_design function

Creates the design matrix for a quasi-symmetry design