You could use either the regress() function or the LinearModel.fit() function.
doc regress
or
doc LinearModel.fit
for details.
In each case, you'll need to combine your X and Y values into one array (predictor variables) and use Z for the response variable.
