How to fit multiple independent slopes and intercepts with the GLM module?

If you coded the linear equation with interaction between group and slope, it should gives you what you want: y = x1 + x2 + x1:x2