I\'m using Mathematica for a project, in which I would like to fit a multiple linear regression model to my data (which is made up of both numerical and categorical variable